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PREFACE 


This  volume  ,  is  the  second  of  six  volumes.,  that  describe  the  LMI  Sortie- 
Generation  Model  System.  Volume  I,  Executive  Summary,  discusses  the  problem 
the  system  is  designed  to  address  and  provides  an  overview  of  the  principal 
parts  of  the  system.  Volume  II,  Sortie-Generation  Model  User's  Guide,  pro¬ 
vides  sufficient  information  to  allow  a  user  to  run  the  Sortie-Generation 
Model  (SGM) .  Volume  III,  Sortie-Generation  Model  Analyst's  Manual,  describes 
the  mathematical  structures,  derivations,  assumptions,  limitations,  and  data 
sources  of  the  SGM  at  a  very  detailed  level.  Volume  IV,  Sortie-Generation 
Model  Programmer's  Manual,  specifies  the  details  of  the  computer  programs, 
file  structures,  job  control  language,  and  operating  environment  of  the  SGM. 
Volume  V  describes  the  maintenance  subsystem  and  explains  the  construction  of 
the  maintenance  input  file  to  the  SGM.  Volume  VI  describes  the  spares  sub¬ 
system  and  shows  a  user  how  to  build  the  spares  file  that  is  used  by  the  SGM. 

Potential  users  are  cautioned  that  no  volume  is  intended  to  provide,  by 
itself,  all  of  the  information  needed  for  a  comprehensive  understanding  of  the 
operation  of  the  SGM. 
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1 .  MODEL  DESCRIPTION 


INTRODUCTION 

The  Sortie-Generation  Model  (SGM)  is  a  hybrid  analytic/simulation  model 
that  estimates  the  expected  maximal  number  of  sorties  that  can  be  flown  by  a 
specified  aircraft  type  in  a  wartime  scenario.  This  estimate  is  based  on 
aircraft  characteristics,  maintenance  manpower  and  recoverable  spares  levels, 
and  user  inputs  that  describe  the  scenario  of  interest.  The  purpose  of  this 
chapter  is  to  describe  the  basic  structure  of  the  SGM.  The  reader  who  is 
interested  in  the  technical  aspects  of  the  SGM  is  referred  to  Volume  III. 
THE  STRUCTURE  OF  THE  MODEL:  STATES  AND  PROCESSES 

The  SGM  consists  of  a  collection  of  aircraft  states ,  processes  that  cause 
transitions  between  states,  and  logic  that  governs  those  processes.  The  SGM 
simulates  the  transition  of  aircraft  between  these  states  throughout  a  daily 
flying  schedule  that  is  specified  by  the  user.  The  definitions  of  the  states, 
the  logic  of  the  state  transitions,  and  the  interaction  of  these  transitions 
with  the  flying  schedule  determine  the  basic  structure  of  the  SGM. 

Aircraft  States 

There  are  five  aircraft  states  in  the  SGM: 

1)  Mission-capable 

2)  Maintenance 

3)  Not  mission-capable,  supply  (NMCS) 

4)  Combat  loss 

5)  Reserve 

These  states  are  mutually  exclusive  and  collectively  exhaustive;  i.e.,  every 
aircraft  is  in  one  and  only  one  state.  The  states  are  described  below. 
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Mission  Capable.  An  aircraft  is  considered  mission-capable  if  it  is 
capable  of  flying  a  combat  mission.  It  is  not  mission-capable  if  it  is  under¬ 
going  essential  corrective  maintenance,  or  is  missing  a  mission-essential 
part.  There  is  no  explicit  representation  in  the  SGM  of  aircraft  that  are 
partially  mission-capable. 

Maintenance .  An  aircraft  is  in  maintenance  status  if  it  requires 
unscheduled,  on-aircraft  repair  that  is  essential  to  the  performance  of  its 
mission.  This  repair  may  or  may  not  be  due  to  failure  of  a  part;  however,  in 
this  model,  an  aircraft  is  not  allowed  to  enter  maintenance  until  all  needed 
parts  have  been  obtained  from  supply  or  repair. 

NMCS .  An  aircraft  is  not  mission-capable,  supply  if  the  aircraft  is 
missing  an  essential  part.  In  the  SGM,  only  mission-essential  Line  Replace¬ 
able  Units  (LRUs)  can  cause  an  aircraft  to  become  NMCS. 

Combat  Loss.  An  aircraft  is  counted  as  a  combat  loss  if  it  does  not 
return  from  a  sortie.  Once  an  aircraft  has  been  lost  it  can  never  be  re¬ 
covered.  Battle-damaged  aircraft  that  return  from  a  sortie  are  not  considered 
in  this  model. 

Reserve.  Reserve  aircraft  consist  of  mission-capable  aircraft  that 
are  used  to  replace  combat  losses.  The  user  specifies  an  initial  number  of 
aircraft  that  are  held  in  reserve  at  the  beginning  of  the  scenario;  these 
reserve  aircraft  replace  combat  losses  at  the  end  of  each  day,  until  all 
reserves  have  been  exhausted.  Aircraft  are  allowed  to  leave  this  reserve 
state,  but  no  aircraft  can  enter  it. 

Processes  -  Transitions  Between  States 

There  are  eight  processes  in  the  SGM  which  cause  transitions  between 
aircraft  states: 

1)  Ground  aborts 
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2)  Breaks 


3)  Aircraft  repairs 

4)  Parts  demands 

5)  Parts  repair 

6)  Cannibalization 

7)  Attrition 

8)  Commitment  of  reserves 

Figure  1-1  depicts  the  relationships  among  the  various  states  and  processes. 
THE  LOGIC  OF  THE  MODEL:  EVENTS  AND  THEIR  SEQUENCING 

The  events  that  occur  in  the  SGM  are  related  to  a  flying  schedule  with 
user-specified  characteristics.  The  flying  schedule  consists  of  a  number  of 
periods  or  cycles  each  of  which  is  divided  into  three  segments  of  lengths  T^, 
Tp,  and  Ty,  respectively.  During  the  last  period,  the  segment  is  replaced 
by  an  overnight  recovery  period.  The  user  specifies  the  first  and  last  take¬ 
off  times  of  the  day;  the  time,  T^,  which  is  the  average  minimal  length  of 
time  required  to  launch  a  sortie  given  a  mission-capable  aircraft;  the  sortie 
length,  Tpj  and  the  number  of  periods  per  day.  The  time,  Tw>  is  then  computed 
by  the  SGM  program.  The  flying  schedule  is  the  same  each  day  except  for  the 
number  of  aircraft  to  be  flown  each  period,  which  the  user  can  vary.  A 
typical  flying  schedule  is  portrayed  in  Figure  1-2. 

Figure  1-3  portrays  two  segments  of  a  flying  day;  the  flying  period  on 
the  left  is  intended  to  be  typical  and  the  one  on  the  right  to  be  the  last 
period  of  the  day.  The  events  that  occur  in  the  SGM  are  denoted  by  circled 
numbers  placed  under  the  figures  at  the  appropriate  positions  on  the  time 
line.  Each  of  those  events  is  described  here. 
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FIGURE  1-1 

STATES  AND  PROCESSES 


Event  1 


All  mission-capable  aircraft  are  prepared  for  launch.  Any  aircraft 
that  is  not  mission-capable  at  this  time  (i.e.,  before  takeoff)  cannot  be 
flown  during  this  cycle  because,  by  definition,  is  the  minimal  time  re¬ 
quired  to  launch  an  aircraft  that  is  mission-capable. 

Event  2 

Aircraft  that  are  repaired  during  the  period  of  length  leave 
maintenance  and  become  mission-capable  but  are  not  available  for  flight  during 
this  cycle. 

Event  3 

All  aircraft  that  were  prepared  for  takeoff  are  subjected  to  the 
probability  of  ground  abort.  A  ground  abort  is  defined  as  an  unsuccessful 
attempt  by  an  aircrew  to  fly  an  aircraft.  The  aborted  aircraft  enter  mainte¬ 
nance.  No  parts  demands  are  generated  by  ground  aborts. 

Event  4 

The  remaining  aircraft  that  were  prepared  for  takeoff  fly  sorties. 

Event  5 

Each  aircraft  that  flies  is  subjected  to  the  probability  of  attri¬ 
tion  and,  for  each  combat  loss,  an  aircraft  is  deducted  from  the  current 
strength  of  the  organization. 

Event  6 

Aircraft  that  are  repaired  during  the  period  of  length  T_  leave 

r 

maintenance  and  become  mission-capable. 

Event  7 

Each  aircraft  returning  from  flight  is  subjected  to  the  probability 
of  break,  i.e.,  the  probability  of  requiring  essential  corrective  maintenance 
prior  to  flying  another  combat  mission.  At  the  same  time,  parts  demands  are 
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generated.  Demands  that  can  be  filled  from  stock  on-hand  result  in  issues  of 
that  stock.  Demands  that  cannot  be  filled  from  stock  and  cannot  be  satisfied 


by  cannibalization  from  aircraft  that  are  NMCS  result  in  additional  aircraft 
becoming  NMCS. 

Event  8 

Aircraft  that  are  repaired  during  the  period  of  length  leave 
maintenance  and  become  mission-capable. 

Event  9 

This  event  occurs  only  after  the  last  flight  of  the  day.  It  ac¬ 
counts  for  the  parts  repair  process  by  subjecting  each  part  in  repair  to  the 
probability  that  the  repair  was  completed  during  the  preceding  24  hours. 
Remaining  parts  shortages  are  consolidated  on  as  few  aircraft  as  possible.  If 
the  consolidation  results  in  fewer  NMCS  aircraft  than  before,  the  aircraft 
leaving  NMCS  status  enter  maintenance  at  this  time. 

Event  10 

This  event  also  occurs  only  after  the  last  flight  of  the  day. 
Combat  losses  may  be  replaced  by  available  reserve  aircraft,  if  the  user  so 
specifies.  Any  remaining  reserve  aircraft  after  losses  have  been  replaced  are 
committed  according  to  user  specification  in  the  scenario  input  parameter.  If 
reserves  are  to  be  used  only  as  attrition  fillers,  then  any  remaining  reserve 
aircraft  are  left  i a  the  reserve  pool;  thus,  the  UE  for  the  scenario  will 
never  increase.  If  t..e  user  has  selected  the  reserve  augmentation  mode,  then 
all  reserve  aircraft  will  be  committed  when  they  become  available;  hence,  the 
UE  for  the  scenario  may  actually  increase. 

THE  REPAIR  PROCESS 

The  entry  of  an  aircraft  into  maintenance  results  from  a  ground  abort  or 
a  "break"  during  a  sortie.  In  either  case,  following  the  ground  abort  or 
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sortie,  the  aircraft  is  subjected  to  a  sequence  of  random  draws  that  deter¬ 
mines  the  subset  of  work  centers  that  will  be  involved  in  the  maintenance  on 
that  aircraft.  A  work  center  is  a  set  of  maintenance  personnel  with  a  par¬ 
ticular  skill.  Examples  of  work  centers  are  the  structural  repair  shop,  the 
hydraulic  shop,  and  the  automatic  flight  control  system  shop. 

In  the  construction  of  the  maintenance  data  base  that  supports  the  SGM, 
the  aircraft  repair  times  for  all  work  centers  involved  in  the  repair  of  the 
aircraft  are  measured  from  the  time  of  the  ground  abort  or  landing  of  the 
aircraft.  For  each  work  center  involved  in  the  repair,  a  random  draw  is  made 
from  an  exponential  distribution  of  repair  time  for  that  work  center.  The 
mean  of  i  at  distribution  is  the  reciprocal  of  the  service  rate  contained  in 
the  maintenance  data  base  for  the  work  center  in  question.  All  work  centers 
involved  in  the  repair  are  assumed  to  work  on  the  aircraft  simultaneously; 
thus,  the  recovery  time  of  the  aircraft  is  simply  the  longest  of  the  repair 
times  for  all  the  work  centers  involved  in  the  recovery  of  that  particular 
aircraft. 

In  the  SGM,  once  the  aircraft  leaves  maintenance  and  becomes  mission- 
capable  again,  it  loses  its  identity  and  is  counted  simply  as  another  aircraft 
in  the  mission-capable  pool. 
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2.  INPUT  DATA 


INTRODUCTION 

This  section  provides  a  description  of  the  inputs  required  by  the  Sortie- 
Generation  Model.  There  are  three  sets  of  inputs:  user  inputs  describing  the 
flying  scenario  and  weapon  system  to  be  simulated,  a  maintenance  manpower  file 
describing  the  characteristics  of  the  maintenance  work  centers,  and  a  recover¬ 
able  spares  file  describing  the  characteristics  of  the  LRUs  being  modeled. 
The  manpower  and  spares  files  are  outputs  of  complex  software  subsystems 
described  in  Volumes  V  and  VI. 

SCENARIO  DESCRIPTION 

At  the  beginning  of  each  SGM  run,  the  user  must  specify  various  inputs 
which  describe  the  flying  scenario  to  be  simulated.  The  mechanics  of  this 
input  process  are  described  in  Chapter  4.  A  list  of  the  inputs  and  sample 
default  values  are  shown  in  Figure  2-1.  Note  that  these  inputs  are  divided 
into  two  groups:  those  which  are  fixed  throughout  the  flying  scenario,  and 
those  that  are  allowed  to  vary  by  day  throughout  the  scenario.  The  following 
is  a  list  of  definitions  for  each  of  these  inputs: 

Number  of  Simulations 

The  number  of  independent  replications  of  the  simulation  experiment 
with  the  specified  user  inputs.  The  sortie  profile  is  computed  as  a  simple 
average  of  these  replications. 

Random  Number  Seed 

A  number  used  to  initialize  the  pseudo-random  number  generator  used 
in  the  simulation.  Changing  this  value  produces  a  different  random  number 
stream  and  allows  the  user  to  conduct  independent  experiments.  The  user  can 
reproduce  a  previous  experiment  by  using  the  same  seed  and  the  same  input 
files. 
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THE  CURRENT  VALUES  OF  THE  SCENARIO  INPUTS  ARE  : 


INPUT 

SCENARIO  ITEM 

CURRENT 

CODE 

VALUE 

1 

#  SIMULATIONS 

40 

2 

RANDOM  NUMBER  SEED 

12.3 

3 

UE 

72 

4 

AIRCRAFT  BREAK  RATE 

.20 

5 

INITIAL  NMCM  RATE 

.3 

6 

#  DAYS 

30 

7 

FIRST  TAKEOFF  TIME 

0600 

3 

LAST  TAKEOFF  TIME 

1324 

9 

SORTIE  LENGTH  (HRS) 

1.7 

10 

MINIMAL  RECOVERY  TIME 

(HRS) 

1.4 

11 

INFINITE  MANPOWER  (YES/NO) 

NO 

12 

INFINITE  SPARE  PARTS 

(YES/NO) 

NO 

13 

AUGMENT  RESERVE  AC  (YES/NO) 

NO 

14 

MAX  SORTIES/DAY  FOR  PLOT (OR  0) 

0 

THE  FOLLOWING  ITEMS  MAY  VARY 

BY  DAY ( D  > 

OR  CYCLE/ DAY 

13 

ATTRITION  RATE 

(D) 

*  .0 

16 

GROUND  ABORT  RATE 

( D ) 

»  .04 

17 

#  MASS  LAUNCHES  PER  DAY  <D) 

*  5 

18 

RESERVE  AIRCRAFT 

( D  > 

a  0 

19 

MAXIMUM  LAUNCH-SIZE 

(C/D) 

a  72 

FIGURE  2-1 

SCENARIO  INPUTS  WITH  SAMPLE  DEFAULT  VALUES 

UE  (Unit  Equipment) 

The  initial  number  of  possessed  aircraft  at  the  base  of  interest. 
This  number  does  not  include  the  available  reserve  aircraft. 

Aircraft  Break  Rate 

Probability  that  an  aircraft  returning  from  a  sortie  requires  un¬ 
scheduled  maintenance  in  one  or  more  work  centers  prior  to  further  flight. 

Initial  NMCM  Rate 

Proportion  of  the  possessed  aircraft  that  are  not-mission-capable- 
maintenance  at  the  start  of  the  flying  scenario.  For  example,  if  t&e  user 
specifies  an  initial  NMCM  rate  of  0.3  with  a  UE  of  72,  the  SGM  will  begin  each 
simulation  experiment  with  22  aircraft  undergoing  maintenance.  The  remaining 
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50  aircraft  will  initially  be  either  mission-capable  or  waiting  for  a  recover¬ 
able  spare  part. 

Number  of  Days 

The  number  of  flying  days  in  the  scenario  of  interest. 

First  and  Last  Takeoff  Times 

Scheduled  takeoff  times  for  the  first  and  last  mass  launches  of  the 
day.  These  times  are  specified  using  24-hour  clock  time,  e.g.,  1:00  o'clock 
p.m.  would  be  input  as  1300  hours.  The  remaining  scheduled  mass  launches  are 
equally  spaced  between  the  first  and  last  takeoff  times. 

Sortie  Length 

Fixed  length  in  hours  of  each  sortie. 

Minimal  Recovery  Time 

The  minimal  required  time  between  the  landing  of  the  aircraft  and 
takeoff  for  the  next  sortie,  given  that  no  corrective  maintenance  is  required. 
It  includes  only  the  time  required  to  taxi,  park,  chock,  shut  down,  refuel, 
rearm,  inspect,  restart,  and  launch. 

Infinite  Manpower 

This  input  allows  the  user  to  generate  a  sortie  profile  with  or 
without  maintenance  manpower  constraints.  (See  MAINTENANCE  MANPOWER,  page 
2-5.)  If  the  user  inputs  YES  as  the  current  value  of  scenario  item  16, 
INFINITE  MANPOWER  (YES/NO) ,  the  SGM  assumes  that  there  are  always  sufficient 
maintenance  teams  in  every  work  center  so  that  no  aircraft  ever  has  to  wait  in 
a  queue  to  begin  unscheduled  maintenance.  Thus,  the  SGM  ignores  the  number  of 
maintenance  teams  in  the  maintenance  manpower  input  file  and  uses  an  infinite 
server  assumption.  The  work  center  break  rates  and  service  rates  still  apply. 
If  the  user  inputs  NO,  the  SGM  uses  the  number  of  maintenance  teams  in  the 
input  file  as  a  constraint  on  maintenance  capacity. 
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Infinite  Spare  Parts 

This  input  allows  the  user  to  generate  a  sortie  profile  with  or 
without  spares  constraints.  (See  RECOVERABLE  SPARE  PARTS,  page  2-9.)  If  the 
user  inputs  YES  as  the  current  value  of  scenario  item  17,  INFINITE  SPARE  PARTS 
(YES/NO),  the  SGM  assumes  that  there  are  always  sufficient  spare  parts,  i.e., 
that  throughout  this  flying  scenario,  no  aircraft  is  NMCS.  Thus,  no  parts 
demands  or  repairs  are  simulated  in  this  run.  If  the  user  inputs  NO,  the  SGM 
uses  the  asset  data  from  the  spares  input  file  as  a  constraint  on  inventory. 

Augment  Reserve  Aircraft 

This  input  allows  the  user  to  specify  how  the  reserve  aircraft  are 
to  be  committed.  If  the  user  inputs  YES,  then  all  reserve  aircraft  are  com¬ 
mitted  on  the  day  they  become  available;  hence,  the  UE-size  of  the  force  may 
actually  increase  during  the  scenario.  If  NO  is  input,  then  the  reserves  are 
used  only  as  attrition  fillers  to  replace  combat  losses;  thus,  not  all  re¬ 
serves  may  be  committed  on  the  day  they  become  available. 

Maximum  Sorties/Day  For  Plot 

This  parameter  is  used  to  set  the  maximum  vertical  scale  on  the 
sorties-per-day  plot  of  the  SGM  results.  For  example,  if  the  user  wanted 
plots  of  a  series  of  SGM  runs,  he  would  use  this  parameter  to  ensure  all  the 
plots  are  on  the  same  scale.  If  zero  is  input,  the  scale  is  determined  from 
the  maximum  sorties  per  day  that  occur  in  the  SGM  results. 

Attrition  Rate 

Probability  that  an  aircraft  does  not  return  from  a  sortie  due  to 
combat  attrition.  This  rate  may  be  different  for  earh  day  of  the  scenario. 
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Ground-Abort  Rate 


Probability  that  an  aircraft  undergoes  some  failure  right  before 
takeoff  that  requires  unscheduled  maintenance  that  renders  it  not-mission- 
capable.  This  rate  may  be  different  for  each  day  of  the  scenario. 

Number  of  Mass  Launches  Per  Day 

The  number  of  equally  spaced  flying  periods  for  each  day  in  the 
scenario.  This  number  is  allowed  to  vary  for  each  day  of  the  scenario. 

Reserve  Aircraft 

The  number  of  aircraft  in  reserve  that  are  available  to  augment  the 
current  UE  of  the  scenario.  As  described  previously,  the  user  may  specify 
whether  these  reserves  are  committed  on  the  day  they  become  available  or  are 
to  be  used  only  as  attrition  fillers  to  replace  combat  losses.  The  number  of 
reserve  aircraft  arriving  on  the  scene  may  be  specified  for  each  day  of  the 
scenario . 

Maximum  Launch  Size 

The  maximum  number  of  aircraft  to  be  scheduled  on  any  particular 
wave.  The  default  value  normally  is  equal  to  the  UE  value.  This  default  is 
used  because  of  the  maximal  surge  scenario  of  the  SGM;  it  forces  all  mission- 
capable  aircraft  to  fly  each  period. 

MAINTENANCE  MANPOWER 

The  SGM  represents  aircraft  maintenance  as  a  queueing  process.  When  an 
aircraft  breaks  (malfunctions  and  cannot  take  off  for  the  next  scheduled 
sortie),  on-aircraft  unscheduled  maintenance  is  performed  by  independent, 
parallel  work  centers.  Each  work  center  consists  of  individuals  with  a  spe¬ 
cific  skill  type.  Skill  type  is  defined  by  an  individual's  Air  Force 
Specialty  Code  (AFSC) .  A  work  center,  besides  possessing  individuals  with  a 
specific  skill  type,  is  represented  by  three  other  characteristics.  These 


2-5 


characteristics  (break  rate,  number  of  servers,  and  service  rate)  are  used  to 
represent  the  work  center  in  the  queueing  model. 

The  Sortie-Generation  Model  system  consists  of  the  SGM,  a  spares  sub¬ 
system,  and  a  maintenance  subsystem.  The  purpose  of  the  maintenance  subsystem 
is  to  translate  standard  Air  Force  maintenance  performance  and  authorization 
data  into  the  queueing  model  inputs  just  described.  Generally,  the  modified 
Common  Data  Extraction  Program  (CDEP)  portion  of  the  maintenance  subsystem 
selects  records  (tasks)  based  on  the  type  of  maintenance  and  the  performing 
work  center  (AFSC).  The  Maintenance  Manpower  and  Performance  Analyzer  (MMPA) 
then  reads  the  task  data,  AFSC  authorizations  from  the  unit  manning  document, 
and  the  number  of  sorties  flown  during  the  same  time  period  to  which  the 
maintenance  data  apply.  After  reading  the  data,  the  MMPA  combines  tasks  into 
jobs  and  manipulates  the  job  data  along  with  the  number  of  sorties  and  man¬ 
power  authorizations  to  produce  the  work  center  (AFSC)  data  required  by  the 
SGM. 

Use  of  CDEP  requires  specification  of  the  workcenters  (AFSCs)  of  inter¬ 
est.  The  decision  of  what  workcenters  to  include  (model  explicitly  in  the 
SGM)  depends  on  the  work  center's  function  and  how  that  function  coincides 
with  the  SGM's  use  as  a  budgeting  tool.  Aircraft  maintenance  is  performed  by 
four  work  center  types: 

1.  Those  work  centers  assigned  to  the  Aircraft  Generation  Squadron 
(AGS)  which  perform  unscheduled  on-aircraft  maintenance. 

2.  Those  work  centers  in  the  Component  Repair  Squadron  (CRS)  or  Equip¬ 
ment  Maintenance  Squadron  (EMS)  which  perform  unscheduled  on- 
aircraft  maintenance  as  a  primary  function  or  in  support  of  the  AGS 
when  needed. 

3.  CRS  and  EMS  work  centers  whose  primary  function  is  performance  of 
scheduled  or  off-equipment  maintenance. 

4.  Those  CRS  and  EMS  work  centers  whose  on-aircraft  unscheduled  mainte¬ 
nance  can  be  deferred. 


Examples  of  the  four  work  center  types  are: 

1.  325X0  -  Automatic  Flight  Control  System 

325X1  -  Instrument  Systems 

2.  423X0  -  Electrical  System 

423X1  -  Environmental  System 

423E3  -  Fuel  System 

423E2  -  Egress  System 

3.  426R2  -  Jet  Engine  Intermediate  Maintenance  Shop 
431R1  -  Inspection  Section 

4.  427R1  -  Corrosion  Control 

Only  work  center  types  1  and  2  are  explicitly  modelled  in  the  SGM.  Work 
center  types  3  and  4  are  not  modelled  because  we  assume  that  most  scheduled 
and  deferable  maintenance  will  not  be  performed  during  a  maximal  surge  effort. 

Break  Rates 

Two  types  of  break  rates  are  used  in  the  SGM.  They  are  an  aircraft 
break  rate  and  individual  work  center  break  rates.  The  aircraft  break  rate  is 
specified  by  the  user  in  specifying  the  scenario  and  represents  an  estimate  of 
the  proportion  of  sorties  flown  in  surge  exercises  that  result  in  the  need  for 
essential  corrective  maintenance  prior  to  further  flight.  A  surge  estimate  is 
used  because  it  reflects  the  purpose  of  the  SGM;  i.e.,  estimation  of  maximal 
wartime  sortie  capability.  In  other  words,  the  surge  break  rate  incorporates 
the  tendency  for  pilots  to  hip-pocket  (not  report)  non-grounding  failures 
until  after  the  last  sortie  of  the  day.  The  result  is  a  lower  break  rate  than 
reported  during  normal  peacetime  (training)  operations. 

A  sample  F-4E  input  file  for  Seymour  Johnson  AFB  is  shown  in  Figure 
2-2.  The  reader  should  refer  to  this  figure  when  reading  the  following  work 
center  input  definitions. 

Given  a  broken  aircraft,  the  work  center  break  rate  reflects  the 
probability  that  corrective  maintenance  is  required  in  a  specific  work  center 
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******************************************** 
**********  AIRCRAFT  MAINTENANCE  ********** 
******************************************** 


BREAK 

TOTAL 

SERVICE  RATE 

c  # 

AFSC 

RATE 

SERVERS 

(ACFT/HOUR) 

1 

321X2 

0.2873 

27.77 

0. 1417 

2 

325X0 

0. 1515 

15.06 

0. 1384 

3 

323R3 

0. 1062 

31.07 

0. 1273 

4 

323X0 

0.2010 

13.  36 

0. 1769 

5 

328X4 

0. 1506 

9.57 

0.2507 

6 

404X1 

0.0225 

12.00 

0. 1510 

7 

423E2 

0. 1699 

9.95 

0.0682 

3 

423E3 

0.0608 

8.31 

0. 1043 

9 

423X0 

0.1138 

12.28 

0. 1327 

10 

423X1 

0.0793 

6.  57 

0. 1571 

11 

423X4 

0.0836 

11.91 

0. 1365 

12 

426X2 

0.0508 

10.81 

0. 1585 

13 

427R0 

0.0379 

4.56 

0.3955 

14 

427X5 

0. 1633 

14.89 

0.2584 

15 

431E1 

0. 0335 

10.73 

0.0857 

16 

431X1 

0. 0527 

131.49 

0.5356 

17 

462X0 

0. 1641 

7.27 

0.5434 

FIGURE 

2-2 

SAMPLE  MANPOWER  INPUT  FILE 
(SEYMOUR  JOHNSON  AFB,  1980:  F-NE) 

prior  to  further  flight.  Maintenance  can  be  required  in  more  than  one  work 


center  at  the  same  time. 

Estimation  of  the  work  center  break  rate  is  based  on  the  probability 


that  one  or  more  jobs  (on  the  same  aircraft)  require  corrective  maintenance  in 
the  same  work  center  following  a  sortie.  A  job  is  a  collection  of  related 
corrective  tasks.  For  example,  a  job  may  consist  of  the  following  types  of 
tasks:  troubleshoot  the  break,  remove  and  replace  the  failed  item,  and  verify 


operation  of  the  new  item. 
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Work  Center  Service  Rates 

The  work  center  service  rate  is  the  rate,  in  aircraft  per  hour, 
unscheduled  on-aircraft  maintenance  is  performed  by  a  maintenance  team  (de¬ 
fined  later  in  detail) .  The  service  rate  is  the  inverse  of  the  expected 
service  time  which  represents  the  expected  time  required  to  complete  all  jobs 
on  an  aircraft  which  has  broken  into  the  work  center.  The  job  time  begins  as 
soon  as  the  broken  aircraft  lands  or  ground  aborts,  and  includes  the  time  to 
perform  the  individual  tasks  as  well  as  the  idle  time  between  tasks.  The  job 
time  ends  when  the  last  task  in  the  job  is  completed. 

Number  of  Servers 

The  number  of  servers  represents  the  number  of  maintenance  teams  in 
a  work  center  that  are  available  to  perform  unscheduled  on-aircraft  mainte¬ 
nance.  A  server,  or  maintenance  team,  is  composed  of  one  or  more  individuals. 
Calculation  of  the  number  of  servers  is  given  by  equation  (1). 


Number  of  Servers  = 


_ Funded  Authorizations 

Expected  Number  of  Men  Working 
On  The  Aircraft 


ID 


Funded  authorizations  for  a  work  center  (AFSC)  are  obtained  from  the 
Unit  Manning  Document.  Only  type  1  and  2  work  centers,  as  described  earlier, 
are  included  in  counting  the  authorizations  for  an  AFSC.  The  CRS  and  AGS 
authorizations  are  added  together  when  the  same  AFSCs  are  found  in  both 
squadrons.  This  total  gives  the  number  of  individuals  available  to  perform 
unscheduled  on-aircraft  maintenance  in  a  surge  environment.  The  number  of  men 
working  on  the  aircraft  is  an  expected  value  that  accounts  for  the  possibility 
of  more  than  one  job  being  worked  on  (i.e.,  more  than  one  crew  at  work)  in  a 
work  center  at  once. 
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The  number  of  servers  calculated  by  equation  (1)  and  shown  in  Figure 
2-2  is  divided  by  two  in  the  SGM  to  give  the  number  of  servers  used  in  the 
queueing  model.  Division  by  two  reflects  the  availability  of  people  twelve 
hours  per  day.  Thus,  we  assume  two  twelve-hour  shifts,  each  with  the  same 
number  of  people. 

RECOVERABLE  SPARE  PARTS 

In  order  to  simulate  the  failure  and  repair  of  recoverable  LRUs  and  their 
impact  on  sortie-generation  capability,  the  SGM  needs,  for  each  LRU  that  is 
installed  on  the  aircraft  of  interest: 

a.  The  configuration  of  the  component’s  installation  on  the  aircraft 

b.  The  removal  rate 

c.  The  asset  position  at  the  start  of  the  scenario 

d.  The  average  base  resupply  time  and  the  average  depot  resupply  time 

e.  The  BNRTS  (Base  Not  Repairable  This  Station)  percentage 

The  installation  configuration,  the  BNRTS  percentage,  and  the  removal 
rate  are  obtained  directly  from  the  Air  Force  D041  data  base.  This  data  base 
contains  data  for  all  aircraft  recoverable  parts  in  the  Air  Force  inventory. 
Maintained  and  updated  quarterly  by  the  Air  Force  Logistics  Command,  this  data 
base  is  typically  used  for  making  budget  projections  and  requirements 
computations . 

The  starting  asset  position,  the  average  base  repair  time,  and  the 
average  depot  resupply  time  are  computed  by  the  LMI  Aircraft  Availability 
Model.  The  budget  levels  for  procurement  and  repair  of  components  must  be 
established  by  the  user  through  the  Interactive  Budget  Allocation  Program 
prior  to  running  the  SGM.  A  partial  listing  of  a  recoverable-spares  input 
file  is  shown  in  Figure  2-3. 
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INITIAL  RESUPPLY  TINES 


NSN 

REMOVAL 

RATE 

QPA 

FAP 

INITIAL 

STOCK 

NO.  IN 
RESUPPLY 

BASE 

NRTS 

(DAYS) 
BASE  DEPOT 

14300104546998F 

.01786 

1 

1.0 

5 

8.156 

0.08 

6.0 

26.2 

1 430010387 0388F 

.01700 

1 

1.0 

6 

7.211 

0.05 

6.0 

27.4 

2620000884523 

.01708 

2 

0.9 

90 

28.649 

1.00 

0. 

18.0 

5865001994210EN 

.00109 

4 

1.0 

1 

1.405 

0.90 

6.0 

31.0 

14300103992448F 

.01285 

1 

0.3 

1 

1.744 

0.08 

6.0 

27.7 

6610004629837BF 

00404 

1 

1.0 

6 

3.091 

0.85 

5.0 

14.7 

14300106103508F 

.00607 

1 

1.0 

3 

2.226 

0. 

6.0 

0. 

1630004463778 

.01676 

2 

1.0 

33 

9.226 

0.06 

5.0 

14.0 

1270010568980 

.00483 

1 

1.0 

4 

0.937 

0.45 

6.0 

11.5 

5826010395000 

.00600 

1 

0.3 

1 

1.373 

0.30 

9.0 

22.1 

5326010401785 

.00692 

1 

0.3 

2 

1.650 

0.40 

8.0 

21.0 

143001 0387055BF 

.00484 

1 

1.0 

3 

1.742 

0. 

6.0 

0. 

5826010183511 

.00168 

2 

1.0 

9 

3.878 

0.60 

4.0 

34.6 

1430002356325BF 

.01023 

1 

0.9 

6 

2.282 

0.09 

3.0 

17.1 

2S40006717414PL 

.00084 

2 

1.0 

2 

2.431 

0.76 

6.0 

28.6 

5865000233292EU 

.00200 

l 

1.0 

0 

0.005 

1.00 

0. 

6.8 

586500371334401 

.00133 

4 

1.0 

3 

0.723 

0.97 

3.0 

13.3 

1270000641997 

.00441 

1 

1.0 

6 

1.345 

0.66 

4.0 

10.9 

6115008681999EU 

.00184 

5 

0.2 

I 

0.092 

0.43 

2.0 

78.8 

586500099934801 

.00103 

5 

1.0 

2 

0.548 

0.98 

2.0 

12.3 

FIGURE  2-3 

SAMPLE  RECOVERABLE  SPARES  INPUT 


Installation  Conf iguration 

Not  all  components  are  installed  one  per  aircraft.  On  some  aircraft 
in  which  a  component  is  used,  more  than  one  may  be  installed.  Furthermore, 
not  every  aircraft  of  a  given  type  will  be  in  a  standard  configuration.  Two 
elements  of  input  data  are  used  to  define  a  component's  installation  config¬ 
uration,  the  quantity  per  aircraft  (QPA)  and  the  future  application  percentage 
(TAP).  The  QPA  defines  the  component's  standard  configuration  with  respect  to 
a  specific  model/design/series  aircraft.  The  FAP  is  used  for  a  variety  of 
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purposes  in  Che  D041  system  but  is  interpreted  by  the  SGM  to  denote  the  pro¬ 
portion  of  aircraft  with  the  component  installed  in  standard  configuration. 
Thus,  component  flying  hours  may  be  computed  by  multiplying  aircraft  flying 
hours  times  QPA  times  FAP. 

Removal  Rate 

The  component  removal  rate  is  defined  as  the  expected  number  of 
component  removals  per  component  flying  hour.  It  is  computed  in  a  straight¬ 
forward  manner  using  the  daily  demand  rate  (DDR)  from  the  D041  data  base  and 
the  peacetime  flying  hours  for  all  the  aircraft  types  which  use  the  component. 
An  important  assumption  of  the  SGM  system  is  that  the  expected  number  of 
component  removals  is  directly  proportional  to  the  number  of  flying  hours 
flown. 

Starting  Asset  Position 

The  starting  asset  position  for  a  component  is  defined  by  the  ex¬ 
pected  number  in  resupply  at  the  start  of  the  scenario  and  the  authorized 
number  of  spares.  By  resupply  we  mean  base  repair,  depot  repair,  or  in 
transit.  If  the  authorized  number  of  spares  is  greater  than  the  starting 
number  in  resupply,  then  there  will  be  spares  on  hand  at  the  start  of  the 
scenario.  If  the  authorized  number  of  spares  is  less  than  the  starting  number 
in  resupply,  then  there  will  be  unsatisfied  (backordered)  demands  at  the  start 
of  the  scenario.  The  starting  number  in  resupply  is  determined  independently 
for  each  replication  run  by  the  SGM.  The  model  does  this  by  drawing  at  random 
a  starting  number  from  a  Poisson  probability  distribution  whose  mean  is  the 
expected  number  in  resupply. 

Resupply  Times 

The  average  base  resupply  time  is  the  average  number  of  days  that  a 
component  is  in  the  base  repair  shop  before  it  is  returned  to  servicable 
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status.  This  equals  the  shop  repair  time  plus  the  average  time  awaiting  SRUs 
(if  any).  The  average  depot  resupply  time  is  the  average  number  of  days 
between  when  an  order  for  a  component  is  made  to  the  depot  and  when  the  com¬ 
ponent  is  received  at  the  base.  This  equals  the  order  and  ship  time  plus  the 
average  delay  at  the  depot  awaiting  an  available  spare  to  ship. 

BNRTS  Percentage 

The  percentage  of  demands  which  are  not  base  repairable.  A  demand 
which  is  not  base  repairable  may  be  condemned  or  repaired  at  the  depot  but  in 
either  case  an  order  will  be  placed  for  depot  resupply. 

The  spares  subsystem  excludes  many  components  from  the  SGM  inputs  on 
the  basis  that  they  will  not  have  any  impact  on  the  estimated  sortie- 
generation  capability.  Components  are  excluded  if  their  removal  rates  are 
Less  than  .0005  per  hour.  All  components  whose  demands  are  not  flying-hour 
dependent  are  also  excluded. 
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3 .  OUTPUTS 


INTRODUCTION 


The  Sortie-Generation  Model  outputs  consist  of  a  summary  of  the  user- 
specified  scenario,  a  sortie  profile  showing  the  average  number  of  sorties 
flown  for  each  period  of  the  scenario,  and  a  graph  of  the  daily  sortie  produc¬ 
tion. 

SCENARIO  SUMMARY 

A  summary  of  the  user  inputs  describing  the  scenario  being  simulated  is 
printed  at  the  beginning  of  each  SGM  run.  A  sample  of  this  summary  is  shown 

in  Figures  3-1  and  3-2.  The  first  section  describes  dimensions  of  the 

♦♦♦■a*************************************************************** 
****************************  SGM  RUN  **************************** 
******************************************************************* 


SIMULATION  -  REPLICATIONS  -  40  RANDOM  NUMBER  SEED  ■  12.3 

AIRCRAFT  -  UE  =  72  RESERVES  *  24  MAXIMUM  LAUNCH-SIZE  =  72 


FLYING  SCHEDULE  - 


WAVES 

DAYS  PER  DAY 

TAKEOFF 

FIRST 

TIMES 

LAST 

MINIMAL 

TURNAROUND 

SORTIE 

LENGTH 

WAIT 

TIME 

OVERNIGHT 

RECOVERY 

30  5 

0600 

1824 

1.40 

1.70 

0.00 

8.50 

RATES  - 

INITIAL 
NMCM  RATE 

ATTRITION 

AIRCRAFT 

BREAK  RATE  GROUND-ABORT 

0.  150 

0.01 

0. 

2000  0. 

0400 

LRU  TYPES  -  262 


FIGURE  3-1 
SCENARIO  SUMMARY 
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'I 


***4H»*#*«*************«*«******»***4MHHHHH»** 

****■*■*■#**#  AIRCRAFT  MAINTENANCE  #**♦*•»**** 


wc  * 

AF3C 

BREAK 

RATE 

TOTAL 

SERVERS 

SERVICE  RATE 
<  ACFT /HOUR ) 

1 

321X2 

0.2878 

27.77 

0. 1417 

2 

325X0 

0. 1515 

15.06 

0. 1384 

3 

32SR3 

0. 1062 

31.07 

0. 1273 

4 

323X0 

0.2010 

18.36 

0. 1769 

5 

328X4 

0. 1506 

9.57 

0.2507 

6 

404X1 

0.0225 

12.00 

0. 1510 

7 

423E2 

0. 1699 

9.95 

0.0682 

3 

423E3 

0.0608 

8.31 

0. 1043 

9 

423X0 

0.1188 

12.28 

0. 1327 

10 

423X1 

0.0793 

6.57 

0. 1571 

11 

423X4 

0.0836 

11.91 

0. 1365 

12 

426X2 

0.0508 

10.81 

0. 1585 

13 

427R0 

0.0379 

4.56 

0.3955 

14 

427X5 

0. 1633 

14.39 

0.2584 

15 

431E1 

0.0335 

10.73 

0.0857 

16 

431X1 

0.0527 

131.49 

0.5356 

17 

462X0 

0. 1641 

7.27 

0.5434 

FIGURE  3-2 
WORK  CENTER  SUMMARY 


simulation,  numbers  of  aircraft,  flying  schedule,  and  various  rates  and  prob¬ 
abilities.  The  second  section  lists  the  work-center  characteristics  provided 
in  the  maintenance-manpower  input  file.  Definitions  of  these  scenario  and 
work-center  inputs  are  provided  in  Chapter  2. 

SORTIE  PROFILE 

The  sortie  profile  is  a  summary  of  the  average  numbers  of  aircraft  in  the 
various  states  modeled  by  the  SGM:  Mission-capable,  Maintenance,  NMCS,  Combat 
Loss,  or  Reserve.  The  numbers  of  aircraft  in  each  SGM  state  are  collected  at 
the  beginning  of  each  sortie  period  for  each  flying  day,  and  an  average  is 
computed  from  the  total  for  all  simulation  replications.  A  sample  of  a  sortie 
profile  is  shown  in  Figure  3-3.  A  short  description  of  the  information  pro¬ 
vided  by  each  output  column  follows: 
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FIGURE  3-3 
Sm  SORTIE  PROFILE 
(FIRST  FIVE  DAYS) 
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Day  and  Period 

The  DAY  and  PER  columns  indicate  the  flying  day  and  period, 
respectively,  for  which  the  averages  were  collected. 

Sorties  per  Period 

The  SORTIES/PERIOD  column  indicates  the  average  number  of  sorties 
flown  for  this  day  and  period. 

Sorties  per  Day 

The  SORTIES/DAY  column  reflects  the  total  number  of  sorties  flown  on 
each  flying  day.  A  cumulative  sortie  total  for  all  days  thus  far  is  also 
printed  immediately  below  the  total  for  each  day. 

Sorties  per  Aircraft 

The  SORTIES/AC  column  provides  an  indication  of  the  average  number 
of  sorties  being  flown  by  each  on-the-scene  a: rcraft  for  this  particular 
flying  day.  An  on-the-scene  aircraft  is  defined  as  an  aircraft  which  is 
either  mission-capable,  in  maintenance,  or  NMCS;  i.e.,  it  is  not  in  the  COMBAT 
LOSS  or  RESERVE  state.  This  statistic  is  computed  by  dividing  the  average 
number  of  sorties  flown  each  day  by  the  average  number  of  on-the-scene  air¬ 
craft  for  that  day.  This  on-the-scene  average  is  computed  by  totaling  the 
number  of  on-the-scene  aircraft  at  the  start  of  each  flying  period  for  the 
day . 

Maintenance 

The  NMCM  column  reflects  the  average  number  of  aircraft  that  were 
unable  to  fly  each  period  because  they  were  in  maintenance  at  the  start  of  the 
pre-takeoff  period.  This  pre-takeoff  period  begins  at  a  time  T^  before 
scheduled  takeoff  (see  Figure  1-3).  This  average  includes  any  aircraft  which 
may  have  been  repaired  during  the  pre-takeoff  period  and  it  does  not  include 
any  aicraft  which  may  have  ground-aborted  before  takeoff. 
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NMCS 


This  is  the  average  number  of  aircraft  that  were  unable  to  fly  each 
period  because  they  were  in  the  NMCS  state  at  the  start  of  the  pre-takeoff 
period . 

Cumulative  Combat  Losses 

The  CUM.  LOSSES  column  indicates  the  cumulative  number  of  combat 
losses  at  the  beginning  of  this  flying  period. 

Remaining  Reserves 

The  RES.  LEFT  column  is  the  average  number  of  reserve  aircraft 
remaining  at  the  start  of  each  flying  day.  Since  reserve  aircraft  are  com¬ 
mitted  only  once  at  the  end  of  each  flying  day  to  replace  combat  losses,  this 
average  is  only  printed  for  this  first  flying  period  of  the  day  to  indicate 
the  remaining  reserves. 

SORTIE  PLOTS 

In  addition  to  the  sortie  profile  for  each  SGM  run,  graphs  of  the  average 
sorties  per  aircraft  per  day  and  average  sorties  flown  are  printed.  Samples 
of  these  plots  are  shown  in  Figures  3-4  and  3-5. 
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SGM  SORTIE  PLOT:  SORTIE  RATES 
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DAY  OF  SCENARIO 

FIGURE  3-5 

SGM  SORTIE  PLOT:  SORTIES  PER  DAY 
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4.  RUN  INSTRUCTIONS 


ENVIRONMENT 

The  Sortie-Generation  Model  (SGM)  has  been  developed  on  System  C,  an 
unclassified  computer  system  located  at  the  Pentagon  and  supported  by  the  Air 
Force  Data  Services  Center  (AFDSC).  This  system  operates  on  a  Honeywell  G-635 
computer  under  the  series  600/6000  GCOS  Time-Sharing  System.  Access  to  the 
system  is  possible  on  remote  terminals  by  a  dial-up  procedure. 

The  SGM  has  been  written  in  the  Honeywell  600/6000  FORTRAN  programming 
language,  the  only  version  of  FORTRAN  available  on  the  system.  The  run  pro¬ 
cess  has  been  designed  so  that  the  model  may  be  run  in  either  the  remote-batch 
or  time-sharing  modes.  There  are  advantages  and  disadvantages  to  both  pro¬ 
cedures.  If  System  C  is  carrying  a  light  load  (i.e.,  only  a  few  users  are 
signed  on),  then  a  time  sharing  run  is  significantly  faster;  however,  through¬ 
out  the  simulation  the  terminal  cannot  be  used  for  any  other  purpose  and  it  is 
not  possible  to  direct  the  output  elsewhere.  Once  a  job  has  been  submitted 
interactively  to  be  run  as  a  batch  job,  the  user  is  free  to  make  other  runs, 
use  the  terminal  for  some  other  purpose,  or  even  to  log-off  the  computer. 

For  a  detailed  description  of  System  C,  the  Time-Sharing  System,  and 
FORTRAN  600/6000,  the  user  is  referred  to  the  following  manuals: 

a.  AFDSC  User's  Handbook,  Volume  III  -  G-635  Computer  Systems 

b.  Honeywell  Series  600/6000  Time-Sharing  System  General  Information 

Manual 

c.  Honeywell  Series  600/6000  FORTRAN  Manual 
RUN  STEPS 

There  are  three  steps  in  obtaining  an  SGM  run  (see  Figure  4-1):  initiate 
communication  with  the  Time-Sharing  System  (log-on  System  C);  specify  scenario 
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SIERS 


FIGURE  4-1 
RUNNING  THE  SGM 


inputs  for  the  run;  and  initiate  execution  of  the  SGM  in  either  the  remote- 
batch  or  time-sharing  mode.  Figure  4-2  is  a  list  of  the  exact  commands  to  be 
typed  by  the  user  for  a  batch  or  time-sharing  run.  Appendix  A  provides  list¬ 
ings  of  the  Job  Control  Language  (JCL)  and  Run-Command  files  which  are 
initiated  by  these  listed  commands. 
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-  Log-on  System  C  - 

RUNY  LA61A/SUBMLR 

RUNC  0S29/N232D/S6M/SETPARAM 

-  Set  Scenario  Parameters 

RUNC  0S29/N232D/SGM/ RS6MTSS 


BATCH  m  COMMAND  SEQUENCE 

-  Log-on  System  C  - 

RUNY  LA61A/SUBMIT,R 

RUNC  0S29/N232D/SGM/SETPARAM 

-  Set  Scenario  Parameters 

RUN  0S29/N232D/SGM/RSGMBTCH 


FIGURE  A-2 

SGH  RUN  INSTRUCTIONS 


The  remainder  of  this  section  provides  a  description  of  each  step  of  the 
SGM  run  process. 

Log-on  System  C 

To  initiate  communication  with  the  System  C  Time-Sharing  System,  the 
user  must  connect  an  appropriate  remote  terminal  to  System  C  and  enter  a 
unique  user  identification  code,  a  verification  password,  and  the  appropriate 
project  identification  code.  These  codes  are  assigned  by  AFDSC.  A  complete 
description  of  the  log-on  procedures  for  the  various  types  of  remote  terminals 
is  provided  in  the  Honeywell  Series  600/6000  Time-Sharing  System  General 
Information  Manual  referenced  previously. 

Specifying  Scenario  Inputs 

The  next  step  for  either  a  batch  or  time-sharing  SGM  run  is  to 
specify  the  scenario.  The  scenario  inputs  are  defined  in  Chapter  2  of  this 
guide.  Entry  of  the  scenario  parameters  is  performed  via  an  interactive 
program  with  a  question-and-answer  format.  The  following  commands  initiate 
this  interactive  program: 

RUNY  LA61A/SUBMIT,R 

RUNC  0S29/N232D/ SGM/ SETP ARAM 

The  user  is  immediately  prompted  for  the  aircraft  type  which  is  used  to  select 
the  appropriate  file  containing  the  default  scenario  values  for  that  aircraft. 
(Default  files  currently  exist  for  the  following  aircraft  types:  A-10,  F-4, 
F-15,  F-16,  and  F-lll.)  Next,  the  user  is  asked  for  a  random-number  seed  to 
initialize  the  random-number  generator.  Then,  the  default  scenario  parameters 
are  loaded,  and  the  user  is  allowed  to  list  or  change  these  scenario  param¬ 
eters.  If  any  changes  are  desired;  the  user  specifies  a  parameter  to  be 
changed  and  the  corresponding  new  value.  This  process  continues  until  the 
user  has  made  all  the  desired  changes;  the  scenario  program  is  then  terminated 
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and  the  SGM  can  be  run  in  either  a  batch  or  time-sharing  mode.  Figure  4-3 
provides  an  example  of  an  interactive  scenario-input  session.  The  user 
responses  for  this  session  are  underlined  and  numbered,  and  an  explanation  of 
each  response  is  provided  in  Figure  4-4. 

Time-Sharing  Run 

To  run  the  SGM  in  this  mode,  the  user  types:  RUNC  OS29/N232D/ 
SGM/RSGMTSS.  This  program  asks  the  user  for  the  maintenance-manpower  file  and 
recoverable-spares  file  (see  Chapter  2  for  description  of  these  files).  The 
SGM  then  begins  execution  and  the  results  (see  Chapter  3  for  description  of 
outputs)  are  printed  at  the  user's  remote  terminal. 

Batch  Run 

To  run  in  batch  mode,  the  user  types:  RUN  OS29/N232D/SGM/RSGMBTCH . 
The  program  asks  for  the  maintenance-manpower  file  and  recoverable-spares  file 
as  described  above.  After  the  job  is  submitted  to  System  C  batch,  a  five¬ 
digit  identification  number  is  printed  to  enable  the  user  to  track  the  job 
while  it  is  being  processed.  The  SGM  output  is  directed  to  the  printer  sta¬ 
tion  specified  by  the  JCL  run  file  (currently  set  for  the  LMI  remote  printer). 

Log-off  Procedure 

Termination  of  an  SGM  run  session  is  a  two-step  procedure.  First, 
the  user  types  DONE  when  prompted  by  an  "="  after  the  end  of  the  SGM  run;  this 
terminates  the  interactive  submission  program.  Then,  the  user  types  BYE  to 
disconnect  the  terminal  from  System  C  Time-Sharing. 
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FIGURE  4-4 

EXPLANATION  OF  SCENARIO- INPUT  SESSION 


APPENDIX  A 


SGM  JOB  CONTROL  LANGUAGE  (JCL) 

This  appendix  provides  listings  of  the  JCL  and  Run-Command  files  used  to 
run  the  Sortie-Generation  Model  (SGM): 

0S29/N232D/SGM/RSGMBTCH 

This  is  the  JCL  file  for  running  the  SGM  in  the  Batch  mode.  Output  is 
directed  to  remote  printer  at  LMI . 

0S29/N232D/SGM/RSGMTSS 

This  is  the  Run-Command  file  for  running  the  SGM  in  the  Time-Sharing  mode 
under  the  LMI  SUBMIT  subsystem. 

0S29/N232D/SGM/SETPARAM 

This  is  the  Run-Command  file  for  setting  the  scenario  parameters  (in 
Time-Sharing  mode)  under  the  LMI  SUBMIT  subsystem. 

0S29/N232D/SGM/ZD.F4 

This  is  the  file  which  provides  the  default  values  for  initializing 
scenario  parameters  for  the  F-4. 


i 

! 


0S29/N232D/SGM/RSGMBTCH 
100##S, R ( XL )  : >Si 16, 58 

110*: NOTE:**  MIKE  **  0S29/N232D/SGM/RSGMBTCH 

120*s I DENT  :  0S20 1 1 N24 1 D  , OS29UGOODWIN 

1 30*: OPT I ON : FORTRAN , NOMAP 

140*: SELECT: 0S29/N232D/SGM/CSGM 

ISO*: EXECUTE 

160*: LIMITS: 14,27K,  ,  IK 

170*: DATA: 01 , NCKSUM, COPY 

ISO*: SELECTD: 0S29/N232D/SGM/PARAMS 

190*: ENDCOPY 

200* : PRMFL : 02 , R , S , 0S29 / N24 1 D/CDEP / SGM I NPT2 / &MANPOWERB ASE . 

210*: PRMFL: 04, R, S, LA61 A/SLAY/DATA/&AC-TYPE. / &SP ARESF I LE . 

220*: FILE: 03, A3S 

230*: FILE: 07, A IS 

240* : OPT I ON : FOR  TRAN , NOMAP 

250*: SELECT: 0S29/N232D/SGM/CPL0T 

260*: EXECUTE 

270*: LIMITS: 1 , 13K, , 2K 

280*: FILE: 07, AIR 

290*: ENDUOB 


0S29/N232D/SGM/RSGMTSS 


REMO  CLEARFILES 

GET  0S29/N232D/SGM/PARAMS"01 " . R 

GET  0S29/N24 1 D/CDEP/ SGM I NPT2/ &MANPOWERBASE .  ,,02",R 
GET  LA6 1  A/ SLAY /DATA/  & AC-TYPE .  /  &SP ARESF I LE .  "04H>R 
TEMP  07503 

RUNY  0S29/N232D/SGM/CSGM. R 
RUNY  0S29/N232D/SGM/CPL0T.R 
REMO  CSGM5 075 035 CPL0T5 015025 04 


0S29/N232D/SGM/SETPARAM 

REMO  CLEARFILES 
TEMP  01 

GET  0S29/N232D/SGM/ZD. &AC-TYPE. 
RUNY  0S29/N232D/SGM/HZDATA,R 
PERM  01  ?  0S29/N232D/SGM/PARAMS 
REMO  PARAMS*  HZDATA! 08 
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“OS’SR 
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0S29/N232D/SGM/ZD. F4 

*ZDATA, 

INFMAN='NO", 

NONORS='NO-', 

NSIM=',40-'. 

UE='72'. 

MAXFLY="72', 

RES='Oy  j 
ATTRIT='.0' , 
ANYBRK='.20% 
ANYGA='.04', 

RNMCM= ' .3' , 

NUMDAY='30' , 

NCYCLE=  5', 
FTOTYM='0600', 

LTOTYM=' 1324 
PREFLT='1.4' , 
SRTLTH='1.7 
SCALE= "O' » 

I AUGMT  =  •"  NO  "  i 

* 


APPENDIX  B 


SGM  PROGRAM  LISTING 


This  appendix  provides  a  listing  of  the  main  program,  subroutines,  func¬ 
tions,  and  block  data  subprograms  comprising  the  Sortie-Generation  Model.  The 
main  program  is  listed  first.  The  subroutines,  function,  and  block  data,  each 
beginning  on  a  new  page,  are  then  listed  in  alphabetical  order. 
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Cwhwh  0S29/N232D/SGH/NEWSGH 

CH<HWWHWWHHWHHHHHHWHHHWHHHHWW<HW<HW 

Cm  MAIN  PROGRAM 

CHIHHHHIWH<HHHW<H<HIHW<H<HIH<HHHHH<H . . 

C++  MAIN  -  MAIN  PROGRAM  FOR  LMI  SORTIE-GOCRATICN  MODEL  (SGM). 
Cm  THIS  IS  TIC  MAIN  PROGRAM  FOR  THE  LMI  SORTIE  GENERATION 
Cm  MODEL  (SGM).  A  LIST  OF  ALL  COMMON  BLOCKS  Mffi  PARAMETER  STATEMENTS 
Cm  USED  IN  THE  MODEL  IS  PROVIDED  AT  THE  BEGINNING  OF  THIS  MAIN 
Cm  PROGRAM.  T>€  PROCESSING  SEQUENCE  IS  AS  FOLLOWS  -  FIRST. 

Cm  ALL  INPUTS  ARE  LOADED  ttffl  NECESSARY  INITIALIZATION  PERFORMED. 

Cm  THEN,  T>£  ACTUAL  SIMULATION  IS  PERFORMED,  AND  FINALLY  T>€  RUN 
Cm  RESULTS  ARE  PRINTED  TO  THE  STANDARD  OUTPUT  FILE. 

Cm 

Cm  INPUT  FILES  - 

Cm  01  -  SCENARIO  INPUT  PARAMETERS 

Cm  02  -  WORK  CENTER  INPUT  DATA 

Cm  03  -  SCRATCH  FILE  USED  FOR  DAILY  SCENARIO  PARAMETERS 

Cm  04  -  SPARES  INPUT  DATA 

Cm  OUTPUT  FILES  - 


Cm  06  -  STANDARD  OUTPUT  FILE  (RUN  RESULTS) 
Cm  07  -  SORTIE  RESULTS  FOR  PLOT  PROGRAM 
Cm  PARAMETERS  - 


Cm 

MAXAC 

Cm 

MXWC 

Cm 

MAXBIT 

Cm 

MAXPRT 

Cm 

HAXVEC 

Cm 

Cm 

Cm 

Cm 

Cm 

MAXDAY 

Cm 

MAXCYC 

Cm 

MAXSTAT 

Cm 

Cm 

LFLD 

Cm 

Cm 

Cm 

Cm 

Cm 

NPERURD 

Cm 

Cm 

Cm 

MXINWC 

Cm 

Cm 

Cm 

Cm 

IFSCEN 

Cm 

IFWC 

Cm 

IFPRT 

-  MAXIMUM  ALLOWABLE  UE-STRENGTH  (•  AIRCRAFT) 

-  MAXIMUM  ALLOWABLE  NUMBER  OF  WORK  CENTERS 

-  NUMBER  OF  BITS  IN  A  COMPUTER  WORD  ON  THIS  SYSTEM 

-  MAXIMUM  ALLOWABLE  NUMBER  OF  PART-TYPES. 

-  MAXIMUM  ALLOWABLE  LENGTH  (IN  COWUTER  WORDS)  OF 
AIRCRAFT  BIT-VECTORS.  A  BIT-VECTOR  MUST  BE  AT 
LEAST  'MAXAC'  BITS  LONG,  PLUS  AN  EXTRA  WORD 

TO  STORE  THE  AIRCRAFT  COUNT  FOR  THAT  VECTOR. 

HENCE,  HAXVEC  IS  A  FUNCTION  OF  IttXAC  AND  MAXBIT. 

-  MAXIMUM  ALLOWABLE  NUMBER  OF  FLYING  DAYS. 

-  MAXIMUM  ALLOWABLE  NIGER  OF  FLYING  CYCLES  PER  DAY. 

-  CURRENT  NIGER  OF  STATISTICS  COLLECTED  PER 
FLYING  CYCLE  PER  DAY. 

-  LENGTH  OF  BIT-FIELD  USED  IN  TIC  WORK-CENTER 
REPAIR  LISTS.  THIS  BIT-FIELD  MUST  BE  LARGE  ENOUGH 
TO  HOLD  (MAXAC-1),  THE  TAIL  NIGER  OF  T>€ 

LAST  AIRCRAFT.  THUS,  <2«tFLD)  MUST  BE  GREATER 
THAN  OR  EQUAL  TO  MAXAC. 

-  NIGER  OF  BIT-FIELDS  PER  COMPUTER  WORD  FOR  TfCSE 
WORK-CENTER  LISTS. THUS  N’ERWRD  IS  A  FUNCTION 

OF  LFLD  AND  MAXBIT, 

-  LENGTH  (IN  COMPUTER  WORDS)  OF  TTC  WORK-CENTER  LISTS. 
MXINWC  IS  COMPUTED  SO  THAT  THE  MAXIMUM  ALLOWABLE 
NUMBER  OF  BIT  FIELDS  IN  A  WORK-CENTER  LIST  IS 
EQUAL  TO  MAXAC,  TO  MAXIMUM  NUMBER  OF  AIRCRAFT. 

-  FILE  NUMBER  OF  SCENARIO  INPUT  FILE 

-  FILE  NIGER  OF  WORK  CENTER  INPUT  FILE 

-  FILE  NIGER  OF  SPARES  INPUT  FILE 


C«m»Hninm<mi»i<<ii<i<»mm<mmnimim»»»m<<»*m 

C— 


PWNETER  MAXACsl06iHAXWC=25,HAXBIT=36,HAXPRT*304, 
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O  <?  CP  <~>  o  o  «.  o  o  «•  o  o  *•  o  o  <p  o  *•  «*  «•  no**  on  oo*»  or>  o  r>  **  <-><-» 


l  HMS€Ca2>  ( MAX AC- 1 )  /NAXB  IT 

PWWCTER  LFUW.NPERWRIHIAXBIT/LFLDiMXINUCsl+IHAXAC-n/fPERHRD 
PflRAJCTEF  NAXDAY*30,MAXCYC=10.MAXSTAT»5 
PARAfCTER  IFSC0W1,  IFUC=02,  IFPRT=04 
LOGICAL  INFNAN, IlfPART 

-/ACSTATE/  -  AIRCRAFT  BIT- VECTORS. 

COWON  /ACSTATE/  LENGTH.  NACVCIMAXVEC),  IRYVC(NAXVEC). 

HAINVC(HAXVEC),  NORSVC(HAXVEC),  LOSTVC ( MAXVEC) 

-/ALIASC/  -  TABLES  FOR  PART-TYPE  SAMPLING. 

COWON  /ALIASC/  FRACTIHAXPRT),  IALIASOWXPRT),  FPMTS 

-/BITS/  -  BIT  MANIPULATION  TABLES. 

COWON  /BITS/  HASK0,MASK(35)<  HL£FT0,MSKLFT(36)> 

IZCOUT,  ICOJNT<43) 

-/DEMAND/  -  MEAN  NO  VttINCE  FOR  TOTAL  P«T  DEIWOS. 

COWON  /OENAM)/  ACfEAN.  AC VAR,  NPERAC 

-/INPUT/  -  FLYING  SCEJMRIO  PARAMETERS. 

COWON  /INPUT/  INITUE,  NAC.  PATTRIT,  IRES,  RNMCM,  INFPWT, 
MAXFLY(MAXCYC),  HPVW,  ISCALE,  IAUGWT 

-/PARTS/  -  PART  CHARACTERISTICS. 

COWON  /PARTS/  l*WTS,  IQPA(MAXPRT),  NBACKO(MAXPRT), 

BRPRATE(MAXPRT),  DRPRATE(MAXPRT),  INITSJ(NAXPRT), 
RESUPP(NAXPRT),  BWTS/WXPRTK  NBASEOttXPRT), 
NDEPOT(HAXPRT) 

■/RSEED/  -  SEED  FOR  RAWOM  NllfflER  GENERATOR. 

COWON  /RSEED/  SEED 

•/STATS/  -  CUMULATIVE  STATISTICS  FOR  SIMULATION  RESULTS. 

COWON  /STATS/  EXPECT  ( HAXSTAT ,  HAXCYC ,  MAXDAY  > , 

WESRV,  I ZDAY ,  ITOTRES (MAXDAY > ,  LOSSTOT 

■/TI(C/  -  FLYING  CYCLE  TIMES  AND  SIMULATION  PARAMETERS. 

COMMON  /TIME/  PRERITE,  SORTLGTH,  MAITCYC,  TYWITE, 

NSIH,  IS1M,  NUMDAY,  IDAY,  NCYdfS,  ICYCLE 

-/HC8RK/  -  WORK  CENTER  BREAK  RATES. 

COWON  /UCBRK/  PACBRK,  PACGABT,  PBRKWC(MAXUC),  PUCPROD, 
PBRKSEQ(2,MAXUC),  HCXUCIMAXMC) 

•/HCINPUT/  -  UORK  CENTER  INPUTS. 

COWON  /UCINPUT/  NUC,  NCREWS(MAXHC),  SRATE(MAXHC) 

/UCMAIN17  -  AIRCRAFT  UORK  CENTER  LISTS. 

COWON  /UCMAINT/  LISTRP(MXINHC.MAXHC),  IffiEFRfMXWO 

•COLLECT  STARTING  CPU-TH€  AND  CORE-fCNORY  REQUIREMENT 
CALL  PTIIC(START) 
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Cfl-L  fdSIZ(KSIZE) 


C — 

C —  H.QAD  AM)  INITIALIZE  SCENARIO.  WORK  CENTER  AW)  PARTS  DATA 
CALL  INIT ( IFSCEN? IFWC. IFPRT) 

C— 

C—  »RUN  TIC  ACTUAL  SIMULATION 
CALL  SIMULA 

C — 

C —  «PRINT-OUT  TIC  RESULTS  OF  THE  SIMULATION 
C#1  PRINTO 

C — 

C—  ♦PRINT  MEMORY  AM)  CPU-TIIC  USED 
CALL  PTIHE(FINISH) 

NRITE(6,900i ) (FINISH-START)#60.  .KSIZE 

C — 

STOP 

9001  FORMAT! //.'OCPU  TIIC  USED  =*,F6.2,*  MIN1,/, 

&  -OfCMORY  USED  =M6.*  K  UORDS") 

EM) 
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CimHIIIHHIIIIIIHllimimillHItHHimimWHIlHIHHIHH 

SUBROUTINE  ALIAStN.FRACT. 1ALIAS) 

CiiiimiiiHHiniiHinmiititwHwwmiiHiiiniiiiiiiiiini 

C++  ALIAS  -  INITIALIZE  TABLES  NEEDED  FOR  ‘M.IAS*  SAWLING  METHOD. 
CM*  ALIAS  IS  A  FORTRAN  SUBROUTINE  WHICH  INITIALIZES  THE 
Cm  TABLES  USED  BY  THE  ALIAS  fCTHOD  FOR  GENERATING  RANDOM 
CM*  VARIABLES  FROM  A  DISCRETE  DISTRIBUTION.  SEE  -  "ON  THE 
Cm*  ALIAS  fETHOD  FOR  GENERATING  RNOOM  VARIABLES  FROM  A  DISCRETE 
CM*  DISTRIBUTION*  IN  THE  AMERICAN  STATISTICIAN.  NOV  1979,  VOL  33, 

CM*  NO  4,  PP  214-218,  FOR  A  DESCRIPTION  OF  THIS  fETHOD  M  Tf€ 

C*»  ALGORITHM  USED  IN  THIS  ROUTINE  TO  CREATE  THE  NECESSARY  TABLES. 
Cm*  THO  TABLES  <«E  NEEDED  FOR  THIS  METHOD  -  A  TABLE  OF 
CM*  FRACTIONAL  CUTOFF  VALUES  AND  ANOTHER  FOR  TIE  CORRESPONDING 
CM*  ALIASES.  THE  PROCEDURE  USED  TO  GEfERATE  THESE  TABLES  IS  A 
C***  SINGLE-PASS,  LINKED-UST  PROCEDURE. 

Cm* 


Cm*  IffUT  - 

CM*  N  -  NUMBER  OF  fttSS  POINTS  OF  THE  DISCRETE  DISTRIBUTION 

Cm*  WHICH  IS  BEING  SAfPLED. 


C*** 

Cm* 

Cm* 

Cm* 

Cm* 

CM* 

Cm# 


I  WITT /OUTPUT  - 

FRACT(I)  -  IRON  INPUT,  FRACT(I)  IS  THE  PROBABILITY 
DISTRIBUTION  OF  A  RANDOM  VARIABLE,  R. 
FRACT(I)=PROBABILITY(  R  =  I),  1=1, 2, ...,N  . 

UPON  OUTPUT  FROM  THIS  SUBROUTINE,  FRACT  CONTAINS 
TIC  TABLE  OF  FRACTIONAL  CUTOFF  VALUES  USED  BY  THE 
ALIAS  METHOD. 


Cm*  OUTPUT  - 

Cm*  IALIAS(I)  -  TABLE  OF  ALIASES  USED  BY  ALIAS  METHOD.  1=1, ...,N 


CimHHWHIHIHHHWHHtWHHWHHUIHtHHIHHIIHHim 


C— 

DIfENSION  FRACT(N),  IALIAS(N) 


C— 

C —  *INITIALIZE  LIST  )£ADERS  TO  NO  ENTRIES 

LHEAD  =  0 
NCAD  =  0 

C — 

C —  *DO  FORtEACH  POINT  OF  Tf€  PROBABILITY  DISTRIBUTION) 

FLOATN  =  FLOAT (N) 

DO  600  1=1, N 

C — 

C—  *INITIALIZE  FRACTIONAL  CUTOFF  VALUE  FOR  THIS  POINT 

FRACT ( I )  =  FLOATN  ♦  FRACT ( I ) 

C— 

C—  »IF( THIS  INDEX  BELONGS  IN  TW  'LESS*  LIST,  I.E.  THOSE 

C—  INDICES  FOR  WHICH  FRACT ( I )  IS  LESS  THAN  1.01THEN 

IF(FRACTU).GE.l.O)  GO  TO  100 

C— 

C —  «ADD  THIS  I«EX  TO  f€AD  OF  ■LESS’  LIST 

IALIASU)  =  U€AD 
LHEAD  =  I 

C— 

C—  *€LSE  (INDEX  BELONGS  IN  'MORE*  LIST) 

GO  TO  200 
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oo  oo  o  o  oo  oo  o  o  <-j  <r>  o  o  o  <->  oo  oo  oo 


CONTINUE 


100 

—  *ADD  INDEX  TO  HEAD  OF  “MORE'  LIST 

I  ALIAS!  I)  =  MEAD 
ffCAD  =  I 

—  *€W)  IF  (WHICH  LIST  TEST) 


200 

CONTINUE 

300 

♦DO  WNILEIBOTH  LISTS  ARE  NOT  EMPTY) 

CONTINUE 

IF(MHEAD.EQ.O)  GO  TO  500 

IF(DCAD.EQ.O)  GO  TO  500 

— 

♦REMOVE  NEXT  INDEX  FROM  'LESS'  LIST 

LNEXT  =  L)€AD 

LHEAD  =  IALIAS(LHEAD) 

: — 

♦SET  ALIAS  FOR  THIS  INDEX  TO  NEXT  ENTRY  IN 
lALIAS(LNEXT)  =  N4EAD 

■MORE'  LIST 

: — 

♦UPDATE  CUTOFF  VALUE  FOR  THIS  ‘MORE"  ENTRY 

FRACT(ffCAD)  =  FRACTOKAD)  -  (1.0-FRACT1  LNEXT)) 

; — 

♦  IF (THIS  INDEX  NO  LONGER  BELONGS  IN  ‘MORE’ 
IF  (FRACT<  ffCAD  ).G£.  1.0)  GO  TO  400 

LISTITHEN 

♦REMOVE  INDEX  FROM  ’MORE*  LIST  AND  ADD  IT  TO  'LESS' 
LNEXT  *  HHEAD 

MHEAD  =  IALIASIMHEAD) 

IALIAS(LNEXT)  =  LHEAD 

LHEAD  =  LNEXT 

400 

♦END  IF  (SNITCH  LISTS  TEST) 

CONTINUE 

500 

♦END  DO  (LISTS  LOOP) 

GO  TO  300 

CONTINUE 

600 

♦END  DO  (INDEX  LOOP) 

CONTINUE 

700 

c— 

RETURN 

END 

♦ADJUST  FRACT(I)  TO  SAVE  TIME  IN  MNOH  SUBROUTINE. 
DO  700  1=1, N 

FRACT(I)  =  FRACT(I)  +  (1-1) 

CONTINUE 
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SUBROUTINE  ATTRIT(F10ST,NT0FLY, iaYVC, LOSTVC, N.OST) 


C++  ATTRIT  -  SIMULATE  ATTRITION  PROCESS  DINING  A  SORTIE  PERIOD. 
Cm  ATTRIT  IS  A  FORTRAN  SUBROUTINE  WHICH  SIMULATES  THE 
C*«  EFFECTS  OF  ATTRITION  DURING  A  SORTIE.  ATTRIT  DRAWS  A  RANDOM 
Cm  SAMFlf  FROM  A  BINOMIAL  DISTRIBUTION  BASED  ON  THE  NUfflER  OF 
Cm  AIRCRAFT  FLYING  THE  SORTIE  AND  THE  PROBABILITY  OF  ATTRITION 
Cm  GIVEN  AN  AIRCRAFT  FLIES  A  StRTIE.  ATTRIT  THEN  SELECTS  THE 
Cm  RIGHTMOST  AIRCRAFT  FROM  Tt£  CURRENT  FLYABLE  AIRCRAFT  VECTOR  AS 
Cm  THE  AIRCRAFT  WHICH  WERE  ATTRITED.  THE  RIGHfflOST  ONES  ARE  SELECTED 
Cm  TO  SPEED  UP  COMPUTATION  IN  OTHER  ROUTINES. 

Cm 

Cm  INPUT  - 
Cm  PLOST 
Cm 

Cm  INPUT/OUTPUT 


PROBABILITY  THAT  AN  AIRCRAFT  ATTR1TS  GIVEN  THAT 
aiES  A  SORTIE. 


IT 


NTOFLY 

IFLYVC 


Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 
Cm 

Cm  OUTPUT 
Cm  NLOST 


LOSTVC 


-  NO.  OF  A/C  TO  ay  THIS  PERIOD. 

-  aYABLE  AIRCRAFT  STATUS  VECTOR.  INDICATES  THOSE 
AIRCRAFT  WHICH  ARE  STILL  aYABLE  DURING  THE  CURRENT 
aYING  CYCLE.  I.E.  THOSE  AIRCRAFT  WHICH  WERE  aYABLE 
AT  THE  START  OF  PREaiGHT  AND  HAVE  NOT  GROUND-ABORTED 
ATTRITED.  OR  BROKEN  THUS  FAR  IN  THE  CYCLE. 

TVE  FIRST  WORD,  iaYVC(l),  CONTAINS  THE  TOTAL 
NUMBER  OF  AIRCRAFT  STILL  aYABLE  THUS  FAR  IN 
THE  CURRENT  aYING  CYCLE.  T>€  REMAINDER  OF  THE 
ARRAY  IS  A  BIT  VECTOR  WITH  EACH  BIT  REPRESENTING 
AN  AIRCRAFT.  A  1-BIT  INDICATES  THE  AIRCRAFT  IS 
STILL  aYABLE.  NOTE  THAT  iaYVC(l)  fiLSO  INDICATES 
THE  NUMBER  OF  1-BITS  IN  THIS  BIT  VECTOR. 

-  ATTRITED  AIRCRAFT  VECTOR.  INDICATES  THOSE  AIR¬ 
CRAFT  WHICH  HAVE  ATTRITED  THUS  FAR  IN  THE  SIMULATION 
AND  NOT  BEEN  REPLACED  BY  RESERVES.  THE  FIRST  WORD, 
LOSTVCU),  CONTAINS  THE  TOTAL  NUMBER  OF  AIRCRAFT 
WHICH  HAVE  BEEN  LOST  AND  NOT  REPLACED  BY  RESERVES. 

THE  REMAINDER  OF  THE  ARRAY  IS  A  BIT  VECTOR  WITH 
EACH  BIT  REPRESENTING  AN  AIRCRAFT.  A  1-BIT  INDICATES 
TVC  AIRCRAFT  MS  BEEN  ATTRITTED.  NOTE  THAT 

LOSTVC ( 1 )  ALSO  INDICATES  THE  NUMBER  OF  1-BITS  IN 
THIS  BIT  VECTOR. 

-  NUMBER  OF  AIROWT  LOST  ON  THIS  SORTIE 


C — 

C—  *  DETERMINE  NUMBER  OF  ATTRITIONS  BY  SAMaiNG  FROM  THE 

C —  APPRECIATE  BINOMIAL  DISTRIBUTION 

M.OST  =  NBINOM(aOST.NTOaY) 

C — 

C —  *IF(ANY  AIRCRAFT  WERE  ATTRITEDJTHEN 

IFIfCOST  .E9.  0)  GO  TO  1000 

C— 

C—  +REDUCE  NO.  OF  A/C  CAPABLE  CF  aYING  THIS  PERIOD 

NToay  =  NToay  -  nlost 


BLOCK  DATA 


CHHIHIIIIIimHHUlimHIHWIInmimiHUHIHIHmHHHI 

C++  BLOCK  DATA  -  INITIALIZES  COWUN  TABLES  FOR  BIT  MANIPULATIONS. 

C++*  THIS  SUBPROGRAM  INITIALIZES  TIC  TABLES  CONTAINED  IN 
Cm  THE  /BIT/  COMMON  BLOCK.  THIS  INITIALIZATION  IS  DONE 
C***  DURING  COFILAriONi  Tf€  SUBPROGRAM  CONTAINS  NO 


Cm  EXECUTABLE  STATEICNTS.  THE  /BIT/  COMMON  BLOCK  CONTAINS 
Cm  TWEE  SETS  OF  TABLES  WHICH  ARE  USED  FOR  ACCESSING  BITS  AND 
Cm  BIT  FIELDS  WITHIN  A  COfUTER  WORD.  NOTE  THAT  THE  FOLLOWING 
Cm  PROGRAMMING  TECWIQUE  IS  USED  IN  EACH  OF  THESE  TABLES  -  AN 
Cm  EXTRA  WORD  IS  PLACED  BEFORE  TIC  BEGINNING  OF  EACH  TABLE.  THIS 
Cm  EXTRA  WORD  REPRESENTS  TABLE(O).  I.E.,  THE  OTH  INDEXED  WORD  IN 
Cm  THE  TABLE.  THUS  .  TIC  TABLE  IS  ACTUALLY  IWJEXED  0,1,2.... 

Cm  THIS  TECHNIQUE  OF  REFERENCING  THE  OTH  WORD  OF  AN  ARRAY  IS 
Cm  NOT  STANDARD  FORTRAN  AND  MAY  NOT  WORK  WITH  OTHER  FORTRAN  COMPILERS. 
Cm  TICSE  TABLES  REMAIN  FIXED  THROUGHOUT  THE  SIMULATION. 


Cm 

Cm  COWON  TABLES  - 

Cm  MASK(I)  -  1=0,1,. ..,35.  MASK  IS  TIC  BIT  ACCESSING  TABLE 
Cm  USED  IN  THE  SGM.  THE  BITS  IN  TIC  COMPUTER  WORD 

Cm  ARE  NUMBERED,  LEFT  TO  RIGHT,  0,1,2 . 35, 

Cm  AND  MSK(I)  HAS  A  1-BIT  IN  THE  ITH  POSITION  AND 

Cm  ZEROES  ELSEWHERE.  THIS  TABLE  IS  USED  TO  MASK-OFF 

Cm  T>€  ITH  BIT  IN  A  COFUTER  WORD. 

Cm  MSKLFT(I)  -  1=0,1,...  >36  .  MSKLFT  IS  USED  TO  MASK-OFF  THE 
C»*  LEFTMOST  BITS  IN  A  COHVTER  WORD.  THE  FIRST 

Cm  (LEFTMOST)  I  BITS  OF  NSKLFTd)  ARE  1-BITS 

Cm  AND  TIC  REMAINING  BITS  ARE  ZERO.  THUS,  FOR 

Cm  EXAMPLE,  MSKLFT<0)  WOULD  BE  ALL  OS  AND 

Cm*  MSKLFT  (36)  WOULD  BE  IS. 

Cm  I  COUNT  ( I )  -  1=0,1,. ..,63.  THIS  IS  A  TABLE  WHICH  IS  USED  TO 

Cm  COUNT  THE  NUMBER  OF  1-BITS  IN  ANY  GIVEN  6-BIT 

Cm  FIELD.  IN  A  6-BIT  FIELD,  THERE  ARE  2m6  =  64 

Cm  POSSIBLE  BIT  PATTERNS  -  TIC  BINARY 

Cm*  REPRESENTATIONS  OF  THE  INTEGERS  0,1, 2,... 63. 

Cm  ICOUNT(I)  CONTAINS  TIC  NUMBER  OF  1-BITS  IN  THE 

Cm  BINARY  REPRESENTATION  OF  I,  E.G.,  IC0UNT<3)=2  . 

Cm  THIS  TABLE  IS  USEFUL  IN  COUNTING  THE  NUMBER  OF 

Cm  l -BITS  REPRESENTING  AIRCRAFT  IN  TIC  VARIOUS 

Cm  AIRCRAFT-STATUS  BIT-VECTORS.  THIS  TECHNIQUE 

Cm  IS  MUCH  FASTER  THAN  A  BIT-BY-BIT  COUNT. 


C — 

COmON  /BITS/  HASK0,MASK(35),  HLEFTO, MSKLFT (36), 

4  1 2COL)T » I  COUNT  ( 63 ) 

C— 

DATA  HASKO/0400000000000  / 

DATA  MASK/  0200000000000,  0100000000000, 

4  040000000000  ,  020000000000  ,  010000000000  , 

4  04000000000  .  02000000000  ,  01000000000  , 

4  0400000000  .  0200000000  ,  0100000000 

&  040000000  ,  020000000  ,  010000000 
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S' 

04000000 

,  02000000 

,  01000000 

k 

0400000 

,  0200000 

,  0100000 

k 

040000 

,  020000 

,  010000 

k 

04000 

,  02000 

,  01000 

k 

0400 

.  0200 

,  0100 

k 

040 

,  020 

,  oto 

k 

04 

,  02 

,  01 

c- 


/ 


Ir 

k 

k 

k 

k 

k 

k 

k 

k 


DATA  rtJTTO/O/ 

DATA  (mn/0400000000000,  0600000000000.  0700000000000, 
0740000000000,  0760000000000,  0770000000000, 
0774000000000,  0776000000000,  0777000000000, 
0777400000000,  0777600000000,  0777700000000, 
0777740000000,  0777760000000,  0777770000000, 
0777774000000,  0777776000000,  0777777000000, 
0777777400000,  0777777600000,  0777777700000, 
0777777740000,  0777777760000,  0777777770000, 
0777777774000,  0777777776000,  077777777700 0, 
0777777777400,  0777777777600,  0777777777700, 
0777777777740,  0777777777760,  0777777777770, 
077777/777774.  0777777777776,  0777777777777  / 


DATA  IIC0UT/07 
DATA  I COUNT/ 
l 
% 

& 

& 

* 

Si 

C — 

END 


i>  ti  2i  It  2*  2»  3>  1>  2i  2s 

3>  2,  3,  3»  4,  i,  2,  2,  3,  2> 

3*  3*  4*  2»  3i  3*  4*  3»  4»  4> 

5»  1»  2»  2*  3»  2>  3»  3»  4»  2» 

3»  3»  4>  3»  4i  4i  5*  2?  3*  3« 

4,  3>  4,  4,  3»  4»  4,  5,  4, 

5,  5»  6  / 
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SUBROUTINE  Bfl£AKIPBREAK»P8RKSEQ> INDXWC.NTOFLY, IFLYVC.NORSVC) 


C++  BREAK  -  SIMULATE  AIRCRAFT  BREAKS  AFTER  A  SORTIE. 

Cm  BREAK  IS  A  FORTRAN  SUBROUTIIC  WHICH  SIMULATES  THE  PROCESS 
C++*  OF  AIRCRAFT  BREAKING  UPON  RETURNING  FROM  A  SORTIE.  GIVEN  A  HUGER 
C++*  OF  FLYABLE  AIRCRAFT  AND  TIC  OVERALL  BREAK  RATE,  THIS  ROUTINE  FIRST 
C++*  DETERMINES  TIC  NUWER  OF  AIRCRAFT  WHICH  BROKE.  IT  THEN  DETERMINES 
C++*  THE  HUGER  AND  DISTRIBUTION  OF  PARTS  DEMANDS  RESULTING  FROM  THESE 
C++*  BREAKS.  THEN  ASSUMING  IftCDIATE  AM)  MAXIMUM  CANNAB1LIZATIGN,  TIC 
C++*  NUMBER  OF  NORS  AIRCRAFT  AMONG  TICSE  BROKEN  AIRCRAFT  IS  DETEM1IICD. 
Cm  THOSE  AIRCRAFT  WHICH  ARE  NOT  NORS  ARE  PROBABILISTICALLY  BROKEN 
Cm  DIRECTLY  INTO  THE  VARIOUS  WORKCENTEFS. 

Cm 


Cm  I  (CUT  - 
Cm  PBREAK 
Cm 
Cm 

C++*  P8RKSE8  - 

Cm 

C++* 

Cm  INDXWC 
C++* 

Cm  INPUT/OUTPUT  - 
Cm  NTOFLY 
Cm  IFLYVC 
Cm 


NORSVC 


PROBABILITY  THAT  A  FLYABLE  AIRCRAFT  BREAKS 
INTO  AT  LEAST  ONE  WORKCENTER  UPON  RETURNING 
FROM  A  SORTIE. 

2-DIMENSIONAL  ARRAY  USED  TO  DETERMINE  TIC 
DISTRIBUTION  OF  ABORTS  INTO  TIC  VttlOUS 
WORKCENTERS. 

AN  INDEX  ARRAY  USED  TO  DETERMIIC  THE  DISTRIBUTION 
OF  BREAKS  INTO  TIC  VARIOUS  WORK  CENTERS. 


-  NO.  OF  A/C  TO  FLY  THIS  PERIOD. 

-  FLYABLE  AIRCRAFT  STATUS  VECTOR.  INDICATES  THOSE 
AIRCRAFT  WHICH  ARE  SHU.  aYABLE  DURING  THE  CURRENT 
aYING  CYCLE.  I.E.  THOSE  AIRCRAFT  WHICH  WERE  aYABLE 
AT  TIC  START  OF  PREaiGHT  AND  HAVE  NOT  GROUND-ABORTED 
ATTRITED,  OR  BROKEN  THUS  FAR  IN  THE  CYCLE. 

THE  FIRST  WORD,  laYVC(l),  CONTAINS  THE  TOTAL 
NUMBER  OF  AIRCRAFT  STILL  aYABLE  THUS  F«  IN 
THE  CURRENT  aYING  CYCLE.  THE  REMAINDER  OF  THE 
ARRAY  IS  A  BIT  VECTOR  WITH  EACH  BIT  REPRESENTING 
AN  AIRCRAFT.  A  1-BIT  IWJICATES  THE  AIRCRAFT  IS 
STILL  aYABLE.  NOTE  THAT  iaYVC(l)  ALSO  INDICATES 
THE  NUMBER  OF  1-BITS  IN  THIS  BIT  VECTOR. 

-  NORS  AIRCRAFT  STATUS  VECTOR.  IWJICATES  THOSE 
AIRCRAFT  WHICH  ARE  NORS  DUE  TO  UNAVAILABLE  PARTS. 

THE  FIRST  WORD,  NORSVC(l),  CONTAINS  TTC  TOTAL 
NUMBER  OF  1-BITS  IN  THE  NORS  STATUS  VECTOR. 

THE  REMAINDER  OF  TIC  ARRAY  IS  A  BIT  STRING  WITH 
EACH  BIT  REPRESENTING  AN  AIRCRAFT.  A  1  INDICATES 
THE  AIRCRAFT  IS  NORS.  NOTE  THAT  NORSVC(l)  fiLSO 
IWJICATES  THE  NUMBER  OF  1-BITS  IN  THIS  BIT  STRING. 


DIMENSION  NORSVCU) 

»IF( THERE  J*E  STILL  Mi  aYABLE  AIRCRAFTJTHEN 
IF(NTOaY.EO.O)  GO  TO  4000 


♦DETERMIIC  NUMBER  OF  AIRCRAFT  BREAKING  INTO  WORKCENTERS 


o  o  oo  o  o  o  oo  o  o  o  o  o  o  o  o  o  o  oo  o  r> 


C —  BY  SAILING  FROM  THE  APPROPRIATE  BINOMIAL  DISTRIBUTION 

NTOTBK  =  NBINOH(PBREAK.NTOFLY) 

—  *IF(THER£  ARE  ANY  BROKEN  AIRCRAFT)THEN 

IF(NTOTBK.EQ.O)  GO  TO  3000 

—  *REDUC£  NO.  OF  A/C  CAPABLE  OF  FLYING  THIS  PERIOD 

NTOFLY  =  NTOFIY  -  NTOTBK 

—  •DETERMINE  NUMBER/DISTRIBUTIQN  OF  PARTS  DEMANDS  RESULTING 

—  FROM  THESE  BROKEN  AIRCfWT  AW)  DETERM  I  f€  WU  NUMBER  OF 

--  NORS  AIRCRAFT  AFTER  IffOIATE  AND  MAXIMUM  C###BILIZATION 

NENNOR  =  NORSBK ( NTOTBK. NORSVC ( 1 ) ) 

NORDIF  =  WUNOR  -  NORSVC(t) 

—  *IF(N0T  ALL  THE  BROKEN  AIRCRAFT  (WE  NORS) THEN 

IFINORDIF.GE. NTOTBK)  GO  TO  1000 

—  *BREAK  T)€  LEFTMOST  FLYABLE  AIRCWWT  INTO  MAINTENANCE 

CALL  tCDIST(NTOTBK-NORDIF.PBRKSEQ. 1NDXUC. IFLYVC) 

—  •END  IF  (ALL  NORS  TEST) 

1000  CONTINUE 

—  *IF(SOME  OF  THE  BROKEN  AIRCRtfT  ARE  NORS) THEN 

IF(NORDIF.LE.O)  GO  TO  2000 

—  *  TRANSFER  LEFTMOST  AIRCRAFT  FROM  FLYABLE  STATUS 

—  VECTOR  TO  THE  NORS  STATUS  VECTOR 
CALL  TBITSL(NORDIF.IFLYVC.NORSVC) 

—  ♦END  IF  (NONZERO  NORS  TEST) 

2000  CONTINUE 

—  «END  IF  (ZERO  BREAKS  TEST) 

3000  CONTINUE 

C — 

C —  »£ND  IF  (ZERO  FLYABLE  AIRCRAFT  TEST) 

4000  CONTINUE 

C — 

RETURN 

END 
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SUBROUTINE  (^SERVUAUGMNT, LOSTVC, tEESRV, NAC) 

ClIHItHIIHIIHIIHHUlimmilllHIHIimillllHIIHlWHHHIl 

C++  CRESERV  -  OMIT  RESERVE  AIRCRAFT. 

Cm  CRESERV  IS  A  FOR^RM  SUBROUTINE  UHICH  WILL  ALLOCATE 
Cm  RESERVE  AIRCRAFT  fO  REPLACE  THOSE  AIRCRAFT  WHICH  HAVE  BEEN 
Cm  LOST  DUE  TO  ATTRITION.  IF  ENOUGH  RESERVES  ARE  LEFT.  ALL  LOSSES 
Cm  ARE  REPLACED?  FENCE  THE  ATTRITION  VECTOR  IS  ZEROED  OUT.  IF 
Cm  TJCRE  ARE  NOT  ENOUGH  RESERVES  TO  COVER  ALL  HE  LOSSES,  THEN 
Cm  THE  ATTRITIONS  ON  TIC  LEFT  OF  THE  ATTRITION  VECTOR  ARE 
Cm  REPLACED.  THIS  ARBITRARY  SELECTION  WILL  ICLP  TO  SPEED  UP  TIC 
Cm  SELECTION  ROUTINES.  NOTE  THAT  ALL  RESERVES  ARE  ASSUMED  TO 
Cm  BE  FULLY  MISSION  CAPABLE  II.E.  FLYABLE)  WHEN  COMMTTED. 


Cm 

Cm  INPUT  — 

Cm  IAUGMT 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm  INPUT/OUTPUT 
Cm  LOSTVC 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

C++* 

Cm 

Cm  NRESRV 

C++*  NAC 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 


FLAG  INDICATING  HCTVER  RESERVES  ARE  TO  BE  USED 
ONLY  AS  ATTRITION  FILLERS  OR  TO  AUGMENT  TIC 
CURRENT  UE-STRENGTH.  IF  IAUGMT=0,  RESERVES 
ARE  USED  OILY  TO  REPLACE  COMBAT  LOSSES!  HENCE 
NOT  ALL  RESERVES  HAY  BE  OMITTED  WHEN  THEY 
BECOME  AVAILABLE.  IF  IAUGMM,  ALL  RESERVES 
ARE  OMITTED  IMMEDIATELY  UPON  BECOMING 
AVAILABLE.  THIS  FLAG  IS  A  USER-SPECIFIED  INPUT 
WHICH  REMAINS  FIXED  THROUGHOUT  TIC  SIMULATION. 

ATTRITED  AIRCRAFT  VECTOR.  INDICATES  THOSE  AIR¬ 
CRAFT  WHICH  HAVE  ATTRITED  THUS  FAR  IN  THE  SIMULATION 
AND  NOT  BEEN  REPLACED  BY  RESERVES.  THE  FIRST  WOW), 
LOSTVC!  1),  CONTAINS  THE  TOTW.  NUMBER  OF  AIRCRAFT 
WHICH  HAVE  BEEN  LOST  AND  NOT  REPLACED  BY  RESERVES. 
THE  REMAINDER  OF  THE  ARRAY  IS  A  BIT  VECTOR  WITH 
EACH  BIT  REPRESENTING  AN  AIRCRAFT.  A  1-BIT  INDICATES 
THE  AIRCRAFT  HAS  BEEN  ATTRITTED.  NOTE  THAT 
LOSTVC!  1 )  ALSO  INDICATES  TIC  NUMBER  OF  i-BITS  IN 
THIS  BIT  VECTOR. 

NUMBER  OF  AIRCRAFT  CURRENTLY  IN  RESERVE. 

CURRENT  UE-STRENGTH.  IF  THE  AUGMENT  FLAG 
IS  SET,  THE  RESERVES  WHICH  REMAIN  AFTER 
REPLACING  COMBAT  LOSSES  WILL  BE  USED  TO 
AUGtCNT  THIS  CURRENT  UE-STRENGTH.  IF  THE 
FLAG  IS  NOT  SET  OR  THERE  ARE  NOT  ENOUGH 
RESERVES  TO  COVER  ALL  THE  COMBAT  LOSSES,  TNEN 
NAC  WILL  REMAIN  UNCHANGED. 


ClHWIIIHHHHHAHIHimiUHIHIHWIHUHIIHHWWWtl) 

c— 

DINENSIGN  LOSTVC! I) 


C —  +REPLACE  AS  MANY  AIRCRAFT  LOSSES  AS  POSSIBLE 
NFILLS  =  MIN0!L0STVCU>,  NRESRV) 

CALL  ZBITSLOf  ILLS.LOSTVC) 

NRESRV  *  (DESRV  -  NEILLS 

C — 

C —  +IF! THERE  ARE  STILL  REMAINING  RESERVES  AND  EXCESS  RESERVES 
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C —  ARE  TO  BE  AU«€NTED>THEN 

IF(NRESRV.LE.O)GO  TO  100 
IF<  I  AUGttfT.  £8.0)00  TO  100 

C — 

C —  ’INCREASE  UE  BY  AUGKNTING  REMAINING  RESERVES 

NAC  =  NAC  +  NRESRV 
CALL  UEUPDAT(NAC) 

C — 

C—  ♦SET  REMAINING  RESERVES  TO  NONE 
NRESRV  =  0 

C — 

C —  «£ND  IF  (AUGMENTATION  TEST) 

100  CONTINUE 
C — 

RETURN 

00 
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SUBROUTINE  FLYCYC 


C++  RVCYC  -  SIMULATE  AIRCRAFT  FLYING  CYCLE. 

Cw  THIS  ROUTItt  IS  T>€  BASIC  LOGICAL  STRUCTURE  OF  THE  SQM 
Cm*  SIMULATION.  TIC  VARIOUS  DISCRETE  EVENTS  WHICH  CAN  OCCUR, 

C«H>  GROUW-ABORTS.  AIRCRAFT  REPAIRS,  BREAKS,  ETC. ,  ARE 
Cm*  STRUCTURED  ACCORDING  TO  A  USER-SPECIFIED  FLYING  CYCLE 
C««  CONSISTING  OF  A  MINIMAL -RECOVERY  PERIOD.  A  SORTIE-PERIOD, 

Cm*  AW  El  TIER  A  WAIT  OR  AN  OVERNIGHT  PERIOD.  A  SPECIFIED  SEQUENCE 
CM*  OF  TICSE  FLYING  CYCLES  COCRISE  A  FLYING  DAY,  AW  THE 
Cm*  SIMULATION  CONSISTS  OF  A  SEQUENCE  OF  FLYING  DAYS. 

CM*  TIC  NUMEROUS  IWUTS  AW  OUTPUTS  OP  THIS  ROUT  1 1C  ARE 
CM*  ALL  CONTAIfCD  IN  COMM  BLOCKS.  DEFINITIONS  ARE  PROVIDED 
Cm*  IN  TIC  VARIOUS  ROUTIWS  CALLED  BY  FLYCYC  AW  HILL  NOT 
Cm*  BE  REPEATED  HERE. 


C — 

PAfWCTER  MAXAC=106,MAXHC=25.HAX8IT=36> 
h  MftXVEC=2+ ( MAX AC- 1 ) /HAXB I T 

PARAMETER  MAXDAY=30,HAXCYC=10,HAXSTAT=5 
LOGICAL  IICMAN,  IICPART 

COMMON  /ACSTATE/  LENGTH,  NACVC(HAXVEC),  IFLYVC(NAXVEC), 

&  MAINVC(MAXVEC),  NORSVC(MAXVEC),  LOSTVC(WXVEC) 

COftWN  /INPUT/  INITUE,  NAC,  PATTRIT,  IRES,  RWCM.  IICPART, 

I <  MAXFLY(MAXCYC),  IICMAN,  ISCALE,  IAUGMNT 

COWON  /STATS/  EXPECT !  MAXSTAT ,  HAXCYC,  MAXDAY ) , 

l  NRESRV,  IZDIAY,  ITflTRESIMAXOAYK  LOSSTOT 

COmON  /TIME/  PREFLITE,  SORTLGTH,  HAITCYC,  TYMNITE, 

I>  NSIM,  ISIM,  NUMDAY,  IDAY,  NCYCIES,  ICYCLE 

COMMON  /HCBRK/  PACBRK,  PACGABT,  PBRKUC(MAXHC).  PHCPROD. 

k  PBRKSEQ(2,HAXUC),  IWXUC(MAXUC) 

COWON  /HC INPUT/  NUC,  NCREWS(MAXHC),  SRATE(NAXHC) 

C— 

C—  HJPDATE  OVERALL  MAINTENANCE  BIT-VECTOR  FOR  THIS  FLYING  CYCLE 
CALL  MUPDATE(NHC.NAINVC) 

C— 

C —  ‘COMPUTE  NUMBER  OF  FULLY-H I  SSI  ON-CAP  ABLE  AIRCRAFT  AW  NUMBER 
C —  OF  AIRCRAFT  TO  SCHEDULE  FOR  ICXT  SORTIE 

C —  (A  MISSION-CAPABLE  AIRCRAFT  IS  DEFIICD  AS  ONE  NOT  IN 

C —  MAINTENANCE,  NORS,  OR  COIffiAT  LOSS  STATES! 

DO  IOO  L*2, LENGTH 

IFLYVC(L)  =  XOR  (NACVC(L),LOSTVC(L),MAINVC(L),NORSVC(L) > 

IOO  CONTINUE 

iaYVC(l)  =  N1VECT  (IFLYVC) 

NTOFLY  =  MINO  (HAXFLYUCYCLEMFLYVCU)) 

C — 

C—  *AIRCRAFT-REPAIR  EVENT  -  OETERMIIC  AIRCRAFT  REPAIRED  IN 
C—  EACH  HORK-CENTER  DURING  MINIMAL-RECOVERY  PERIOD 
CALL  REPAIRIPRERITE.NHC.NCREWS, SPATE! 

C— 

C —  *GROUND-ABORT  EVENT  -  OETERMIIC  NUWER  OF  GROUND  ABORTS 
CALL  GA80RT  < PACGABT  >  PBRKSEQ,  I WXHC ,  NTOFLY ,  IFLYVC ) 


OOO  C~i  e~)  OOO  OO  OO  O  C“>  O  oo  oooo  ooooo 


c — 

C —  ♦STATISTICS  EVENT  -  UPDATE  CUMULATIVE  STATISTICS 
C —  (STATISTICS  ARE  ALWAYS  COLLECTED  AT  THE  BEGIMUNG 

—  OF  THE  SORTIE  PERIOD.  TIC  MAINTENANCE  STAT 

—  REFLECTS  ONLY  THOSE  AIRCRAFT  IN  MAINTENANCE  AT  THE 

—  START  OF  THE  MINIMAL-RECOVERY  PERIOD  AND  DOES 

—  NOT  ACCOUNT  FOR  GROUND-ABORTS  OR  REPAIRS  DURING 

—  THIS  MINIMAL  RECOVERY  PERIOD) 

EXPECT ( 1 . ICYCLE. IDAY)  =  EXPECT ( 1 > I CYCLE > IDAY )  ♦  NTOFLY 
EXPECT (2i ICYCLEi IDAY)  =  EXPECT<2. ICYCLE. IDAY)  ♦  NAIWC(l) 

EXPECT(3. ICYCLE. IDAY)  =  EXPECT<3. ICYCLE. IDAY)  +  NORSVC(l) 

EXPECTW,  ICYCLE.  IDAY)  =  EXPECTM,  ICYCLE,  IDAY)  ♦  LOSSTOT 
EXPECT(5, ICYCLE, IDAY)  =  EXPECT(5. ICYCLE, IDAY)  ♦  WESRV 

—  *IF ( THIS  IS  NOT  THE  LAST  CYCLE  OF  THE  LAST  DAY (THEN 

—  (NONE  OF  THE  FOLLOWING  WORK  WILL  AFFECT 

—  THE  OUTPUT  RESULTS  IF  THIS  IS  THE  LAST  SORTIE) 

IF(  ( ICYCLE.  GE.NCYCLES).  AM).  (IDAY.  GE.NUMDAYDGO  TO  400 

—  *ATTRITION  EVENT  -  DETERM  M  NUMBER  OF  AT  TR I  TED  AIRCRAFT 

CALL  ATTRIT(PATTRIT,NTOFLY,IFLYVC,LOSTVC,NLOST) 

LOSSTOT  =  LOSSTOT  +  MOST 

•REPAIR  EVENT  -  DETERMINE  NUMBER  OF  AIRCRAFT  REPAIRED  IN  EACH 

—  WORK  CENTER  DURING  THE  SORTIE  PERIOD 
CALL  REPAIR(SORTLGTH,NUC,NCREWS,SRATE) 

—  »IF(IT  IS  NOT  TT€  LAST  SORTIE  OF  THE  FLYING  OAY)THEN 

IF< ICYCLE. EQ.NCYCLES)  GO  TO  200 

—  *  BREAK-EVENT  -  DETERMINE  AIRCRAFT  BREAKS  AND  PART  DEMANDS 

CALL  BREAK (PACBRK.PBRKSEQ, INDXWC, NTOFLY, I FLYVC.NURSVC) 

—  tAIRCRAFT -REPAIR  EVENT  -DETERMINE  AIRCRAFT  REPAIRED  IN 

—  EACH  WORK-CENTER  DURING  THE  WAIT  PERIOD 
CALL  REPAIR(WAITCYC,NWU,NCREWS,SRATE) 

—  *ELSE  (THIS  IS  TM  OVERNIGHT  PERIOD) 

GO  TO  300 
200  CONTINUE 

—  «PARTS-REPAIR  EVENT  -  DETERMINE  SPARE  PARTS  REPAIRED 

—  IN  THE  LAST  24-HOUR  PERIOD 
CALL  PRTREP(24. .PBRKSEQ, INDXWC, NORSVC) 

C — 

C —  * BREAK -EVENT  -  DETERMINE  AIRCRAFT  BREAKS  AND  PART  DEMANDS 

CALL  BREAK(PACBRK, PBRKSEQ,  INDXWC, NTOFLY,  iaYVC.NORSVC) 

C — 

C —  *A I RCRAFT -REPAIR  EVENT  -  DETERMINE  OVERNIGHT  AIRCRAFT  REPAIRS 

CALL  REPAIR(TYHNITE,NWC,NCR£WS,SRATE) 

C— 

C —  ♦COfTIIT-RE  SERVES  EVENT  -  BRING-IN  AVAILABLE  FULLY-MISSION- 

C—  CAPABLE  RESERVE  AIRCRAFT  TO  REPLACE  ANY  COMBAT  LOSSES 
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CALL  CRESERV(IAUGWT.LOSTVC.II«SRV,NAC) 


C — 

C —  too  IF  (OVERNIGHT  PERIOD  TEST) 

300  CONTINUE 
C — 

C —  *00  IF(LAST-OCt£-CF-iAST-DAY  TEST) 
400  CONTINUE 
C — 

RETURN 

END 


B-17 


SUBROUTINE  GAB0RT(PA80RT.PBRKSEQ,  INDXWC. NTORY. IRYVC) 
C»tminiiiimiiHH»mmninmmnnHinnHHmHHHHH 
C++  GABORT  -  SIMULATE  AIRCRAFT  GROUND-ABORT  PROCESS. 

C++*  GABORT  IS  A  FORTRAN  SUBROUTINE  WHICH  SIMULATES  THE  PROCESS 
C»*  OF  AIRCRAFT  GROUND-ABORTING  INTO  WORKCENTERS  AT  THE  END  OF  PRERIGHT. 
C++*  GIVEN  A  NUMBER  OF  RYA8LE  AIRCRAFT  AM)  THE  OVERHJ.  GROUND-ABORT 
C*«  RATE.  THIS  ROUTINE  CALCULATES  THE  TOTH.  NUMBER  OF  GROUND-ABORTS. 

Cm*  AND  TVEN  DETERMINES  WHICH  WORKCENTERS  THESE  AIRCRAFT  BROKE  INTO. 

Cm* 

C*m  IfEUT  - 

C*m  P ABORT  -  PROBABILITY  THAT  A  RYABLE  AIRCRAFT  GROUND-ABORTS 

Cm*  INTO  AT  LEAST  ONE  WORKCENTER  DURING  PRERIGHT. 

C*m  PBRKSES  -  2-DIMENSIONAL  ARRAY  USED  TO  DETERMINE  TIE 
Cm*  DISTRIBUTION  OF  ABORTS  INTO  THE  VARIOUS 

Cm*  WORKCENTERS. 

Cm*  INDXWC  -  HI  IffflEX  ARRAY  USED  TO  DETERMINE  THE  DISTRIBUTION 
Cm*  OF  BREAKS  INTO  THE  VARIOUS  WORK  CENTERS. 

Cm*  INPOT/OUTPUT  - 

C*m  NTORY  -  NO.  OF  A/C  TO  RY  THIS  PERIOD. 

Cm*  IRYVC  -  RYABLE  AIRCRAFT  STATUS  VECTOR.  IffflICATES  THOSE 

C*m  AIRCRAFT  WHICH  ARE  STILL  RYABLE  DURING  THE  CURRENT 

Cm*  FLYING  CYCLE.  I.E.  THOSE  AIRCRAFT  WHICH  HERE  RYABLE 

Cm*  AT  THE  START  OF  PRERIGHT  AND  HAVE  NOT  GROUND-ABORTED. 

C*m  ATTRITED.  OR  BROKEN  THUS  FAR  IN  THE  CYCLE. 

C*m  THE  FIRST  WORD.  IRYVC(i).  CONTAINS  THE  TOTH. 

CM*  NUT8ER  OF  AIRCRAFT  STILL  RYABLE  THUS  FAR  IN 

Cm*  THE  CURRENT  FLYING  CYCLE.  THE  REMAINDER  OF  THE 

CM*  ARRAY  IS  A  BIT  VECTOR  WITH  EACH  BIT  REPRESENTING 

Cm*  AN  AIRCRAFT.  A  1-BIT  INDICATES  THE  AIRCRAFT  IS 

C*m  STILL  RYABLE.  NOTE  THAT  IRYVC(I)  ALSO  IfffilCATES 

CM*  THE  NUWER  OF  1-BITS  IN  THIS  BIT  VECTOR. 

CM  l«iHIM*t*M*****M***§M*  ************  I  HM>H  Mil  IMMMMMMM* 

C— 

C —  *DETERMIIC  NUMBER  OF  AIRCRAFT  GROUND-ABORTING  INTO  WORKCENTERS 

C—  BY  SAMPLING  FROM  Tit  APPROPRIATE  BINOMIAL  DISTRIBUTION 

NTOT8K  =  NBINOMIPABORT, NTORY) 

C — 

C —  ♦ IF (THERE  ARE  ANY  BROKEN  AIRCRAFT )T1£N 

IF(NTOTBK.Efi.O)  GO  TO  1000 

C— 

C—  *REDUCE  NO.  OF  A/C  CAPABLE  OF  RYING  THIS  PERIOD 

NTORY  =  NTORY  -  NTOTBK 

C — 

C—  +  BREAK  Tit  LEFTMOST  RYABLE  AIRCRAFT  INTO  MAINTENANCE 

CALL  WCDISTINTOTBK.PBRKSEQ,  INDXWC,  IRYVC) 

C — 

C —  «€ND  IF  (ZERO  BREAKS  TEST) 

1000  CONTINUE 

C — 

RETURN 

EM) 
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SUBROUTINE  INITf IFSCEN, IFWC, IFPRT) 
CiMnmnininnnt»n«»wmHimniiminnmnwin 
C++  INIT  -  INITIALIZE  SON  SIMULATION. 

Cw  THIS  ROUTINE  READS  AND  INITIM.IZIS  T>€  VARIABLES  FOR  AN 
C++*  SOM  RUN.  IT  PERFORMS  THE  FOLLOWING  SERIES  OF  STEPS  - 
Cm*  1)  LOAD  AND  SET  TIE  VARIOUS  PARAMETERS  DESCRIBING  THE  RUN  SCENARIO 
CM*  2)  LOAD  AIRCRAFT  MAINTENANCE  MANPOWER  INPUTS.  3)  LOAD 
Cm#  SP<«E  PARTS  INFORMATION,  AND  4)  SET  MISCELLANEOUS  PARAMETERS 
Cm*  FOR  MODEL  USE.  EACH  INPUT  FILE  IS  CLOSED  IMMEDIATELY  AFTER 
Cm*  ALL  OF  ITS  INFORMATION  MS  BEEN  READ. 

Cm* 

Cm*  IMFUTS  - 

Cm*  IFSCEN  -  INPUT  FILE  CONTAINING  SCENARIO  PARAMETERS 

Cm*  IFWC  -  INPUT  FILE  CONTAINING  MAINTENANCE  NNPOUER  INPUTS 

Cm*  IFPRT  -  INPUT  FILE  CONTAINING  SPARE  PARTS  DATA 

Cm*  COMMON  INPUTS  - 

Cm*  INFPWT  -  LOGICAL  FLAG  IWJICATING  WHETfER  INFINITE  PARTS 
Cm*  ASSUMPTION  HOLDS. 

Cm*  INFMAN  -  LOGICAL  FLAG  INDICATING  WHETHER  INFINITE  (WNPOWER 
Cm*  ASSUMPTION  IS  BEING  MADE. 

Cm*  SORTLGTH  -  LENGTH  (IN  HOURS)  OF  EACH  SORTIE. 

Cm*  PACBRK  -  AIRCRAFT  BREAK  RATE. 

Cm*  NAC  -  CURRENT  UE  STRENGTH. 

Cm*  COMMON  OUTPUT  - 

Cm*  EXPECT(I.J.K)  -  CUMULATIVE  STATISTICS  ARRAY. 


C — 

PARMETER  IWXCYC=10.  MAXSTAT=5,  MAXDAY=30,  MAXWC=25 
COMMON  /INPUT/  INITUE,  WC.  PATTRIT,  IRES.  RMCM,  INFPART , 
It  HAXFLY(MAXCYC),  INFMAN,  ISCALE,  IAUGfWT 

COMMON  /STATS/  EXPECT ( MAXSTAT , MAXCYC , HAXDAY ) , 

Si  fPESRV,  I ZDAY ,  I TOTRES (HAXDAY ) ,  LOSSTOT 

COMMON  /THE/  PRERITE,  SORTLGTH,  WAITCYC,  TYMNITE. 

V  NSIH,  ISIM,  NUTCAY,  IDAY,  NLYCLES,  ICYCLE 

COHION  /WCBRK/  PACBRK,  PACGABT,  PBRKWC(HAXWC),  PWCPROD, 

l  °BRK$EQ(2,NAXWC),  INDXWC(HAXWC) 

C— 

C—  +LOAD  AND  SET  SCENARIO  INPUT  PARAMETERS 
CALL  INITSCN< IFSCEN) 

C — 

C —  *LOAD  (M  INITIALIZE  SP«E-PARTS  DATA 

CALL  INITPRT ( IFPRT , IMPPART , SORTLGTH, PACBRK ) 

C— 

C —  *READ  AIRCRAFT  MAINTENANCE  IWPOHER  INPUTS 
CALL  INI !WC( IFUC. IMPMAN, NAC) 

C— 

C —  *ZERO-OUT  CUMULATIVE-STATISTICS  ARRAY 

CALL  ZERO(EXPECT,HAXSTAT*MAXCYC*MAXDAY) 

C — 

RETURN 

EM) 
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SUBROUTINE  INITBO(NAC) 

CHHHtlimlHHHIIHHIHIIIIIIIlllmHtllHIHHWtH(HUmH 

C++  INITBO  -  INITIALIZE  PARTS  IN  RESUPPLY  AT  START  OF  SIMULATION. 
Cm  INITBO  INITIALIZES  THE  NUMBER  OF  BACKORDERS 
C++*  FOR  EACH  PART  TYPE  AT  TIC  START  OF  EACH  SIMULATION 
Cm*  REPLICATION.  FOR  EACH  PMT  TYPE.  A  RANDOM  SWLE  IS  DRAUN 
Cm  FROM  TIC  APPROPRIATE  POISSON  DISTRIBUTION  TO  DETERMINE  THE 
Cm  NUMBER  IN  RESUPPLY  AND  THE  NUMBER  OF  BACKORDERS  IS  COMPUTED 
C««  USING  THIS  RESLPPLY  NUMBER  AND  THE  INITIAL  STOCK 
Cm  LEVEL.  THE  ICAN  OF  TIC  POISSON  FOR  EACH  PART  IS  TIC 
Cm  PIPELINE  FOR  EACH  TYPE.  "NAC*  INDICATES  NUMBER  OF  ON-Htfffl 
Cm  AIRCRAFT  AT  THE  STMT  OF  THE  SIMULATION. 

Cm 

Cm  INPUTS  - 

Cm  NAC  -  CURRENT  UE-STRENGTH. 


Cm  comm  INPUTS  - 


Cm 

NPARTS 

Cm 

RESUPP(K) 

Cm 

Cm 

Cm 

Cm 

INITSJ(K) 

Cm 

C++* 

IQPA(K) 

Cm 

Bf«TS  <K) 

Cm 

Cm 

-  NUMBER  OF  PART-TYPES  BEING  MODELED. 

-  <K=1.... NPARTS)  EXPECTED  NUMBER  OF  TYPE-*  PARTS 
IN  RESUPPLY  AT  TIC  START  OF  THE  SCENARIO.  USED 

AS  THE  MEAN  OF  A  POISSON  DISTRIBUTION  TO  GENERATE 
A  SAMPLE  OF  TYPE-K  PMTS  INITIALLY  IN  RESUPPLY. 

-  (K=1,...,W>ARTS)  INITIAL  BASE  STOCK  LEVEL  OF  KTH 
PART  TYPE. 

-  <K=1.... .NPARTS)  QPA  OF  KTH  PART  TYPE. 

-  IK=1,...,W>ARTS)  BASE -NOT -REPAIRED- THIS-STATION 
RATE.  INDICATES  PROPORTION  OF  TYPE-K  FAILURES 
UHICH  ARE  REPAIRED  AT  TIC  BASE. 


Cm  COMMON  OUTPUTS  - 

<K=1,..., NPARTS)  NUMBER  OF  BACKORDERS  FOR  KTH 
PART-TYPE.  BACKORDERS  ARE  DEFIICD  AS 
(i  IN  RESUPPLY)-! INITIAL  STOCK  LEVEL) 

<K=1..... NPARTS)  NUMBER  OF  TYPE-K  PARTS  IN  BASE 
RESUPPLY. 

(K=l,...,fCAflTS)  NUMBER  OF  TYPE-K  PARTS  IN  DEPOT 
RESUPPLY. 


Cm 

NBACKO(K) 

Cm 

Cm 

Cm 

NBASE(K) 

Cm 

Cm 

NDEPOT(K) 

C**t 

»  M  III 

c— 

PARAMETER  HAXPRT=304 
COmON/RSEED/  SEED 

COIfICN  /PARTS/  ICARTS.  I  QPA  ( MAXPRT ) .  NBACKO  < MAXPRT ) . 

&  BRPRATE  (MAXPRT ) ,  DRPRATE  ( MAXPRT  ). INITSJ( MAXPRT ) .  RESUPP  (MAXPRT ) . 

8-  BNRTS  ( MAXPRT ) .  NBASE  ( MAXPRT ) .  NDEPOT  ( MAXPRT ) 

C — 

C —  +DO  FOR(EACH  PART  TYPE) 

DO  200  K= I .NPARTS 


C— 

C —  +ORAN  SAMPLE  FROM  POISSON  DISTRIBUTION  FOR  NUMBER  OF 

C —  PARTS  IN  RESUPPLY 


NRESUPP=IPOISSON(RESUPP(K) i SEED) 

C— 

C~  tCOMPUTE  INITIAL  BACKORDERS 

NBACXO(K)=)#£SUPP  -  INITSJ(K) 


i 
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o  o  o  o  «•  o  o  oo 


—  *IF  (IF  BACKORDERS  GREATER  THAN  PARTS  ON-HAND) 

IF(N8ACXO(K).L£.NAC*]QPA(IO)  GOTO  100 

—  ♦PRINT  WARNING  MESSAGE  AND  TRUNCATE  NBACKOOO 

WRITER,  9001) 

K,  NBACKOOO,  NAC,  IQPA<K)»Mtt.SUPP,  INITSJ(K),R£SUPP(K) 
NBACKO<K)=NAC*IQPA(K) 

NR£SUPP=NBACKO(K)+INITSJ(K) 

—  *END  IF( TRUNCATE  BACKORDERS  AT  IWXIHUH  AVAILABLE) 

100  CONTINUE 

—  ♦ALLOCATE  THESE  PARTS  BETWEEN  BASE  AND  DEPOT  RESUPPLY 

RBASE=0.0 

IF(8RPRATE(K).GT.0.0)  RBAS£=( l-BMTTS(K) ) /BRPRATE(K) 
RDEP0T=0.0 

IF(DRPRATE(K).GT.O.O)  RDEP0T=8NRTS(K)/DRPRATE(K) 
NDEP=NBINOf« RDEPOT/ ( RBASE+RDEPOT) ,  MIESUPP ) 
NDEPOT(K)=NDEPOT(K)+NDEP 
NBASE(K)=NBAS£(K)+NRESUPP-NDEP 

C — 

C —  ♦END  DO  (PARTS  LOOP) 

200  CONTINUE 
C— 

RETURN 

9001  FORMAT! ’(>*$*$$***  INITBQ  ERROR  -  TOO  MANY  PARTS  IN  RESUPPLY’,/, 
t>  "  ««<«<«««  K=’,I3, *  NBACK0(K)=*.I5,’  NAC=’,I3,*  IQPA(K)=’,I3, 

i  /,’  ***$*«$  NRESUPP,  INITSJ(K),  RESUPP(K)  =  ’,2I5,F10.3) 

END 


i 

I 


A 


D 

SUBROUTINE  INITPRT ( IFILE. I1FPART.S0RTLGTH.PAC8RK) 
CwmiiiiiHmnmunmmmwiimummiHHimiH 
C++  INITPRT  -  LOAD  A®  INITIALIZE  SPWC-PARTS  DATA. 

C++*  THIS  ROUTINE  LOADS  TIC  SPARE-PARTS  IfCUT  DATA  AND 
Cm  INITIALIZES  THE  STATISTICS  AND  TABLES  (CEDED  FOR 
C++*  SAMPLING  TOTAL  P<WT  DEJWfflS  HQ  dSO  DETERMINING  PART-TYPE 
C«+  FOR  A  GIVEN  BROKEN  PART. 


K  Ti  TT *  *  I  * 


c — 

PARAMETER  MAXAC=108, MAXBIT=3A, MAXVEC=2+<MWXAC-1 ) /MAXBIT 
PARAMETER  MAXPRT=304 
LOGICAL  INFPART 
CHARACTER  CNSN+18 

COMMON  /PARTS/  fCARTS . I QPA ( HAXPRT ) . NBACKO ( HAXPRT ) . 

St  BRPR ATE ( MAXPRT ) .  CRPRATE (MAXPRT ) » INITSJ< HAXPRT ) . RESUPP (HAXPRT ) . 

St  BNRTS(HAXPRT)>NBASE(MAXPRT) iNDEPOT(MAXPRT) 

COWION  /ALIASC/  FRACTiMAXPRT). IALIAS(MAXPRT).FPARTS 
COMMON  /DEMAND/  ACMEAN,  ACVAR.  NPERAC 

C— 

C—  *IF< INFINITE  PARTS  ARE  NOT  ASSUCD,  I.E.  NORS  AIRCRAFT  ARE 
C —  TO  BE  MODELED! THEN 

I F  <  I NFPART )  GO  TO  900 

C— 

C —  +READ-IN  PARTS  DATA  AND  PERFORM  ERROR  CHECKS 

NPARTS=1 
100  CONTINUE 

READI  IFILE, END=200)  CNSN.FRACT(NPARTS) .  J CPA  WARTS)  .FAP. 

St  INITSJ  WARTS).  RESUPP  WARTS)  .BNRTS  WARTS),  BDAYS.DDAYS 

IF( (FRACT(NPARTS).GT.O.Ol.AND. (FRACT(NPARTS).LE.l.O)) 

St  GO  TO  50 

IF(  IQPA  WARTS )  .GT.  0)8)  TO  50 
IF( (FAP.GT.O.O).AND. (FAP.LE. 1.0) >G0  TO  50 
IFUNITSJ  WARTS).  GE.OIGO  TO  50 
IF(RESUPP(NPARTS).GE.O.O)GO  TO  50 
IF(  (BNRTSWARTS).GE.O.O).  AND.  (BNRTS  WARTS).  LE.  1.0)) 
i  GO  TO  50 

IF(BDAYS.GE.O.O)GO  TO  50 
IF(DDAYS.GE.O.O)GO  TO  50 

WRITER,  9003)  CNSN.FRACT  WARTS), I8PA(NPARTS).FAP, 

It  INITSJWARTS) ,  RESUPP(NPARTS) ,  BNRTS  WARTS) ,  BDAYS,  DOAYS 

GO  TO  100 

50  CONTINUE 

BRPRATE ( NPARTS ) =0. 0 

IF(BDAYS.GT.O.O)  BRPRATEWARTS)=1.0/(24.0*BDAYS> 

DRPRATE 1  fPARTS ) =0. 0 

IF ( DDAYS. GT. 0. 0)  CRPRATE ( NPARTS )=1 . 0/ ( 24.0+DDAYS ) 

FRACT  <  1CARTS )  =FRACT  WARTS  >  ♦Fff’tSORTLGTH 
IF(FRACT(NPARTS).LE.O.O)GO  TO  100 
NPARTSWARTS+l 
IF(NPARTS.LE.MAXPRT)GO  TO  100 
READ! IFILE, END=200)CNSN 
WRITE(6, 9004 (MAXPRT 


oo  o  o  o  o  o  <r>  oo  o  o  o  <r>  cp 


200  CONTINUE 

NPARTS=NPARTS-1 
WRITE(6. 9005)NPARTS 

—  CLOSE-OUT  SPARES  IffVT  FILE 

CALL  FCLOSE(IFILE) 

—  COMPUTE  MEAN  MB  VARIANCE  OF  RMDOH  VARIABLE  -  TOTAL 

—  -PART-DEMANDS  PER  BROKEN  AIRCRAFT 

CALL  PSTATiPACBRK.NPARTS.IQPA.FRACT.ACMEAN.ACVAR.fPERAC) 

—  CONVERT  PART-DEMANDS-PER-FLYING-HOUR  TO  A  PDF 

CALL  MAKEPD(NPARTS.IQPA.FRACT) 

—  *SET-UP  TABLES  NEEDED  FOR  ALIAS  METHOD  OF  SASHING  PART-DEMANDS 

CALL  ALIAS1NPARTS.FRACT. IALIAS) 

FPARTS  =  FLOAT (N5  ARTS) 

—  *ELS£  (INFINITE  PARTS  ASSUMED) 

GO  TO  950 
900  CONTINUE 

—  *PRINT  MESSAGE  INDICATING  INFINITE  SPARE  PARTS  ASSUTPTIGN 

WRITE(6.9002) 

—  CND  IF  (INFINITE  SPARE-PARTS  TEST) 

950  CONTINUE 

C — 

RETURN 

9002  FORMAT*  1H0,7X,"HFINITE  SPME  PARTS  ASSUMED  FOR  THIS  SOM  RUN,  *, 

!>  /,-I.E. .  NO  AIRCRAFT  EVER  WAITS  FOR  A  SPARE  PART.*) 

9003  FORMAT COM******  INITPRT  ERROR  -  INVALID  PMT  CHARACTERISTIC*. 


S' 

/, 

•  <$«$*$$* 

NPARTS.  CNSN  =  M3.IX.AI8, 

& 

/, 

*  )<»»» 

FRACT,  IQPA,  FAP  =  *.F6.4, I4.F5.2, 

S' 

/, 

"  *«**«*«« 

INITSJ.  RESUPP.  BNRTS  =  M4.2F8.3, 

S' 

/, 

*  tttttm 

BDAYS.  DDAYS  =  *,2F10.2> 

9004  FORMAT (*<»$**$*«  INITPRT  ERROR  -  TOO  MANY  LRU  TYPES"./. 
S,  *  fttttttt  MAXPRT  =  ".15) 

9005  FORMAT* 1H0,3X,*LRU  TYPES  -  *, 14) 

EM) 
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SUBROUTIIC  INITREP 

CHHIU[|HmHHimmilllllll[HII(imilllHHHHHmHHIHIl 

C++  INITREP  -  INITIALIZE  VARIABLES  FOR  A  SIMULATION  REPLICATION. 
C**+  THIS  ROUTINE  PERFORMS  THE  LENGTHY  INITIALIZATION  NEEDED 
C***  EACH  SIMULATION  REPLICATION  OF  THE  SGM.  THIS  PROCESS 
C++*  IS  ORGANIZED  IN  THE  FOLLOWING  MWNER.  FIRST,  MISCELLANEOUS 
C+*+  OPERATIONS  ARE  PERFORfED,  THEN  SPWES  INITIALIZATION,  AND 
C+++  FINALLY,  WORK  CENTER  INITIALIZATION.  THE  NUMEROUS  COMMON  OUTPUTS 
C+++  OF  THIS  ROUTINE  ARE  NOT  DEFINED  HERE,  BUT  THE 


C**+  OPERATIONS  BEING  PERFORMED  SHOULD  BE  CLEAR  FROM  THE  PROGRM- 
C***  DESIGN  LANGUAGE  (PDL)  CORRESPONDING  TO  EACH  OPERATION. 

CHIIHHH»HimHHI[HIHHHIWI»mil|[IIHHHIimumill 

c— 


& 


& 

k 

k 

& 

k 

k 

k 


PARAMETER  MAXAC=108,MAXWC=25,lttXBIT=36,fttXPRT=304, 


MAXVEC=2+(MAXAC-l)/fWXBIT 

PARAMETER  LFU)=7, W5ERWRD=MAXBIT/LFLD,MXINWC=1+(MAXAC-I l/NPERWRD 

PARAMETER  MAXDAY=30,fttXCYC=10,MAXSTAT=5 

COMMON  /ACSTATE/  LENGTH,  MCVCUWXVEC),  IFLYVCOttXVEC), 

MAINVC(MAXVEC),  NORSVC(IWXVEC),  LOSTVC(MAXVEC) 
COMMON  /INPUT/  INITUE,  FMC,  PATTRIT,  IRES,  RWCM,  INFPART, 
MAXFLY(MAXCYC),  INFMAN,  I SCALE,  IAUGftiT 


COMMON  /PARTS/  H^TS,  ItfA(MAXPRT),  NBACKO(MAXPRT), 

BRPRATE(MAXPRT),  DRPRATE(MAXPRT),  INITSJ(MAXPRT), 
RESUPP(MAXPRT),  BNRTS(HAXPRT),  NBASEfMAXPRT', 
NDEPOT(MAXPRT) 

CCrtWN  /STATS/  EXPECT ( MAXSTAT , MAXCYC , MAXDAY ) , 

fUESRV,  I ZDAY , ITOTRES ( MAXDAY ) , LOSSTOT 
COrtlON  /UC8RK/  PACBRK,  PACGABT,  PBRKWC(fWXWC),  PWCPROD, 
P8RKSEQ(2,HAXWC)>  INDXUC(MAXWC) 

COMMON  /UCINPUT/  NWC,  NCREWS(HAXWC),  SRATEOWXWC) 

COMMON  /WCMAINT/  LISTRP<MXINWC,fWXWC),  INREPR(HAXWC) 

LOGICAL  INFPART 


•INITIALIZE  RESERVE  AIRCRAFT  COUNTS 
CALL  ZERO! ITOTRES, MAXDAY) 

IZDAY=0 

NRESRV=0 


C — 

C—  •RESET  INITIAL  UE  FOR  THIS  REPLICATION 
NAC  =  INITUE 
CALL  UEUPDAT(NAC) 

C— 

C —  •INITIALIZE  CUMULATIVE  AIRCRAFT  LOSSES  TO  NONE 
LOSSTOT  =  0 

C — 

C —  ♦CLEW  AIRCRAFT-STATUS  BIT-VECTORS 
CALL  ZERO(LOSTVC,MAXVEC. 
k  NORSVC.MAXVEC. 

I,  IWEPR.mXWC) 

C— 

C—  »SET  FIRST  m  BITS  OF  UE  FLYING  BIT-VECTOR  TO  FLYABLE 


DO  100  1=1. LENGTH 
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<r>  o  no  o  o  <->  o  oo  oo 


100  IFLYVC(1)=NACVC(I) 

—  *IF ( INFINITE  PARTS  MOT  ASSUMED  )T>€N 

IF<  I^PARTIGO  TO  200 

—  CLEAR  BASE  AM)  DEPOT  RESUPPLY  COUNTS 

CALL  ZERO ( NBASE . NAXPRT ,  NDEPOT » NAXPRT ) 

—  CALCULATE  INITIAL  BACKORDERS  FOR  EACH  PfKT-TYPE 
CALL  INITBO(NAC) 

—  ‘INITIALIZE  NUMBER/DISTRIBUTION  OF  NORS  AIRCRAFT 

NORS  -  NORSACOFARTS.IfiPA.NBACKO) 

CALL  TBITSL(NORS,IFLYVC,NORSVC) 

—  CND  IF  (INFINITE  SPARES  TEST) 

200  CONTINUE 

—  ‘INITIALIZE  NUHBER/DISTR1BUTI0N  OF  AIRCRAFT  IN  MAINTENANCE 

N8RKAC  =  INT(RNHCM*FLOAT(NAC)) 

CALL  UCDIST(NBRKAC,PBRKSEQ,INDXUC,IFLYVC) 

RETURN 


i 
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SUBROUTINE  INITSCN(IFSCEN) 

ClllHmHHHHWWHmHilimHHtHHH»WI«HIHIIIIIII»IH 

C++  INITSCN  -  READ  AND  INITIALIZES  SCEWWIO  INPUTS. 

Cm*  THIS  ROUTINE  LOADS  TIC  SCEMRIO  PARAMETERS  SPECIFIED 
CM*  BY  TIC  USER.  IT  ALSO  PREPARES  TMC  SCRATCH  FILE  (FILE  03) 

Cm*  WHICH  IS  USED  TO  WRITE  A  COPY  OF  THOSE  SCENARIO 
Cm*  PARAMETERS  WHICH  ARE  M10WED  TO  VMY  ON  A  DAILY  BASIS 
Cm*  THROUGHOUT  THE  SCENMIO.  A  LIST  OF  VALUES  IS  WRITTEN  TO  THIS 
CM*  SCRATCH  FILE  FOR  EACH  SIMULATION  DAY.  Tt£N.  WHEN  EACH 
Cm*  FLYING  DAY  BEGINS  (FOR  EACH  SIMULATION  REPLICATION),  TIC 
Cm*  PARtfCTER  VALUES  FOR  THAT  DAY  ARE  LOADED. 

ClHHIIIIIIIHIIIHWtHtHlllimHIIHIHIIIHHmHHIIIIIIHIII 

C — 

PARAMETER  MXM>25,  MAXV(WY=5,  HAXCYC=10 
CHARACTERS  FTOTYM.LTOTYM 
CHARACTER*20  CHSEED 
CHARACTER+80  NEXTLINE 

LOGICAL  VARY ( MAXVARY ) ,  VARYSU,  IICPART ,  IICMAN 
COMMON  /RSEED/  SEED 
COMMON  /TIIC/  PREFLITE,  SORTLGTH,WAITCYC, 
i  TYIWITE.NSIH. ISIM.NUMDAY, IDAV.NCYCLES, ICYCLE 

COMMON  /IMVT/  INITUE,NAC,PATTRIT,  IRES.RfiCM.INFPART, 
t,  MAXFLY(HAXCYC),  INFMAN,  ISCALE,  IAUGWT 

COTTON  /UCBRK/  PACBRK,  PACGABT,  PBRKWC(TMXHC),  PWCPROD, 
k  PfiRKSEQ(2,MAXWC),  IKDXUC(HAXUC) 

C — 

C —  *READ  STORED  IMMJT 
5  FORMAT(V) 

READ(IFSC£N,5)  (VARY<I),I=l,lttXVARY) 

READ(  IFSCEN.5)  CHSEED,  FTOTYM.LTOTYH,  INFMAN,  IICPART, NSIM,NAC 
i  .PACBRK, RNMCH.NUTfflAY, PREFLITE, SORTLGTH 

READ< IFSCEN.5)  ISCALE 
NEXTLIIC-  ' 

INITUE=NAC 

WRITE(06,9000)  NSIH, CHSEED, NAC 
VARYSW=.F. 

READUFSCEN.5)  PATTRIT, PACGABT, NCYCLES,  IRES, WAITCYC.TYWITE 
READ! IFSCEN.5)  (MAXFLYI I > , I = 1 , ICYCLES ) 

IF  ( VARY  ( 4 ) )  GOTO  100 
WRITE (06, 9008)  IRES 
GOTO  200 
100  CONTINUE 

WRITE(06.9007) 

CALL  C0NCAT(ICXTLINE,28,  .'VARY  BY  DAY', 1,11) 

VARYSU.  T. 

200  CONTINUE 

IF  <VARY<5) )  GOTO  300 

WRITEI06.9010)  MAXFLY(l) 

GOTO  400 
300  CONTINUE 

WRITE(06,9009) 

VARYSU*. T. 
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CALL  CONCATOCXTLINE, 45, 'VARIES  BY  CYCUE/DAYM,19> 

400  CONTINUE 

WRITE<06»9011)  NEXTLINE 
IF  (VARY(3) )  GOTO  500 

WRITE<06,9012)  NUHDAY, NCYCLES,  FTOTYN.LTOTYH, PREFLITE. SGRTLGTH, 
&  HAITCYC,TYHinE 

GOTO  600 
500  CONTINUE 

HRITE(06,9013)  NUW»Y,FTOTYN,LTOTYH,PREFLITE,SORTLGTH 
VARYSW=.T. 

600  CONTINUE 

URITE(06.9014)  RtfOl.PACBRK 
NEXTUI€='  ' 

IF  ( VARY ( 1  > )  GOTO  700 
NRITE(06,9015)  PATTRIT 
GOTO  800 
700  CONTINUE 

WRITE(06,9016) 

VARYSJKT. 

CALL  C0NCAT(fEXTLI)€>19,'BY  DAY',1,6) 

800  CONTINUE 

IF  (VARY12))  GOTO  900 

HRITE(06,9017)  PACGABT,  )€XTL I NE 
GOTO  1000 
900  CONTINUE 

CALL  CONCAT (NEXTLINE) 46, 'BY  DAY',1,6) 

WRITE106.9018)  fEXTLIIt 
VARYSM=.T. 

1000  CONTINUE 

DO  1200  I DAY=1, NUHDAY 

WRITEI03)  PATTRIT . PACGABT,  NCYCLES, IRES, WAITCYC, TYHNITE 
WRITE(03)  (HAXFLY<J),J=1, NCYCLES) 

IF  (.NOT.VARYSW)  GOTO  1100 
URlTE(06i9001)  'DAY  =', IDAY 
IF  (VARY(D)  WRITE(06,9002)  PATTRIT 
IF  ( VARY ( 2 ) )  URITE(06i9003)  PACGABT 
IF  (VARY13))  WRITE (06. 9004)  NCYCLES, WA I TCYC , TYNNI TE 
IF  { VARY ( 4 > )  WRITE (06. 9005)  IRES 
IF  (VARY<5>)  WRITEI06.9006)  'CYCLE', (J,J=1, NCYCLES) 

IF  (VARY(5) )  WRITE(06,9006)  'HAX-FLY', 

4  (MAXFLY(J),,M>  NCYCLES) 

1100  CONTINUE 

IF  (IDAY.EQ.NUMJAY)  GOTO  1200 

READ(IFSC£N,5)  PATTRIT . PACGABT , NCYCLES. IRES, WAITCYC. TYNNITE 
READUFSCEN.5)  (NAXFLY(J),J=1, NCYCLES) 

1200  CONTINUE 
C — 

C—  CONVERT  USER  SEED  TO  A  REAL  NUtOER 
DECODE (CHSEED, 5)  SEED 

C— 

C —  CLOSE  SCENARIO  INPUT  FILE 
C —  CALL  FCLOSE ( IFSCEN ) 

C— 
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RETURN 

9001  FORMAT! 'O', A5, 13) 

9002  FORMAT!'  ATTRITION  RATE  =',F6. 4) 

9003  FORMAT!'  GROUND-ABORT  RATE  =',F6.4) 

9004  FORMAT!'  HAVES  PER  DAY  =',I3/'  WAIT  TIME  =',F4.2, 

4  /'  0VERN1TE  RECOVERY  =',F5.2> 

9005  FORMAT! '  RESERVES  =',I3) 

9006  FORMAT! '  ',A7,10!2X,I3)( 

C 

9000  F0RMAT('r///4X'«mn»miimiiimin  minim 


4  '*******',///, n 

4  'SIMULATION-  REPLICATIONS  =',I4,5X, 

4  'RANDOM  NUMBER  SEED  =  ',A8/ 'O', 3X,' AIRCRAFT  -',3X, 

4  'UE  «M3) 

9007  FORMAT! '*',28X,  RESERVES') 

9008  F0RlttT<'*',27X, 'RESERVES  =',I3> 

9009  FORMAT ! '+' ,44X> 'MAXIMUM  LAUNCH-SIZE') 

9010  FORMAT! '♦',44X, 'MAXIMUM  LAUNCH-SIZE  =',I3> 

9011  FORMAT!'  ',A80/'0',3X,'FLYING  SCHEDULE  -'/'O', 10X, 'HAVES', 3X, 

4  'TAKEOFF', 

4  '  TIMES', 6X, 'MINIMAL', 3X, 'SORTIE', 2X,'HAIT',2X, 'OVERNIGHT'/ 

4  4X, 'DAYS', 2X,'PB?  DAY',3X, 'FIRST', 3X, 'LAST', 4X, 

4  'TUW1AR0UND' ,  2X ,  'LENGTH' ,  2X>  'TIfC' ,  2X,  'RECOVERY' ) 

9012  F0RMAT!'0',3X,I3,5X,I2,6X,A4,4X,A4,6X,F5.2,5X,F5.2,3X,F4.2, 

4  3X.F5.2) 

9013  FORMAT! 'O', 3X, 13, 4X, 'VARY',5X,A4,4X,A4,6X,F5.2,5X,F5.2,2X, 

4  'VARIES', 2X, 'VARIES'/'  ',9X,'BY  DAY',39X,'BY  DAY',2X,'BY  DAY') 

9014  FORJWT! 'O', 3X, 'RATES  -'/'0',6X> 'INITIAL', 17X, 'AIRCRAFT'/ 

4  6X,'NMCM  RATE',3X, 'ATTRITION', 3X, 'BREAK  RATE',3X, 

4  'GROUND-ABORT '/'O', 7X.F5. 3, 7X,13X,F6. 4, ) 

9015  FORMAT! '+', 17X.F6.2) 

9016  FORMAT! '♦',18X, 'VARIES') 

9017  FORMAT! '+',43X,F6. 4/'  ',A80) 

9018  FORMAT! '+',45X, 'VARIES'/'  ',A80) 

END 
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SUBROUTINE  INITWC(IFILE,INFMAN,NAC) 

CllHHmiHtHIIHIIHIH)HmHII»»HIHIIIHtlHIHHHHHHHI 

C++  INITUC  -  LOAD  AM)  INITIM.IZE  MAINTENANCE  WORK  CENTER  DATA. 

Cm  THIS  ROUTINE  INITIALIZES  THE  INFORMATION  NEEDED 
Cm  FOR  MODELING  THE  AIRCRAFT  MAINTENANCE  WORK-CENTERS. 

Cm  IT  READS  D£  MAINTENANCE  (WHFONER  INPUT  FILE.  PRINTS 
Cm  A  LISTING  OF  THESE  INPUTS,  AM)  COMPUTES  ADJUSTED  BREAK  - 
Cm  RATE  ARRAYS  FOR  WORK-CENTER  BREAKS 
Cm 

Cm  INPUTS  - 
Cm  IFILE 

Cm 
Cm 
Cm 

Cm  INFMAN 

Cm 
Cm 
Cm 
Cm 

Cm  NAC 

Cm 


-  UNIT  NUMBER  OF  THE  UFVT  FILE  FROM 
WHICH  DC  WORK-CENTER  INPUTS  WE  READ. 

THIS  FILE  IS  CLOSED-OUT  AFTER  THE  IfCUTS 
ARE  READ 

-  LOGICAL  VARIABLE  INDICATING  WFtTHER  INFINITE 
MANPOWER  IS  ASSUMED  FOR  ft!  WORK-CENTERS.  IF 
INFMAN=TRUE  DO  NUMBER  OF  SERVERS  FOR  EACH  WC 
IS  SET  EQUAL  TO  DC  FBXIMUM  ALLOWABLE  NUMBER 
OF  AIRCRAFT. 

-  UE  (UNIT  EQUIPMENT);  NUMBER  OF  AIRCRAFT 
POSSESSED  BY  THE  BASE  OF  INTEREST. 


Cm  COMMON  INPUTS  - 

Cm  PACBRK  -  USER-IfPJT  AIRCRAFT  BREAK-RATE.  PROBABILITY 
Cm  THAT  AN  AIRCRAFT  RETURNING  FROM  A  SORTIE  REQUIRES 


Cm 

Cm 

Cm 

Cm 

Cm* 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 


UNSCHEDULED  MAINTENANCE  IN  AT  LEAST  1  WORK-CENTER. 

PACGABT  -  PROBABILITY  THAT  AN  AIRCRAFT  GROUND-ABORTS  DURING 
THE  PRE-TAKEOFF  PERIOD. 

COWON  OUTPUTS  - 

NNC  -  NUMBER  OF  WORK-CENTERS  TO  BE  MODELED 

NCREWS  -  NUMBER-OF-SERVERS  ARRAY.  NCREWS(I)  IS  THE 
NUMBER  OF  SERVERS  IN  THE  ITH  WORK-CENTER. 

SRATE  -  SERVICE-RATES  ARRAY.  SRATE(I)  IS  THE  SERVICE- 
RATE  (IN  AIRCRAFT  PER  HOUR)  OF  THE  SERVERS 
FOR  THE  ITH  WORK-CENTER. 

PBRKSEQ  -  WORK-CENTER  BREAK  ARRAYS  FOR  THE  SEQUENTIAL 

SAFFLINB  PROCESS  OF  DETERMINING  WHICH  WORKCENTERS 
AIRCRAFT  BREAK  INTO. 

INDXWC  -  SORTED  ARRAY  OF  WORK-CENTER  INDICES  USED  WITH 
SEQUENTIAL-SAM1ING  PROCESS. 


C»*HimHimnmniwinnmmtwimnninnminnnuiiii 

C— 

PARAMETER  NAXWC=25,  MAXAC=108,  LFLD=7 

COfflON  /WCINPUT/  NWC,  NCREWS (MAXWC),  SRATE (MAXWC) 

COFflON  /WCBRK/  PACBRK.  PACGABT,  PBRKWC( MAXWC),  PWCPROD, 
fc  PBRKSEQ(2,  MAXWC),  1NDXWCIMAXWC) 

LOGICAL  IMMW 

C— 

C—  *READ  AM)  ECHO-PRINT  MAINTENANCE  MAFfOWER  IlffUT  FILE 
CALL  WCREAD(  IFILE, MAXWC, NWC, PBRKWC.NUSWS, SRATE) 


C— 


C—  *IF  INFINITE  HWOWER  ASSUMED  -  RESET  NUMBER  OF  SERVERS 
C —  PER  SHIFT  TO  MAX  NUMBER  OF  AIRCRAFT;  THUS  NO  AIRCRAFT 


i 
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C—  WILL  EVER  MIT  FOR  A  SERVER 

IF < INFMAN)  CALL  SPRAY(MAXAC,NCREHS,NUC> 
IF(WHAN)WRITE(6,9001) 

C— 

C—  ♦INITIALIZE  WORK-CENTER  BREAK  ARRAYS  FOR  GROUND-ABORTS  AND 
C—  BREAKS?  NOTE  THAT  TIC  SANE  MRAYS  ARE  CURRENTLY  USED  FOR  BOTH 
CALL  WCPRQB(NWC>PBRKWC,PBRKSEQ, INDXWC.PWCPROD) 

C — 

C —  iPERFORM  ERROR  CHECK  TO  ENSURE  LFU)  PARAMETER  LARGE  ENOUGH 
C —  SO  TMT  THE  BIT-FIELD  CAN  STORE  THE  MAX  AC  « 

IF<HAXAC.GT.2«LFLD)MITE(6,9002)LFLD,HAXAC 

C— 

RETURN 

9001  FORMAT  (1H0»7X,*WIN1TE  fWPOWER  ASSUICD  FOR  THIS  SOM  RUN 

i  7X,*I.E.,  THERE  ARE  NEVER  ANY  AIRCRAFT  QUEUES  IN  MAINTENANCE.") 

9002  FORMAT! ’Ott******  INITWC  ERROR  -  LFLD  PARAItTER  TOO  SMALL’,/, 

k  ■  m»»l  LFLD,  MAX  AC  =  *,215) 

END 
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C++  1  POISSON  -  GENERATE  RANDOM  SWPLE  FROM  A  POISSON  DISTRIBUTION. 

C+*+  THIS  ROUTINE  GENERATES  A  RANDOM  SAfPLE  FROM  A 
Cm  POISSON  DISTRIBUTION  WITH  A  GIVEN  MEAN.  T>£  EXPONENTIAL-DRAW 
Cm  fETHOO  IS  USED  FOR  DISTRIBUTIONS  WITH  SHALL  fEANS.  AM) 

Cm  A  NORMAL  APPROXIMATION  IS  USED  FOR  LARGER  (CANS  (  >20). 

Cm 

Cm  IfEUT  - 

Cm  WEAN  -  MEAN  OF  POISSON  DISTRIBUTION  FROM  WHICH  SAMPLE 
Cm  IS  TO  BE  GENERATED. 

Cm  KPUT/OUTPUT  - 

Cm  SEED  -  SEED  OF  RANDOM  NUMBER  GENERATOR. 


C — 

C —  »IF( INPUT  PARAMETER  IS  A  LEGITIMATE  MEAN  FOR  A  POISSON) 

IFIRMEAN.LT. 0.0)G0  TO  400 

C — 

C—  *IF(MEAN  IS  NOT  TOO  LARGE) 

1FIRMEAN  .GT.  20.0)  GO  TO  200 

C — 

C —  «USE  EXPONENTIAL  DRAW  fETHOD  FOR  POISSON  SAfPLE 

IP0ISSC#=-1 
PROD=i.O 
TEST=£XP(-fifEAN) 

100  CONTINUE 

IP0ISSGN=IP0ISSWH1 

PR0D=PR0D*UNIFH1(SEED) 

IF(PROO.GE.TEST)GO  TO  100 

C— 

C —  *ELSE  (LARGE  fEAN) 

GO  TO  300 
200  CONTINUE 

C — 

C —  *OSE  NORMAL  APPROXIfWTION  TO  POISSON 

IPOISSON=MAXO(Oi  INT(XN0RM(RtEAN,S8RT(RMEAN),SEED)+.5)) 

C— 

C—  «€ND  IF  (SIZE  OF  MEAN  TEST) 

300  CONTINUE 

C— 

C —  *ELSE  (MEAN  IS  LESS  THAN  ZERO) 

GO  TO  500 
400  CONTINUE 

C— 

C—  »SET  RETURN  VALUE  TO  ZERO  AND  PRINT  ERROR  MESSAGE 

IPOISSON  =  0 
WRITE(6i9001)RtEAN 

C— 

C— 

500 
C— 


RETURN 


*END  IF  (LEGITIMATE  !€AN  TEST) 
CONTINUE 


r 


9001  FORMAT IPOISSQN  ERROR  -  NEGATIVE  MEM 
k  *  »«««  RMEAN  =  "iFlO.5) 

END 
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INTEGER  RNCTION  LBITSIIWURD.NBITS) 


CiinmiiiwHiimwmHinnHHni[ininmnnm«itntHm 

C++  LBITS  -  HASK-OFF  LEFTMOST  1-BITS  IN  A  COMPUTER  WORD. 

C*«  LBITS  IS  A  FORTRAN  FUNCTION  WHICH  WILL  SELECT  A  GIVEN 
Cm*  NUMBER  OF  1-BITS  FRCH  THE  LEFTMOST  PORTION  OF  A  GIVEN  INPUT 
Cm*  WORD.  LBITS  RETURNS  A  WORD  CONSISTING  OF  THESE  SELECTED 
Cm*  1-BITS  WITH  0'S  EVERYWHERE  ELSE.  NOTE  THAT  THE  INPUT 
Cm*  WORD  SHOULD  CONTAIN  AT  LEAST  AS  MANY  1-BITS  AS  THE  NUMBER  TO  BE 
Cm*  SELECTED,  (BITS'  . 

Cm*  THIS  ROUTINE  IS  SPECIFIC  TO  A  COUNTER  WITH  36-BIT  WCRDS 
C*m  SINCE  IT  WORKS  BY  EXTRACTING  6-BIT  FIELDS. 

Cm* 


C*M 

Cm* 

Cm* 

Cm* 

C*M 

Cm* 

C*M 

Cm* 

C*** 

C*** 

Cm* 

Cm* 

C*** 

C*** 

Cm* 

C*M 


i  (furs  - 

I WORD  -  WORD  FROM  WHICH  THE  1-BITS  ARE  TO  BE  SELECTED. 

THIS  WORD  SHOULD  CONTAIN  AT  LEAST  AS  MANY  1-BITS 
AS  REOUESTED.  IF  MORE  THAN  THAT  ARE  REQUESTED, 

THIS  ROUTINE  WILL  RETURN  AN  EXACT  COPY  OF  TIC  INPUT. 

MBITS  -  NUMBER  OF  1-BITS  TO  BE  SELECTED  FROM  'IWORD'. 

OUTPUT  - 

LBITS  -  A  COPY  OF  THE  PORTION  OF  THE  IFFUT  WORD  CONTAINING 
THE  SPECIFIED  NUMBER  OF  LEFTMOST  1-BITS. 

COMMON  TABLES  USED  - 

I  COUNT ( I )  -  NUMBER  OF  1-BITS  IN  Tit  BINARY  REPRESENTATION  OF  THE 
INDEX  I.  1=0, 1,2,. ..,63 

MSKLFT(I)  -  MASK  FOR  WHICH  THE  LEFTMOST  I-BITS  ARE  1-BITS, 

AND  THE  REMAINDER  OF  THE  WORD  IS  ZERO.  1=1, 2, 3,. ..,36 

(WSK(I)  -  CONTAINS  A  1  IN  THE  1TH  BIT  (COUNTING  FROM  THE  LEFT) 
AND  O'S  EVERYWHERE  ELSE.  I=0,l,2,...,35 


Cm************************************************************m***m 


C— 


PARAMETER  MAXBIT=36.LFIELD=6 

COmON  /BITS/  NASK0,HASX(35),  fl£FT0,NSKLFTI36)> 

&  IZCOUT, IC0UNTI63) 

C— 

C—  *IF(NO  BITS  ARE  REQUESTED )TT£N 

IF(NBITS.GT.O)  GO  TO  1000 

C— 

C —  *RETURN  A  VECTOR  OF  ALL  ZEROES 

LBITS  =  0 

C— 

C—  *aSE  (SELECT  LEFTMOST  BITS) 

GO  TO  5000 
1000  CONTINUE 
C — 

C—  MSEARCHING  FROM  LEFT  TO  RIGHT,  FIND  TIC  6-BIT  FIELD 

C—  CONTAINING  THE  LAST  BIT  TO  BE  SELECTED 

C — 

C—  ^INITIALIZE  DO 

IBIT  =  0 
I  FOUND  =  0 

C— 

C—  *00  UNTIL  (APPROPRIATE  6-BIT  FIELD  IS  FOUND) 
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2000  CONTINUE 

C— 

C—  HPDATE  NUMBER  OF  1-BITS  FOUND  SO  FAR 

I  FOUND  =  I  FOUND  +  ICOUNT(FLD(IBIT,lFIELD,IMORD)) 

C— 

C —  UPDATE  FIELD  COUNTER 

IB1T  =  1B1T  +  LFIELD 

C— 

C—  *END  DO  (FIELD  LOOP) 

IF( (IFOUND.LT. NBITS)  .MD.  (IBIT.LT.MAXBIT) )  GO  TO  2000 

C — 

C —  ♦COMPUTE  NUMBER  OF  EXTRA  1-BITS  INCLUDED 

NEXTRA  =  I FOUND  -  NBITS 

C — 

C—  ♦PERFORM  ERROR  CHECK  TO  ENSURE  PROPER  NUMBER  OF  IS  FOUND 

IF(NEXTRA.LT.0)WRITE(6, 9001 ) NBITS, IFOUND, NEXTRA 

C— 

C—  ♦DO  WHIL£< THERE  ARE  EXTRA  l'S  TO  ELIMINATE) 

3000  CONTINUE 

IF(NEXTRA.LE.O)  GO  TO  4000 

C— 

C—  ♦DECREMENT  BIT  COUNTER 

IBIT  =  IBIT  -  1 

C— 

C—  ♦DECREMENT  EXTRA-BIT  COUNTER  IF  THIS  BIT  IS  A  1 

IF(AND( IMORB, MASK (IBIT) t.NE.O)  NEXTRA=NEXTRA-1 

C — 

C—  ♦END  DO  (EXTRA  l'S  LOOP) 

00  TO  3000 

4000  CONTINUE 

C— 

C—  ♦RETURN  PORTION  OF  INPUT  UP  TO,  BUT  NOT  INCLUDING 

C —  THIS  LAST  BIT 

LBITS  =  AND(  IMORD,  NSKLFT<  IBIT)) 

C— 

C—  ♦END  IF  (ZERO  BITS  REQUESTED  TEST) 

5000  CONTINUE 
C— 

RETURN 

9001  FORMAT ( "0$***t*$»  LBITS  ERROR  -  TOO  FEW  1-BITS  TO  MASK*,/, 

&  *  «****$$  NBITS,  IFOUND,  fCXTRA  =  *,315) 

END 
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C++  MAKEPD  -  CONVERTS  PARTS  LOW®  ARRAY  INTO  A  PuF. 

C++*  TO  ALIAS  METHOD  REQUIRES  A  LEGITIMATE  PROBABILITY 

C«rt  DISTRIBUTION  AS  M  TOUT.  THIS  ROUTINE  TAKES  TO  QPA 

C**+  AND  PROBI DEMAND 3  FIGURES  FOR  EACH  PART,  MD  FORMS  A  PDF 
C***  FROM  THEIR  PRODUCTS.  TORE  ARE  N  PART  TYPES. 


C— 

DIMENSION  IQPAIN).DEMAND(N) 

C— 

SUM  =  0.0 
DO  30  K=1.N 

DEMAND(K)  =  DEMANDS)  •  FLOAT!  IQPA(K) ) 
SUM  =  SUM  +  DEMAND(K) 

30  CONTINIE 

RECIP  =  1.0  /  SUM 
SUM  =  0.0 
DO  60  K=l.N-l 

DEHAND(K)  =  DEfttNDOO  *  RECIP 
SUM  =  SUM  +  DEMAND(K) 

60  CONTINUE 

DENAND(N)  =  1.0  -  SUM 

C — 

RETURN 

END 
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Cnnii«imn[imiiiwmmn»iniiiiii»»n»iniiimiHiiiimi 

INTEGER  FUNCTION  MNOM  (DUMMY) 

C++  (t«)  -  GENERATE  HULTINOniAL  SAMPLE  FOR  POT  DEMWffi  TYPE. 

Cm  THIS  FUNCTION  GENERATES  A  MULTINOMIAL  SAMPLE  INDICATING 
Cm  WHICH  PART  TYPE  (WS  BROKEN.  IT  USES  TWO  TABLES  CREATED 
Cm  PREVIOUSLY  BY  SUBROUTINE  'ALIAS'. 

Cm 

Cm  COMMON  INPUTS  - 

Cm  FPARTS  -  FLOATING-POINT  VALUE  OF  NUMBER  OF  PART  TYPES,  N. 

Cm  FRACT(I)  -  TABLE  OF  FRACTIONAL  CUTOFF  VALUES  USED  BY  TIC 

Cm  ALIAS  METHOD.  1=1,2 . N 

Cm  I  ALIAS!  I)  -  TABLE  OF  ALIASES  USED  BY  ALIAS  METHOD.  1=1,  ...,N 
Cm  OUTPUT  - 

Cm  MNOM  -  INDEX  INDICATING  TYPE  OF  POT  WHICH  HAS  BROKEN. 

Cm  MN0M=1,2,...,N  (NUMBER  OF  PART  TYPES) 

mmmitmmimmmmmmmmmmmmmmw 

C — 

PARAMETER  (MXPRT=304 
OXMON  /RSEED/  SEED 

COMMON  /ALIASC/  FRACT(MAXPRT),iALIAS(MAXPRT>,FPOTS 

C — 

C —  +MAKE  'U'  A  UNIFORM  (O.N)  RANDOM  REAL  NUMBER 
U  =  FPARTS  »  UNIFM1 (SEED) 

C —  * 

C—  *MAKE  'IU'  A  UNIFORM  RANDOM  INTEGER  (1,N) 

IU  =  IFIX(U)  +  1 

C — 

C —  «IF  (CCESSARY,  REPLACE  'IU'  BY  ITS  (LIAS 
IF  (U  .GT.  FRACT(IU))  IU  =  IALIAS(IU) 

C— 

C —  RESULT  IS  RETURNED  AS  'MNOM' 

MNOM  =  IU 

C — 

RETURN 

END 
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SUBROUTINE  HUPDATE(NUC,MA1NVC> 

CniiniiiimiiiHiii»innmi*iiiiii»»niiii[ii»iHiHnnmm 

C++  HUP DATE  -  UPDATE  MAINTENANCE  AIRCRAFT-STATE  BIT-VECTOR. 

Cm*  NUPDATE  UPDATES  TIC  OVERALL  MAINTENANCE  BIT-VECTOR.  THIS 
Cm  UPDATING  PROCESS  CONSISTS  OF  GOING  THROUGH  EACH  WORK-CENTER 
Cw  LIST  OF  AIRCRAFT  AND  MARKING  TIC  CORRESPONDING  BIT-POSITION  IN 
Cf«  THE  MAINTEWCE  BIT-VECTOR  FOR  ANY  AIRCRAFT  IN  SUCH  A  LIST;  THUS. 
C++*  mi  AIRCRAFT  IN  AT  LEAST  ONE  WORK-CENTER  LIST  WILL  BE  HARKED  AS 
C»*  BEING  IN  MAINTEMNCE  STATUS. 

Cm*  THE  MAINTENANCE  BIT-VECTOR  IS  UPDATED  AT  THE  BEGNWNG  OF 
Cm*  EACH  FLYING  CYCLE.  IT  IS  NOT  MAINTAINED  DURING  TIC  FLYING  CYCLE. 
Cm*  SINCE  IT  IS  0N.Y  NEEDED  AT  THE  START  OF  PREFLIGHT  TO  DETERM1IC 
CM*  THOSE  AIRCRAFT  WHICH  ME.  NOT  MISSION-CAPABLE  BECAUSE  TICY  ARE 
CM*  IN  AT  LEAST  ONE  WORK-CENTER.  IT  IS  MUCH  FASTER  TO  UPDATE 
CM*  ONCE  EACH  FLYING  CYCLE  RATHER  THAN  UPDATING  IT  EACH  TIME  A 
Cm*  WORK-CENTER  BREAK  OR  REPAIR  OCCURS. 


Cm* 

C*m  I  (PUTS  - 


Cm*  NWC  -  NUMBER  OF  WORK-CENTERS  BEING  SIMULATED 

Cm*  COMMON  INPUTS  - 

Cm*  HASK(I)  -  CONTAINS  A  1  IN  THE  ITH  BIT  (COUNTING  FROM  LEFT) 


C*M 

C«M 

Cm* 

Cm* 

Cm* 

Cm# 

Cm* 

Cm* 

Cm* 

Cm* 

Cm* 

C*m 

C*** 

C*** 

C**« 

Cm* 

Cm* 

Cm* 


AND  ZEROES  EVERYWHERE  aSE.  1=0,..., 35 
LENGTH  -  LENGTH  (IN  WORDS)  OF  AIRCRAFT  BIT-VECTORS 
INREPR(J)  -  NUMBER  OF  AIRCRAFT  IN  WORKCENTER-J. 

LISTRP(I.J)  -  LISTRPI  .  ,J)  IS  A  LIST  OF  AIRCRAFT  NUMBERS 

INDICATING  THOSE  AIRCRAFT  REQUIRING  MAINTENANCE  IN 
THE  J TM  WORK-CENTER  IJ=1,2,...,NWC).  THIS  LIST 
CONTAINS  EXACTLY  INREPR(J)  AIRCRtfT  NWffiERS.  TO  SAVE 
SPACE.  THESE  LISTS  HAVE  BEEN  PACKED  INTO  BIT-FIELDS 
INSTEAD  OF  WORDS.  EACH  NUMBER  IS  STORED  IN  A  FIT-FiaD 
■LFLD*  BITS  WIDE;  HENCE,  IF  'HAXBIT*  IS  THE  LENGTH 
OF  A  OmiTER  WORD  ON  THIS  SYSTEM,  THEN  THERE  ARE 
(MAXBIT/LFLD)  BIT-FIELDS  STORED  PER  WORD.  THE  AIRCRAFT 
NUMBERS  STORED  IN  T(CSE  BIT-FIELDS  INDICATE  A  UNIQUE 
BIT-POSITION  IN  TIE  VARIOUS  AIRCRAFT-STATUS  BIT- 
VECTORS.  THE  AIRCRAFT  ARE  NUMBERED, LEFT-TO-RIGHT , 
0,1,2,...,(MAXAC-1)  .  TO  GET  TIE  ITH  AIRCRAFT  NUMBER 
IN  A  WORK -CENTER  LIST,  THE  CORRESPONDING 
BIT-POSITION  WD  WORD-INDEX  MUST  BE  COMPUTED. 


CM*  OUTPUTS  - 
CM*  MAINVC 
Cm* 

Cm* 

Cm* 

Cm* 


-  MAINTENANCE  AIRCRAFT-STATUS  BIT-VECTOR.  EACH  BIT 
REPRESENTS  m  AIRCRAR.  A  1  INDICATES  THE 
CORRESPONDING  AIRCRAFT  JS  BEING  REPAIRED  IN  AT 
LEAST  ONE  WORK-CENTER,  AND  0  INDICATES  THE 
AIRCRAR  IS  NOT  CURRENTLY  IN  MAINTENANCE. 


PARAMETER  MAXWC-25 

PtfWCTER  NAXAC=108,HAXBIT=36,MAXVEC=2+(NAXAC-1)/MAXBIT 
PARAfCTER  LFLD=7,(PERWRD=«AXBIT/LFLD,flXINWC=I+(HAXAC-l)/NPERWRD 
COftlON  /WCmiNT/  LISTRP(MXINWC,MAXWL>,  INREPfK hAXWC) 


COMMON  /BITS /  MASKO,MASK(35),lt£RO,MSKLR(36), 
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i  izcout.  icami63) 

DIMENSiuN  MAIMVC(MftXVEC) 

C — 

C—  (INITIALIZE  MAINTENANCE  BIT-VECTOR  TO  NO  AIRCRAFT 
CALL  ZEROIMAINVC.MAXVEC) 

C— 

C —  *00  FORtEACH  WORK-CENTER) 

IF(NWC.EQ.O)  GO  TO  400 
DO  300  J=1,NWC 

C — 

C—  *00  WHILEfSTILL  AIRCRAFT  IN  THIS  WORK-CENTER  NAINTENANCE  LIST) 
NUN  =  INREPR(J) 

IF(NUH.EQ.O)  GO  TO  200 
DO  100  1=1. NUN 

C — 

C—  *GET  NEXT  AIRCRAFT  NUMBER  ON  LIST  FROM  tfPROPIATE  BIT-flELD 

I  AC  =  FU)  ( HOD  ( I-l>ffERWRLl)*LFLD,LFLD, 

L  LISTRPd-HI-D/NPERHRD.J)) 

C— 

C —  (COMPUTE  WORD  AND  BIT  POSITIONS  IlfflICATED  BY  THIS  AC  I 

IWORD  =  2  +  IAC/HAXBIT 
IBIT  =  HOD( IAC.HAXBIT) 

C— 

C —  *HARK  CORRESPONDING  POSITION  IN  AIRCRAFT  HUNT  BIT-VECTOR 

NAINVCf IWORD)  =  OR(HAINVC( IWORD),  HASKIIBIT)) 

C— 

C —  *EM?  DO  (WORK-CENTER  LIST  LOOP) 

100  CONTINUE 

200  CONTINUE 

C — 

C—  *£ND  DO  (WORK-CENTER  LOOP) 

300  CONTINUE 
400  CONTINUE 
C — 

C —  *am/TE  TOTAL  AIRCRAFT  IN  NAINTENANCE  BY  COUNTING  1-BITS  IN 
C—  AIRCRAFT  MAINTENANCE  VECTOR 

HAINVCU)  =  NIVECT(HAINVC) 

C — 

RETURN 

END 
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CHmHIHHmHHHIIIHHHHtHHHIIWHHfmHH«HHimHH 

INTEGER  FUNCTION  NIBITS(IWORD) 

CmiiniiWHi»niH«mmn>HHnmnnm«mniHiimnm 

C++  N1BITS  -  COUNT  NltfflER  OF  1-BITS  IN  A  COMPUTER  WORD. 

Cm  N1BITS  IS  A  FORTRAN  SUBROUTINE  WHICH  WILL  RETURN  TTC 
Cm  NUMBER  OF  1-BITS  IN  A  GIVEN  WORD.  THIS  ROUTINE  IS  SPECIFIC 
Cm  TO  A  COfFUTER  WITH  36-BIT  WORDS,  SINCE  IT  WORKS  BY  EXTRACTING 
Cm  6-BIT  FIELDS  FROM  THE  WORD.  IT  USES  EACH  6-BIT  FIELD  EXTRACTED 
Cm  AS  AN  INDEX  INTO  A  TABLE,  M  TTC  ENTRIES  IN  THE  TABLE  CONTAIN  THE 
Cm  CORRESPONDING  NUMBER  OF  1-BITS  FOR  THAT  INDEX. 

Cm 

Cm  INPUT  - 

Cm  1U0RD  -  WORD  FOR  WHICH  THE  1-BITS  ARE  TO  BE  COUNTED 
Cm  OUTPUT  — 

Cm  MBITS  -  NltfflER  OF  1-BITS  IN  THE  GIVEN  INPUT  WORD.  TCNCE, 

Cm  MBITS  RETURNS  AN  INTEGER  BETWEEN  0  AND  36. 

Cm  TABLE  USED  - 

Cm  I  COUNT!  I)  -  NUMBER  OF  1-BITS  IN  THE  BINARY  REPRESENTATION  OF  THE 

Cm  INDEX  I.  1=1,. ..,63 

CHHIHHHHIIIIIHIHIHmilHIIHHHHmiHH»limiHHHHIH 

C — 

COMMON  /BITS/  MASK0,(WSK(35),  NLEFTO,  HSKIFK36), 

*.  IZCOUT,  IC0UNT<63) 

C— 

N1BITS  =  ICOUNT!  FLD(  0.6.IWORD)  ) 

L  +  I  COUNT!  FLD(  6.6.IW0RD)  ) 

L  +  ICOUNT!  FLD(12,6, IWORD)  > 

L  +  ICOUNT!  FLD(18,6, IWORD)  ) 

fc  +  ICOUNT!  FLD(24,6, IWORD)  > 

l<  +  ICOUNT!  FLD(30,6, IWORD)  ) 

C— 

RETURN 

END 
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C++  N1VECT  -  COUNT  NUMBER  OF  1-BITS  IN  A  BIT-VECTOR. 

Cm*  N1VECT  IS  A  FORTRAN  SUBROUTINE  WHICH  WILL  RETURN  THE 
C*«*  NUMBER  OF  1-BITS  IN  THE  WORDS  COMPRISING  A  GIVEN  IMVT  ARRAV. 

Cm*  NOT  INCLUDING  TIC  FIRST  WORD.  THIS  ROUTINE  IS  USED  TO  COUNT 
Cm*  THE  NUMBER  OF  1-BITS  IN  TIC  VARIOUS  AIRCRAFT-STATUS 
Cm*  BIT-VECTORS. 

C*m  THIS  ROUTINE  IS  SPECIFIC  TO  A  36-BIT-WORD  COWUTER,  SINCE  IT 
Cm*  WORKS  BY  EXTRACTING  6-BIT  FIELDS  FROM  THE  ARRAY  WORDS.  IT  USES 
Cm*  EACH  6-BIT  FIELD  EXTRACTED  AS  AN  INDEX  INTO  A  TABLE.  AM)  Tt£ 

Cm*  ENTRIES  IN  THE  TABLE  CONTAIN  THE  CORRESPOM1ING  NUMBER  OF  1-BITS 
C*m  FOR  THAT  INDEX. 

Cm* 

Cm*  INPUT  - 

C*m  I ARRAY  -  ARRAY  FOR  WHICH  THE  1-BITS  ARE  TO  BE  COUNTED. 

Cm*  COrtlON  TABLE  USED  - 

C+m  ICOUNT(I)  -  NUMBER  OF  1-BITS  IN  THE  BINARY  REPRESENTATION  OF  THE 
C+m  INDEX  I.  1=1.. ...63 

Cm*  OUTPUT  - 

C+M  N1VECT  -  NUMBER  OF  1-BITS  IN  THE  GIVEN  IICUT  ARRAY. 

Cm*  EXCLUDING  TIC  FIRST  WORD. 


C— 

PARAMETER  f»XAC=  108,  HAXBI T=36,  fWXVEC=2+ ( HAXAC-1 ) /MAXBIT 
COMMON  /BITS/  MSK0,MASK(35),  ML£FT0»HSKLFT(36)> 

4  I ZCCJUT ,  ICOUNT  <  63  > 

CWflON  /ACSTATE/  LENGTH,  NACVC !  MAXVEC ) ,  IFLYVC(MAXVEC), 

4  MAINVC(MAXVEC)»NORSVC(MAX\CC),  LOSTVCI MAXVEC) 

DIMENSION  I ARRAY ( 1 ) 

C— 

N1VECT  =  0 
DO  1000  1=2, LENGTH 
IWORD  =  I ARRAY ( I ) 

N1VECT  =  N1VECT  +  ICOUNTI  FLD(  0,6, IWORD)  ) 

4  ♦  ICOUNTI  FUJI  6,6,  IWORD)  ) 

&  +  ICOUNTI  FUJI  12,6,  IWORD)  > 

4  +  ICOUNTI  FLDI 18,6, IWORD)  ) 

4  ♦  ICOUNTI  FLDI24.6, IWORD)  ) 

4  +  ICOUNTI  FLDI 30, 6, IWORD)  ) 

1000 
C — 

RETURN 
END 


CONTINUE 


C< 


INTEGER  FUNCTION  NB I NON  (PBINOH. NTRYS) 
C«»HiHiin»iHiminuiinm»n»Himiiiniiimnimn» 

C++  (CINCH  -  GENERATE  RANDOM  SAMPLE  FROM  BINOMIAL  DISTRIBUTION. 

Cm  (ffilNOM  GOERATES  A  RANDOM  SATFLE  FROM  A  BINOMIAL  DISTRIBUTION 
Cm  WITH  THE  GIVEN  INPUT  CHARACTERISTICS.  THIS  ROUTINE  USES  A 
C**+  COMBINATION  OF  TWO  (CTHODS  TO  GEfERATE  THIS  SAMPLE.  FOR 
Cm  BINOMIALS  WITH  RELATIVELY  SMM±  NUMBERS  OF  TRIALS.  TIE 
Cm  STRAIGHTFORWARD  BERNOULLI  TRIALS  METHOD  IS  USED.  FOR  LARGER 
Cm  VALUES.  THE  INVERSE  TRANSFORM  (CTHOD  IS  USED. 

Cm  NOTE  THAT  THE  NUtffiER  OF  FAILURES  IN  A  BINOMIAL  SATFLE  IS 
Cm  THE  COMPLEIENT  OF  TIE  NUMBER  OF  SUCCESSES  IN  THAT  DRAW.  HENCE 
Cm  THIS  ROUTINE  WILL  SAMPLE  FROM  THE  COWLEJENTARY  BINOMIAL 
Cm  DISTRIBUTION  OF  FAILURES  HtfN  TIE  PROBABILITY  OF  SUCCESS  IS 
Cm  GREATER  TWN  .5  . 

Cm 

Cm  INPUTS  - 

Cm  PBINOH  -  PROBABILITY  CHARACTERISTIC  OF  THE  BINOMIAL. 

Cm  PBINOH  ALSO  EQUALS  THE  PROBABILITY  THAT  THE 

Cm  BERNOULLI  VARIABLE  UNDERLYING  THIS  BINOMIAL  EQUALS  1 

Cm  NTRYS  -  NUMBER  OF  BERNOULLI  TRIALS  CHARACTERIZING  THIS 
Cm  BINOMIAL. 


C— 

COMMON  /RSEED/  SEED 

C— 

C—  ^INITIALIZE  SAMPLE  TO  NO  SUCCESSES 
NBINOH  =  0 

C — 

C —  *IF ( THIS  IS  NOT  A  SPECIAL  DISTRIBUTION  TO  BE  HANDLED  SEPERATELY) 
IF( (PBINOM.LE.O.O)  .OR.  ( PBINOH. GE.  1.0) 

fc  .OR.  (NTRYS.LE.*))  GO  TO  3000 

C— 

C —  *DRAU  RANDOM  SAMPLE  FROM  UNIFORM  (0, I)  DISTRIBUTION 

RDRAM  *  UNIFMUSEED) 

C — 

C —  tDETERMINE  WHETHER  TO  SAfFLE  SUCCESSES  OR  FAILURES 

PFAIL  =  AMAXI (PBINOH. l.O-PBINOM) 

PSUCC  =  1.0  -  PFAIL 

C — 

C—  COMPUTE  QUICK  (FPROXIMATION  TO  PROBIO  SUCCESSES) 

PROB  =  1.0  -  FLOAT!  NTRYS  )*PSUCC 
IF(RDRAW.Lf.PROB)  GO  TO  2000 

C — 

C—  COMPUTE  EXACT  PROBABILITY  OF  NO  SUCCESSES 
PROB  =  PFAILMNTRYS 

C — 

C —  »IF(RANDOM  DRAW  DOES  NOT  FALL  WITHIN  THIS  PORTION  OF  THE  CDF ) 

IF(RORAW.LE.PROB)  GO  TO  2000 

C — 

C—  ♦INITIALIZE  LOOP  TO  FIND  APPROPRIATE  PLACE  IN  THE  CDF 

RATIO  =■  PSUCC/PFAIL 
(FLUS1  *  NTRYS  ♦  1 


CDF  =  PROB 

C — 

tDO  UNTIL(APPROPRIATE  CDF  INDEX  IS  FOUtffi) 

1000  CONTINUE 
C— 

C —  *UPDATE  SAMPLE  COUNTER 

NBINOM  =  NBINOM  ♦  1 

C— 

C —  COMPUTE  NEXT  ENTRY  IN  CUMULATIVE  DISTRIBUTION  FUNCTION 

PROB  =  <FLQAT(NPLUSl-f€INOH)/FLQAT(NBINOM>)»RAn(}*PROB 
CDF  =  CDF  ♦  PROB 

C— 

C—  ♦£»  DO  (CDF  LOOP) 

IFURDRAU.GT.CDF)  .AND.  (NBINOM.LT.NTRYS))  GO  TO  1000 

C — 

C —  «£ND  IF  (0  SUCCESSES  TEST) 

2000  CONTINUE 

C — 

C —  COMPLEMENT  RESULT  IF  FAILURES  WERE  SAMPLED 

IFIPBINOM.GT.  .5)  N8IN0M  =  NTRYS  -  WINOM 

C— 

C—  «€LSE  (SPECIAL  CASES) 

GO  TO  7000 
3000  CONTINUE 

C— 

C—  *IFCmiS  IS  A  DEGEMOTIVE  DISTRIBUTIONCTHEN 

IF(  (PBINOM.GT.O).AW.  (PBINOH.LT.  l.Ol.AND.  (NIRYS.GT.O) ) 
fc  GO  TO  4000 

C — 

C—  *SA»>L£  FROM  DISTRIBUTION  (IF  PBINWH),  OR  NTRYS=0 

C—  THEN  WE  (WE  DOE) 

IFIPBINON.GE. 1.0)  NBINOM  =  NTRYS 

C— 

C—  »ELSE  (USE  BERNOULLI  TRIAL  METHOD) 

GO  TO  6000 
4000  CONTINUE 

C — 

C —  *PERFORM  APPROPRIATE  NUMBER  OF  BERNOULLI  TRIALS 

DO  5000  1=1, NTRYS 

IF(UNIFNKSEED).LE.PBINOM)  NBINOM  =  NBINOM  +  1 
5000  CONTINUE 
C— 

C—  *£ND  IF  (DEGENERATIVE  DISTRIBUTION  TEST) 

6000  CONTINUE 

C — 

C —  IF  (SPECIM.  CASES  TEST) 

7000  CONTINUE 

C— 

RETURN 

EM) 
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INTEGER  FUNCTION  NDMDSl  NBRKAC) 


CHiimiiiiiiiiiHiwwimmuimiHiHUHiHiwmmiumi 
C++  MMDS  -  GENERATE  SAMPLE  OF  TOTAL  SORTIE  PART  DEMANDS. 

Cm  NDMNDS  IS  A  FORTRAN  FUNCTION  WHICH  GENERATES  A 
Cm  SAtFLE  NUMBER  OF  PARTS  DEMANDS  ON  A  SORTIE,  GIVEN  THE  TOTAL 
Cm  NUMBER  OF  AIRCRAFT  WHICH  BROKE  ON  THAT  SORTIE.  THE 
Cm  PROBABILITY  DISTRIBUTION  OF  TOTA.  PARTS  DEMANDS  IS  APPROXIMATED 
Cm  USING  EITHER  A  NORMAL  DISTRIBUTION  (IF  MEAN  IS  LARGE  ENOUGH  TO 
Cm  APPLY  TO  CENTRAL  LIMIT  THEOREM)  OR  A  POISSON  DISTRIBUTION. 

C<HHt 

Cm  INPUTS  -- 

Cm  NBRKAC  -  NUMBER  OF  AIRCRAFT  WHICH  BROKE  ON  TO  SORTIE 
Cm  COMMON  INPUTS  - 
Cm  ACf€AN  -  EXPECTED  VALUE  OF  TO  RWfflOH  VARIABLE  REPRESENTING 
Cm  THE  NUMBER  OF  PARTS  DEMANDS  PER  AIRCRAFT,  GIVEN 

Cm  THAT  TO  AIRCRAFT  HAS  BROKEN  IRON  RETURNING  FROM 

Cm  A  SORTIE. 

Cm  ACVAR  -  VARIANCE  OF  TOTAL  PARTS  DEMAND  PER  BROKEN  AIRCRAFT 

Cm  NPERAC  -  TOTAL  NUMBER  OF  PARTS  PER  AIRCRAFT.  THIS  IS  USED  TO 

Cm  ENSURE  THAT  A  LEGITIMATE  SAWLE  IS  GO€RATED. 


c — 

PARAMETER  CUT0FF=O.0 
COMMON  /RSEED/  SEED 
COftWN  /DEMAND/  A0EAN,  ACVAR,  f#>ERAC 

C — 

C —  +INITIALIZE  SAMPLE  TO  NO  PARTS  DEMANDS 

NDMNDS  =  0 

C — 

C —  *IF(TOR£  WERE  ANY  BROKEN  AIRCRAFT) TON 

IF (NBRKAC. Eff.O)  GO  TO  300 

C— 

C—  COMPUTE  MEAN  OF  DISTRIBUTION  OF  TOTft.  DEMANDS 

C —  CORRESPONDING  TO  NUMBER  OF  BROKEN  AIRCRAFT 

FLTAC  =  FLOAT(NBRKAC) 

BMEAN  =  FLTAC  ♦  ACMEAN 

C — 

C—  *IF(EXPECTED  TOT«.  DEMANDS  IS  SMALL) TON 

IF(BI€AN.GT. CUTOFF)  GO  TO  100 

C— 

C—  HJSE  POISSON  APPROXIMATION 

NDMNDS  *  I  POISSON)  BftAN,  SEED) 

C— 

C —  *€LSE 

GOTO  200 
100  CONTINUE 

C— 

C—  +USE  NORMAL  APPROXIMATION 

BSTDEV  =  SORT ( FIT AC+ACVAR ) 

NDMNDS  =  HAX0<0, INT ( XN0RH<  BMEAN, BSTDEV, SEED>+.5) ) 
C— 

C—  *END  IF  (APPROXIMATION  TYPE  TEST) 


B-43 


o  o 


200 


CONTINUE 


♦ENSURE  THAT  A  FEASIBLE  ANSWER  HAS  BEEN  GOtRATED 
MimS  =  MN0()«tffiS.NPERAC*N8RKAC) 

C — 

C—  *END  IF  (ZERO  BROKEN  AC  TEST) 

300  CONTINUE 
C — 

RETURN 

END 
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CHHHHiH»HHHHIHHIIHmiHHHHHHIHHIHHHHIH»WI 

INTEGER  FUNCTION  NORSAC (NPARTS, Ift'AiNBACKO) 


CHHmmmUOHUlimHUmiHHHHHIIIHmHHHHWHH 

C++  NORSAC  -  CALCULATE  INITIAL  NUMBER  OF  NORS  AIRCRAFT. 

C++*  NORSAC  IS  A  FORTRAN  FUNCTION  WICH  CALCULATES  TO  CURRENT 
Cm#  NUMBER  OF  NORS  AIRCRAFT  -  ASSUMING  PERFECT  CANNIBALIZATION. 

C++* 

Cm#  INPUT  - 


Cm# 

Cm# 

C++* 

Cm# 

Cm# 

Cm* 

Cm# 


NPARTS 

IQPA(K) 

NBACXO(K) 


-  TOTAL  NltfflER  OF  PART  TYPES. 

-  NUMBER  OF  TYPE-*  PARTS  INSTALLED  ON  EACH  AIRCRAFT, 

-  NUMBER  OF  BACKORDERS  FOR  PttTS  OF  TYPE-K. IF 
NBACKO(K)  IS  POSITIVE,  THEN  UNFULFILLED  REQUESTS 
FOR  PARTS  OF  THIS  TYPE  HAVE  BEEN  MADE.  IF  IT  IS 
fCGATIVE,  THEN  JfiACKOIK)  INDICATES  THE  NUMBER  OF 
OF  PARTS  ON-THE-SHELF. 


Cm#  OUTPUT  - 

C++*  NORSAC  -  CURRENT  NUMBER  OF  NORS  AIRCRAFT  BASED  ON  THE  GIVEN 
Cm#  BACKORDER  MD  (PA  INFORMATION  AM)  ASSUMING  PERFECT 
Cm#  CANNABILIZATION. 


CMM#*MHIM»IMMmMM»M#M«MMMHMH«M*M 

c— 

DIMENSION  NBACKO(NPARTS),  IQPAdPARTS) 


C — 

C —  *INITIALIZE  NUMBER  OF  NORS  AIRCRAFT  TO  NONE 

NORSAC  =  0 

C— 

C—  *DO  FClRfEACH  PART  TYPE) 

DO  2000  K- l, NPARTS 
C — 

C —  * IF (THESE  PARTS  CAUSE  THE  MAX  NUMBER  OF  NORS  THUS  FAR) 

IF(NORSAC*IQPA(K)  .GE.  NBACKO(K))  GO  TO  1000 

C— 

C —  UPDATE  NUMBER  OF  NORS  AIRCRAFT 

N0RSAC=1NT(FL0AT(NBACK0(K)) /FLOAT* IQPA(K) )  +  .999) 

C — 

C—  *END  IF  (NEW  NORS  MAXIMUM  TEST) 

1000  CONTINUE 

C — 

C —  *END  DO  (PARTS  LOUP) 

2000  CONTINUE 

C— 

RETURN 

END 
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INTEGER  FUNCTION  NORSBK ( NfiRKAC > NOROLO > 


CHiHHniiiniimnninnmmmnmH«nmiitiiniimiuin 

C++  NORSBK  -  DETERMINES  NORS  AIRCRAFT  FROM  A  SORTIE. 

Cm  NORSBK  IS  A  FORTRAN  FUNCTION  WHICH  CALCULATES  THE  NUMBER 
Cm  OF  NORS  AIRCRAFT  RESULTING  FROM  A  SORTIE  WITH  A  SPECIFIED  NUMBER 
Cm  OF  BROKEN  AIRCRAFT  -  ASSUMING  IMMEDIATE  AND  MAXIMUM 
Cm  CAMABILIZATIGN  OF  PfWTS.  NORSBK  DETERMI!€S  THE  TOTAL 
Cm  NUWER  AND  DISTRIBUTION  OF  THE  PARTS  OEfWfflS  RESULTING  FROM  THIS 
Cm  SORTIE.  IT  UPDATES  FOR  EACH  P«T  TYPE  DEMANDED.  THE 
Cm  NUMBER  OF  PARTS  ON-THE-SHELF,  BACKORDERED,  AW)  IN  RESUPPLY. 


Cm 


Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 


INPUTS  - 

NBRKAC  -  NUMBER  OF  AIRCRAFT  WHICH  BROKE  DURING  THE  SORTIE 

NOROLD  -  NUMBER  OF  NORS  AIRCRAFT  BEFORE  THIS  LATEST  SORTIE. 

COMMON  HEUTS  - 


IQPAOO  -  NUMBER  OF  TYPE-K  PMTS  INSTALLED  ON  EACH  AIRCRAFT. 

INFPART  -  LOGICAL  FLAG  IMIICATING  WHETHER  THE  INFINITE  PARTS 
ASSUMPTION  HOLDS.  IF  INFART  IS  TRUE  THEN  THERE 
IS  NEVER  ANY  SHORTAGE  OF  PARTS;  HENCE.  NO  NORS  AC. 

COMMON  INPUTS/OUTPUTS  - 

NBACKO(K)  -  NUMBER  OF  BACKORDERS  FOR  PARTS  OF  TYPE-K. IF 

NBACKO(K)  IS  POSITIVE.  THEN  UNFULFILLED  REQUESTS 
FOR  PARTS  OF  THIS  TYPE  MVE  BEEN  HADE.  IF  IT  IS 
fFGATIVE,  TIEN  NBACKO(K)  INDICATES  THE  NUMBER 


Cm 

Cm  OUTPUT  - 
Cm  NORSBK 
Cm 


OF  PARTS  ON-THE-SHELF. 

-  NUMBER  OF  NORS  AIRCRAFT  AT  THE  END  CF  THIS  SORTIE 
ASSUMING  MAXIMUM  AND  IMMEDIATE  CWKBILIZATION. 


CWUHtlHHUHHIH; 

C — 


PARAMETER  MAXPRT=304,  MAXCYC=I0 

COMMON  /PARTS/  fPARTS , IQPA ( HAXPRT ) , N8ACKQ ( MAXPRT ) , 

&  BflPRATE(MAXPRT) ,DRPRATE(MAXPRT>,INITSJ(MAXPRT),RESUPP(HAXPRT). 

L  BNRTS<  HAXPRT ) . NBASE ( MAXPRT ) , NDEPOT ( MAXPRT ) 

COMMON  /INPUT/  INITUE,  NAC.  PATTRIT,  IRES.  RNMCM,  DEPART, 
i  HAXFLY(MAXCYC),  INFMAN,  ISCALE,  IAUGMNT 

LOGICAL  INFART 

C— 

C —  +INITIALIZE  f€W  NORS  TO  OLD  NUMBER  OF  NORS  AIRCRAFT 

NORSBK  =  NOROLD 

C— 

C —  »IF(THERE  ARE  ANY  BROKEN  AIRCRAFT  AND 


C —  INFINITE  PARTS  NOT  ASSUMED) THEN 

IF(NBRKAC.EQ.O)  GO  TO  5000 
IF(  DEPART)  GO  TO  5000 

C— 

C—  ♦  DETERMINE  TOTW.  NUMBER  OF  PARTS  DOVWDS  FROM  T)€SE  BROKEN  AC 

WENS  =  NDMNDSI  NBRKAC) 

C — 

C—  *IF(ANY  PARTS  WERE  DEMANDED  (THEN 

IF(NDEMS.EQ.O)  GO  TO  4000 

C — 


oo  o  o  o  o  o  o  o  o  ooo*  <->  o  <r>  o  o  o  ooo 


C—  «D0  FOR  (EACH  PMT  DEMAND) 

00  3000  M.NDEMS 

—  *DETERMINE  P/WT-TYPE  FOR  THIS  OEHAK)  BY  SAMPLING 

—  FROM  A  MULTN0MIAL  DISTRIBUTION 
KTYPE  =  lf«M() 

—  *UPDATE  BACKORDERS  FOR  THIS  PART  TYPE 

NBACKO(KTYPE)  =  NBACKO(KTYPE)  +  l 

—  *IF(AN  UNFULFILLED  CGWffl  HAS  OCCURRED) THEN 

IF<NBACK0(KTYPE).L£.0>  GO  TO  2000 

—  *IF(THIS  DEMAND  CAUSES  A  NEW  NORS  ACJTHEN 

IF(NORSBK*IQPA(KTYPE) . GE.NBACKO(KTYPE) ) 

GO  TO  1000 

—  * INCREMENT  NUMBER  OF  NORS  AIRCRAFT 
NORSBK  =  NORSBK  +  l 

1F((NORSBK-NOROLD).GE.NBRKAC)  GOTO  5000 

—  *END  IF  (NEW  NORS  AIRCRAFT  TEST) 

1000  CONTINUE 

—  *END  IF  (UNFULFILLED  DEMAND  TEST) 

2000  CONTINUE 

*£W)  DO  (DEMAND  LOOP) 

3000  CONTINUE 

—  *END  IF  (ZERO  DEMANDS  TEST) 

4000  CONTINUE 

—  #END  IF  (NO  BROKEN  AC  OR  INFINITE  PARTS  TEST) 

5000  CONTINUE 

C — 

RETURN 

END 
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INTEGER  FUNCTION  H$PS(TIMET,NR£PJ,NCRWSJ,SRATEJ> 


C++  fUEPS  -  RANDOM  SAMPLE  OF  AIRCRAR  REPAIRS  IN  A  WORK  CENTER. 
C++*  NREPS  IS  A  FORTRAN  FUNCTION  WHICH  FOURNS  A  SAMPLE  NUMBER 
Cm*  OF  AIRCRAFT  REPAIRED  IN  A  WORKCENTER,  BASED  CM  THE  LENGTH  OF  THE 
C+*+  REPAIR  PERIOD,  NUMBER  OF  SERVERS,  REPAIR  RATE  FOR  EACH  SERVER, 
Cm*  AND  THE  NUMBER  OF  AIRCRAFT  IN  THE  WORKCENTER  AT  THE  START  OF 
C»*  THE  REPAIR  PERIOD.  IT  IS  ASSUMED  THAT  NO  F£W  AIRCRAR  ARRIVE 
Cm*  DURING  THE  REPAIR  PERIOD,  REPAIR  TIMES  ARE  EXPONENTIALLY 
Cm*  DISTRIBUTED,  WITH  THE  SAME  DISTRIBUTION  APPYING  TO  EACH  SERVER 
Cm*  INDEPENDENTLY. 

C*M 

cm*  inputs  - 

C*M  TIMET  -  LENGTH  (IN  H<XRS)  OF  Tl€  REPAIR  PERIOD. 

Cm*  W£PJ  -  NUMBER  OF  AIRCRAR  IN  T>£  WORKCENTER  AT  THE 

Cm*  START  OF  THE  REPAIR  PERIOD. 

CM*  NCRWSJ  -  NUMBER  OF  REPAIR  CREWS  (SERVERS)  FOR  THIS 
C*m  WORKCENTER. 

Cm*  SRATEJ  -  REPAIR  RATE  (AIRCRAR/ HOUR)  FOR  EACH  CREW  IN 
Cm*  THIS  WORKCENTER. 

Cm*  COrtlON  INPUT/OUTPUT  - 

C*m  SEED  -  SEED  FOR  RANDOM  NUMBER  GENERATOR. 

Cm*  OUTPUT  - 

CM*  NREPS  -  NUMBER  OF  AIRCRAR  REPAIRED  IN  THIS  WORKCENTER 
Cm*  DURING  THE  REPAIR  PERIOD.  NR£PS=0, 1,2,...  .«£PJ 


C— 

COMMON  /RSEED/SEED 

C — 

C —  *IF(NUMBER  OF  AC  IN  REPAIR  IS  LESS  THAN  THE  NUMBER  OF  CREWS) 

IF(NR£PJ.GT. NCRWSJ)  GO  TO  1000 

C — 

C —  *0ETERMI(€  NUMBER  OF  AIRCRAR  REPAIRED  BY  SAMPLING 

C —  FROM  TF£  APPROPRIATE  BINOMIAL  DISTRIBUTION 

NREPS  =  hFEPJ  -  NBINGM(  EXP(-SRATEJ*TIMET)  ,  NREPJ  ) 

C — 

C —  *aSE  (MORE  AIRCRAR  THAN  CREWS) 

GO  TO  4000 
1000  CONTINUE 

C — 

C—  ^INITIALIZE  VARIABLES 

NREPS  =  0 
CUMP  =  1.0 

MAXREP  =  «£PJ  -  NCRWSJ  ♦  1 
CWRATE  =  FLOAT ( NCRWSJ ) *SRATE J 
EXPTYM  =  EXP(-CWRATE»TIMR) 

C — 

C —  *00  UNTIL! A  SERVER  BECOMES  IDLE  OR  TTC  NEXT  AIRCRAR 

C —  DEPARTURE  TIFt  EXCEEDS  LENGTH  UF  REPAIR  PERIOD) 

2000  CONTINUE 

C—  *GEF€RATE  W  ACCUMULATE  NEXT  AIRCRW^I  DEPARTURE  FROM 
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C—  THIS  WORKCENTER 

CUMP  =  COT  *  UNIFHKSEED) 

C— 

C—  *EXIT  LOOP,  IF  REPAIR  TI)CS  EXCEED  TIPE  INTERVAL  LENGTH 

IFICUHP  .LT.  EXPTYH)  GO  TO  3000 

C — 

C—  ♦INCREPENT  NUMBER  OF  AIRCRAFT  REPAIRED 

HEPS  =  HEPS  +  1 

C— 

C —  #£ND  DO  (REPAIRED  AIRCRAFT  LOOP) 

IF(MEPS.LT.HAXREP)  GO  TO  2000 

C — 

C—  «A  SERVER  «S  JUST  BECOPE  IDLE,  PERFORM  A  BINOMIAL  ORAN 

C —  TO  DETERMINE  HOU  MANY  MORE  AC  ARE  REPAIRED 

C— 

C—  COMPUTE  TIPE  LEFT  IN  THE  INTERVAL 

C —  (LENGTH  OF  REPAIR  PERIOD)  -  (THE  OF  LAST  REPAIR) 

TLEFT  =  TlfET  +  ALOG(CUMP) /CURATE 

C— 

C —  COMPUTE  PROBABILITY  AN  AIRCRAFT  IS  NOT  REPAIRED 

C—  IN  THE  REMAINDER  OF  THE  INTERVAL 

PNOREP  =  EXP(-SRATEJ»TL£FT) 

C — 

C —  CENERATE  A  BINOMIAL  ORAN  TO  DETERMIIC  NUMBER  OF 

C —  REMAINING  AIRCRAFT  WHICH  ARE  NOT  REPAIRED 

NOTREP=PffiINOM(PNOREP,NCRWSJ-i ) 

C— 

C—  COMPUTE  TOTAL  AIRCRAFT  REPAIRED  DURING  PERIOD 

NREPS=NREPJ-NOTREP 

C— 

C— 

C —  CXIT  FROM  DO  LOOP 

3000  CONTINUE 

C— 

C —  *END  IF  (MORE  AIRCRAFT  THAN  CREWS  TEST) 

4000  CONTINUE 
C— 

RETURN 

END 
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CHiimniiniiiui»tmumiiiimnuimnmniinninnniHin 

SUBROUTINE  PRINTO 

CHHIlHIHHIHIIIIIIIIIIIIHHWmummHHHHHUHIIIIHIH 

C++  PRINTO  -  PRINT-OUT  RESULTS  OF  THE  SIMULATION  RUN. 

C++*  THIS  ROUTINE  PRINTS  THE  RESULTS  OF  T)€  SGH  SIMULATION. 

Cm  THESE  RESULTS  CONSIST  OF  THE  AVERAGE  NUMBERS  OF  AIRCRAFT  IN 
Cm  THE  VARIOUS  POSSIBLE  AIRCRAFT  STATES  AT  TIC  START  OF  EACH 
Cm  SORTIE  PERIOD  FOR  EACH  FLYING  DAY  OF  THE  SCENARIO.  THE  AVERAGE 
Cm  SORTIES  PER  AIRCRAFT  PER  DAY  IS  COfCUTED  AND  ALSO  PRINTED. 

Cm  TICSE  SORTIES  PER  AIRCRAFT  PER  DAY  FIGURES  AMD  THE  TOTH.  SORTIE 
Cm  PRODUCTION  PER  DAY  ARE  PRINTED  TO  A  SCRATCH  FILE  (FILE  07)  TO  BE 
Cm  USED  AS  INPUT  FOR  SORTIE  PLOTS. 

CHHIIIHIHIllimWHIlHIIIIIHHIIHHmHHHHttllHHHHWI 

C— 

PARAMETER  NAXAC=108,MAXWC=25,MAXBIT=36,IWraT=304, 

&  MAX VEC=2+ ( MAX AC- 1 ) /MAXBI T 

PARAMETER  MAXDAY=30,MAXCYC=10,HAXSTAT=5 
COMION  /STATS/  EXPECT (MAXSTAT . MAXCYCi MAXDAY ) . 

St  NRESRVi  I ZDAY  >  ITOTRESt  MAXDAY )  >  LOSSTOT 

COMMON  /TIME/  PRERITE.SORTLGTH.UAITCYC, 

St  TYMNITE.NSIM,  IS1M<NUMDAY<  IDAYtNCYCLESt  ICYCLE 

COMMON  /INPUT/  INITUfc.NAC.PATTRIT,  IRES.RNMCM.DFPART. 

St  MAXFLY(MAXCYC) .  INFMAN,  I  SCALE.  IAUGMNT 

C0M10N  /UCINPUT/  NUC.  NCREUS(HAXUC).  SRATE(HAXUC) 

C — 

C—  +PRINT  EXPECTED  NUMBER  OF  AVAILABLE  AIRCRAFT  FOR  EACH  SORTIE  PERIOD 
WRITE  (6.9005) 

TOTFLY  =  0.0 

WRITE  (7)  FLOAT ( ISCALE) . NUMDAY 
FSIM  =  FLOAT(NSIM) 

DO  6000  J=l, NUMDAY 

SORTYDAY  =  EXPECT(l.l.J) 

OFFSCENE  =  EXPECT(4, 1, J)+£XPECT(5, 1, J) 

WRITE<6,9008)  J. 1 , < EXPECT (M, 1 , J) /FSIM.  1 , 5) 

DO  5000  I =2 ,  MCYCLES- i 

WRITE  (6.9006)  MEXPECT(M,I„J)/FSIH,H=M) 

SORTYDAY  =  SORTYDAY  +  EXPECT(l.I.J) 

OFFSCENE  =  OFFSCEfC  ♦  EXPECT ( 4, I , J )  +  EXPECT(5.I,J) 

5000  CONTINUE 

I  =  NCYCIES 

SORTYDAY  =  (SORTYDAY  ♦  EXPECTS  1 . 1 ,J> » /FSIM 

OFFSCENE  =  (OFFSCOC  +  EXPECTI4.I.J)  ♦  EXPECT (5, I, J>) /FSIM 

AONSCENE  =  NCYCLES+  ( INITUE+ ITOTRES  ( J ) )  -  OFFSCOE 

SORTYAC=NCYCLES+SORTYDAY/AONSCEJC 

WRITE  (6.9009)  I.EXPECTU, I, J)/FSIH. SORTYDAY, 

St  SORTYAC,  (EXPECT(H,  I,J)/FSIM.M=2,4> 

TOTFLY  =  TOTRY  +  SORTYDAY 

C— 

C—  +WRITE  THE  MEAN  SORTIES  PER  DAY  (FOR  EACH  DAY)  TO  A  FILE  THAT 
C—  COULD  BE  USED  BY  'CALLPtT2'  TO  PRODUCE  A  PLOT. 

WRITE  (7)  SORTYAC,  SORTYDAY 
WRITE(6,9003)  TOTRY 
6000  CONTINUE 


WRITE  (6.9010)  TOTFLY 


C— 

RETURN 

9002  FORMAT(V) 

9003  FORMAT  (F26.ll 

9005  FORMAT  ( ' 1 ' /// 10X. ^SORTIES/  SORTIES/  SORTIES/'. 21 I. 

It  'CUM.  RES.'/ 

t>  '  DAY  PER  PERIOD'. 5X. 'DAY', 8X»' AC' »4X» 

k  '  MCH  NMCS  LOSSES  LEFT'//) 

9006  FORMAT  ('  '.2X,I5,F9.1.F29.1.F9.1,F8.1) 

9008  FORMAT  ('  ', 12.15, F9.1,F29.1,F9.1,F8.I,F8.1) 

9009  FORMAT  ('  ",2X, I5.F9. 1.F9. 1.F11.2.2F9. 1.2F8. 1 ) 

9010  FORMAT  <//'  TOTAL  SORTIES  FLOWN  =',F11.1) 

EM) 


0 


SUBROUTINE  PRTHEP(RTIHE,PBAKSEO.D*>ne.NORSVC> 

CHHHWWiWHIWHHWHiHtHHHHHHHIWHHWIWHWHfH) 

C++  PRTREP  -  SIMULATES  PROCESS  OF  REPAIRING  PARTS. 

Cut  PRTREP  IS  A  FORTRAN  SUBROUTINE  WHICH  SIMULATES  DE  PROCESS 
Cm*  OF  REPAIRING  PARTS.  PARTS  REPAIR  IS  ASSUKD  TO  BE 
CM*  EIPOCNTIAL  WITH  AN  INFINITE  NUWER  OF  SERVERS.  I.E.  NO  PART 
Cm*  ESER  HAS  TO  UAIT  TO  BEGIN  SERVICE.  IN  AUDITION  TO  REPAIRING 
CM*  PARTS.  THIS  FUNCTION  CALCULATES  THE  NEW  NUMBER  OF  NORS  AIRCRAFT 
CM*  REMAINING  AFTER  REPAIR  OF  TT£SE  PARTS.  IF  ANY  PREVIOUSLY  NORS 
Cm*  AIRCRAFT  ARE  READY  TO  GO  INTO  MAINTENANCE.  THIS  ROUTDE 
CM*  HILL  DISTRIBUTE  THEM  PROBABILISTICALLY  AMONG  THE  VARIOUS 
CM*  UORKCENTERS. 

Cm* 

Cm*  ItPUT  - 

CM*  RTIIC  -  AVERAGE  THE  (IN  HOURS)  TIESE  PARTS  HAVE  BEEN 
Cm*  IN  REPAIR  SINCE  T>€  LAST  TIIE  PARTS  REPAIR  NAS 

Cm*  SIMULATED. 

Cm*  PBRKSEQ  -  2-DIMENSIONAL  ARRAY  USED  TO  DETERMINE  TIE 
Cm*  DISTRIBUTION  OF  ABORTS  INTO  TIE  VARIOUS 

Cm*  UORKCENTERS. 

C*M  INPUT/OUTPUT  - 

CM*  NORSVC  -  NORS  AIRCRAFT  STATUS  VECTOR.  INDICATES  THOSE 
Cm*  AIRCRAFT  WHICH  ARE  NORS  DUE  TO  UNAVAILABLE  PARTS. 

Cm*  THE  FIRST  WORD.  NORSVC(l).  CONTAINS  THE  TOM. 

C**«  NUKER  OF  1-BITS  IN  THE  NORS  STATUS  VECTOR. 

CM*  ARRAY  IS  A  BIT  VECTOR  UITH  EACH  BIT  REPRESENTING 

Cm*  m  AIRCRAFT.  A  1-BIT  INDICATES  THE  AIRCRAFT  IS 

C**«  STILL  FLYABLE.  NOTE  THAT  IFLYVC(l)  ALSO  IWICATES 

Cm*  TIE  NUWER  OF  1-BITS  IN  THIS  BIT  VECTOR. 

Cm*  COmON  IIEUT  - 

Cm*  DEPART  -  LOGICAL  FLAG  IWICATING  HETHER  THE  DFINITE  PARTS 
Cm*  ASSUETION  HOLDS.  IF  INFPART  IS  TRUE  TIEN  HERE 

CM*  IS  NEVER  ANY  SHORTAGE  OF  PARTS;  FENCE.  NO 

Cm*  NORS  AC. 

Cm*  IEARTS  -  NUMBER  OF  PART  TYPES  BEING  MODELED. 

Cm*  1QPA(K)  -  NUMBER  OF  TYPE-K  PARTS  INSTALLED  ON  EACH  AIRCRAFT. 

Cm*  RPRATE(K)  -  REPAIR  RATE  (PARTS/HOUR)  FOR  TYPE-K  PARTS 

C*m  INITSJ(K)  -  INITIAL  BASE  STOCK  LEVEL  FOR  TYPE-K  PARTS. 

CM*  COTTON  INPUT/OUTPUT  - 

Cm*  NBACKO(K)  -  NUWER  OF  BACKORDERS  FOR  PARTS  OF  TYPE-K.  IF 
Cm*  NBACKO(K)  IS  POSITIVE.  DEN  UNFULFILLED  REOESTS 

Cm*  FOR  PARTS  OF  THIS  TYPE  HAVE  BEEN  MADE.  IF  IT  IS 

Cm*  NEGATIVE.  THEN  WACKO(K)  INDICATES  THE  NUWER  OF 

Cm*  OF  PARTS  ON-THE-SFELF. 


C— 

PARAMETER  HAXPRT*304,  HAXCYOIO 
COMMON  /PARTS/  NPARTS.  IflPA(HAXPRT) .NBACKOi HAXPRT). 
t  BRPRATE  ( HAXPRT ) .  DRPRATE  ( MAXPRT ) .  INI  TSJ  ( HAXPRT ) .  RESUPPl  HAXPRT ) 

i  BNRTS  ( MAXPRT ) .  WASE  ( MAXPRT ) .  WEPOT  ( HAXPRT ) 

COMMON  /IIEUT/  INITUE.  NAC.  PATTRIT.  IRES.  RMCM.  DEPART. 

HAXFLY(MAXCYC).  DEMAN.  ISCALE.  IAUGMNT 


It 


or>  oo  oo  o  o  o  o  o  o  o  o  o  n  n  o  n  on  o  o  o  o  o  o  o  o  o 


DIMENSION  NORSVC(l) 
LOGICAL  INFPART 


—  ♦INITIALIZE  fCW  NlfCER  OF  NORS  AIRCRAFT  TO  MM 

fCUNOR  »  0 

—  *IFI  INFINITE  PARTS  NOT  ASSIMD)THEN 

IF (Iff PART)  GO  TO  5000 

—  *D0  FORIEACH  PART  TYPE) 

DO  3000  K*1.W>ARTS 

—  ♦IF(T)€R£  ARE  ANY  OF  THIS  PART  IN  REPAIR  >T)€N 

INSHPK  *  tCACKO(K)  ♦  INITSJ(K) 

IF(INSHPK.LE.O)  GO  TO  2000 

—  ♦DETERMINE  NUMBER  OF  THESE  WHICH  ARE  NEW  DEMANDS 

NEW  =  INSHPK  -  (NEEPOT(K)  ♦  fCASE(K)) 

—  ♦PERFORM  BINOMIAL  DRAW  TO  DETERMINE  BASET DEPOT  SPLIT 

NENDEP*NB!NON(BNRTS(K)>NEW) 

NOEPOT(K)=NDEPOT(K)+NEWDEP 

NBASE(K)=NBASE(K)+(NE)HeWDEP) 

—  •COMPUTE  PROBABILITY  OF  REPAIR 

PDEP  =  1.0  -  EXP(-DRPRATE(K)tfTItC) 

PBSE  -  1.0  -  EXP<-BRPRATEIK)«RTIME) 

—  *  DETERMINE  NUCER  OF  PARTS  REPAIRED  BY  SMCLING 

—  FROM  THE  APPROPRIATE  BINOMIAL  DISTRIBUTION 
NUWEP  »  WINOM(PDEP,NDEPOT(K)1 

NUMBSE  =  NBINOM(PBSEifCASE(K)) 

—  HJPDATE  NUMBER  IN-SHOP  AND  BACKORDERED 

NBACKO(K)  *  NBACKO(K)  -  (NEIOEP  ♦  NUMBSE) 

NDEPOT(K)  =  WEPOT(K)  -  NUMDEP 
NBASE(K)  *  NBASE(K)  -  NUMBSE 

—  ♦IFITTESE  PARTS  CAUSE  THE  MAX  NUMBER  OF  MORS  THUS  FAR) 

IFUCWNORtlOPA(K)  .GE.  NBACKO(K))  GO  TO  1000 

—  ♦UPDATE  NUMBER  OF  NORS  AIRCRAFT 

NEWNOR»INT(FLOAT(«ACKO<K))/FLOAT(IQPA(K))  ♦  .999) 

—  t€ND  IF  (£W  NORS  MAXIMUM  TEST) 

1000  CONTINUE 

—  ♦€»  IF  (NO  PARTS  IN  REPAIR  TEST) 

2000  CONTINUE 

—  ♦€»  DO  (PARTS  LOOP) 

3000  CONTINUE 

— ■  ♦ IF (ANY  PREVIOUSLY  NORS  AIRCRAFT  ARE  READY  FOR  REPA1R)T>€N 
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i 


NORDIF  *  NQRSVCU)  -  NEMNOR 
IF(NORDIF.Lf.O)  GO  TO  4000 

C— 

C —  tSELECT  LEFTMOST  MORS  AIRCRAFT  TO  ENTER  MAINTENANCE 

CALL  MCOIST(NOROIF.PSRKSEQ.INDXNC.NORSVC) 

C — 

C —  100  IF  <tCU  NONNORS  AC  TEST) 

4000  CONTINUE 
C — 

C —  *00  IF  ( INFINITE  PARTS  TEST/ 

5000  CONTINUE 

C — 

RETURN 

00 
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SUBROUTINE  R£PAIR(TIIET,NWC,NCi<£WS,SRATE) 


C++  REPAIR  -  SIMULATES  PROCESS  OF  WORK  CENTER  AIRCRAFT  REPAIR. 

Cm  REPAIR  IS  A  FORTRAN  ROUTINE  WHICH  SIMULATES  TIE 
Cm  PROCESS  OF  REPAIRING  AIRCRAFT  IN  ALL  WORKCENTERS  DURING  A 
Cm  SPECIFIED  TIME  PERIOD.  A  NUMBER  OF  REPAIRED  AIRCRAFT  IS  GEIERATED 
Cm  FOR  EACH  WORK  CENTER,  BY  SAMPLING  FROM  THE  APPROPRIATE  PROBABILITY 
Cm  DISTRIBUTIONS.  THIS  ROUTIIE  SIIELY  UPDATES  TIE  TOTAL  NO.  OF 
Cm  AIRCRAFT  IN  REPAIR  IN  EACH  WORKCENTER?  IT  DOES  NOT  CONCERN 
Cm  ITSELF  WITH  WHICH  PARTICULAR  AIRCRAFT  IN  A  W/C  ARE 
Cm  BEING  REPAIRED,  NOR  TIE  TOT#.  NO.  OF  AIRCRAFT  IN  MAINTENANCE. 

Cm  IMPLICITLY,  IT  IS  ASSUMED  THAT  WE  ARE  REPAIRING  THE  RIGHTMOST 
Cm  AIRCRAFT  ON  A  LIST  OF  AIRCRAFT  THAT  FEED  WORK  IN  A  GIVEN  W/C, 

Cm  AND  JWT  IF  AIRCRAFT  ARE  PLACED  ON  THAT  LIST  (IN  ROUTIIE 
Cm  'WCDIST')  IN  A  RANDOM  ORDER,  THEN  THIS  IETHOD  OF  REPAIR  IS 
Cm  ALSO  RAMIOH?  I.E.,  IT  DOESN'T  FAVOR  LOW-NUMERED  A/C,  OR 
Cm  HIGH-NUMBERED  A/C,  OR  RECENTLY-BROKEN  A/C,  ETC. 

Cm 

Cm  IIEUT  - 

Cm  TIMET  -  LENGTH  (IN  HOURS)  OF  TIE  REPAIR  PERIOD. 

Cm  NWC  -  TOTAL  NUWER  OF  WORKCENTERS. 

Cm  NCREWS  -  NUMBER  OF  CREWS  IN  WQRKCENTER-J. 

Cm  SRATE(J)  -  REPAIR  RATE  (AIRCRAFT/HOUR)  FOR  EACH  CREW  IN 

Cm  WORKCENTER-J. 

Cm  COMON  INPUTS/OUTPUTS  - 

Cm  INREPR(J)  -  NO.  OF  A/C  IN  MAINTENANCE  IN  W/C  J. 


C— 

PARAMETER  HAXWO25,MAXBIT=36,HAXAC»10e,MAXVEt=2+(IWXAC-l )/NAXBIT 
PARAMETER  HAXDAY=30 

PARAMETER  LFLD=7,NPERWRD=MAXBIT/LFLD,MXINWC=1+(HAXAC-1)/NPERWRD 
DIMENSION  NCREWS  (NWC),  SRATE(NWC) 

COmON  /WCMAINT/  LISTRP(MXINWC,MAXWC),IMEPR(MAXWC) 

C — 

C —  +DO  FOR (EACH  WORKCENTER) 

DO  400  J=1,NWC 

C— 

C —  »IF( THERE  ARE  ANY  AIRCRAFT  IN  REPAIR  IN  THIS  WORKCENTER) 

NACINJ  *  INREPR(J) 

IF(NACINJ.EQ.O)  GO  TO  200 

C— 

C —  +GOERATE  SAMPLE  NUMBER  OF  AIRCRAFT  REPAIRED 

NFIXED  -  NREPS(TIfET,NACINU,NCREWS(J),SRATE(J)) 

C— 

C—  tUPDATE  NO.  OF  A/C  IN  THIS  W/C 

INREPR(J)  -  NACINJ  -  NFIXED 

C— 

C—  «0B  IF  (ZERO  AIRCRAFT  IN  REPAIR  TEST) 

200  CONTINUE 

C— 

C—  *€»  DO  (WORKCENTER  LOOP) 

400  CONTINUE 


RETURN 


1 
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SUBROUTINE  SIMULA 

CnHHininimmm»HHinniiiim«n»iiinni<iiunnmiHii 

C++  SIMULA  -  PERFORM  SIMULATION  REPLICATIONS. 

C++*  THIS  ROUTINE  PROVIDES  THE  BASIC  STRUCTURE  OF  THE  SGM 
C++*  SIMULATION.  THE  LOOPS  WHICH  CONTROL  TIE  SIMULATION  REPLICATIONS. 
C++*  TIE  FLYING  DAYS  FOR  EACH  REPLICATION.  AM)  THE  FLYING  CYCLES 
Cm  FOR  EACH  FLYING  DAY  ARE  CONTAINED  IN  THIS  ROUTINE.  TIE 
Cm  ROUT IIC  BASICALLY  JUST  EXECUTES  A  SPECIFIED  NUTffiER  OF 
Cm  FLYING  CYCLES  EACH  FLYING  DAY. 

Cm  THE  ROUTIIE  READS  A  SUBSET  OF  SCENARIO  PARAMETERS  FOR  EACH 
Cm  FLYING  DAY  OF  TIC  SCENARIO  FROM  A  SCRATCH  FILE  (03)  WHICH 
Cm  MAS  INITIALIZED  IN  THE  INITSCN  SUBROUTINE.  THIS  APPROACH  ALLOWS 
Cm  THESE  PARAfCTERS  TO  VMY  ON  A  DAILY  BASIS  DURING  THE 
Cm  SIMULATION. 

Cm 

Cm  COMMON  HPUTS  - 

Cm  NSIM  -  NUMBER  OF  SIMULATION  REPLICATIONS  TO  BE  PERFORMED 

Cm  NUMOAY  -  NUMBER  OF  FLYING  DAYS  TO  SIMULATE 

Cm  NCYCLES  -  NUMBER  OF  FLYING  CYCLES  FOR  EACH  DAY.  THIS  PARAMETER 

Cm  IS  READ  FROM  THE  SCRATCH  FILE  EACH  SIMULATION  DAY 

Cm  COfttON  OUTPUTS  - 


Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 


ISIH  -  <ISIH=1,...,NSIM)  CURRENT  REPLICATION  NUMBER 
IDAY  -  (IDAY=l>..iNUMDAY)  CURRENT  FLYING  DAY 
ICYCLE  -  (ICYCLE=1,.... NCYCLES)  CURRENT  FLYING  CYCLE 
NRESRV  -  CURRENT  NUTffiER  OF  AIRCRAFT  IN  THE  RESERVE  POOL 
ITOTRES(I)-  <1=0, ....NUHLAY)  CUMULATIVE  NUMBER  OF  AVAILABLE 
RESERVE  AIRCRAFT  UP  TO  INCLUDING  THE  ITH  DAY. 

M  OTH  DAY  REPRESENTS  THE  INITIAL  NUMBER  OF  RESERVE 
AIRCRAFT.  THIS  ARRAY  IS  USED  IN  COTWTING  ON-THE- 
SCETC  AIRCRAFT  FOR  SORTIES/AIRCRAFT/DAY  IN  T)€ 
PRINTO  ROUTINE. 

IRES  -  NUTffiER  OF  RESERVE  AIRCRAFT  WHICH  BECOME  AVAILABLE 
ON  THIS  FLYING  DAY. 

PATTRIT.PACGABT.WAITCYC, TYIKITE,  IWXFLY( I ) 

-  SCENARIO  PARATCTERS  WHICH  ARE  MJ.0UED  TO  VWY 
ON  A  DAILY  BASIS.  TFC.Y  ARE  LOADED  IN  THIS 
ROUTITE  AND  SUPPLIED  TO  THE  FLYC'YC  ROUTINE 


CiiHiiiiHHiHHiHHHWiHWHiiniiiiiiiiimimmiimmiu 


C — 


it 


& 

& 

It 

C— 

c— 


PARAMETER  NAXWC=25,  HAXDAY=30*  NUCYC=10.  NAXSTAT=5 
COMMON  /INPUT/  INITUE.  NAC,  PATTRIT,  IRES.  RNMCM.  ITFPART, 
MAXFLY(MAXCYC),  INFHAN,  I SCALE,  IAUGHNT 
COMNON  /STATS/  EXPECT  <  HAXSTAT , HAXCYC , HAXDAY ) , 

NRESRV,  I ZDAY , I TOTRES ( MAXDAY ) ,  LOSSTOT 
COMMON  /TIME/  PRERITE,  SORTLGTH,  WAITCYC,  TYTfdTE, 

NSIM,  ISIH,  NUMDAY,  IDAY,  NCYCLES,  ICYCLE 
COMMON  /WC8RK/  PACBRK,  PACGABT,  PBRKWC(MAXWC),  PWCPROD, 
PbRKSE0(2,MAXWC),  INDXWC(MAXWC) 


*00  FOR(EACH  SIMULATION  REALIZATION) 
DO  300  ISIH=1»NSIH 


C— 


B-58 


oo  o<->  oo  <_>  o  oo  o  o  o  u  oo  o  o 


C —  *INITIALIZE  VARIABLES  AT  START  OF  EACH  REPLICATION 

CALL  INITREP 

—  *REWIND  DAILY  SCEMRIO  PARAMETERS  FILE 

REWIND  03 

—  *00  FORIEACH  FLYING  DAY) 

DO  200  IDAY=1.NUHDAY 

—  *READ  DAILY  SCENARIO  PARAMETERS 

READ103)  PATTRIT.PACGABT.NCYCLES.  IR£S,UAITCYC.TYt*UTE 
READ(03)  <HAXFLY(J),J=1,NCYCLES) 

—  *UPDATE  RESERVE  AC  POOL  WITH  NEW  AVAILABLE  RESERVES 

fUESRV  *  NRESRV  *  IRES 
ITOTRES(IDAY)  =  IT0TRES< IDAY-1 )  *  IRES 

—  *00  FORIEACH  FLYING  CYCLE) 

DO  100  I CYCLE5 1 iNCYCLES 

—  *SIMUIATE  A  RYING  CYCLE 

CALL  FLYCYC 

—  *£ND  DO  (CYCLE  LOOP) 

100  CONTINUE 

—  *04)  DO  (DAY  LOOP) 

200  CONTINUE 

—  «END  DO  (REPLICATION  LOOP) 

300  CONTINUE 

C — 

RETURN 

END 
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SUBROUTINE  SGRTDSiNRECS.RKEYS, INDEX) 


C++  SORTDS  -  DESCENDING  SORT  OF  A  REAL  ARRAY. 

C«*  THIS  ROUTINE  RETURNS  A  SORTED  LIST  OF  INDICES  ACCORDING 
Cm  TO  T>€  GIVEN  ARRAY  OF  FLOATING-POINT  KEYS.  THIS  IS  A 
Cm  DESCENDING  SORT  II.E.  LARGEST  TO  SMALLEST)  USING  THE 
Cm  SIMPLE  EXCHANGE  SORT  TECHNIQUE. 


Cm 

Cm 

INPUTS  - 

Cm 

uses 

Cm 

Cm 

RKEYS 

Cm 

Cm 

OUTPUTS  - 

Cm 

INDEX 

Cm 

Cm 

Cm 

-  NUMBER  OF  ENTRIES  OR  RECORDS  IN  THE  INPUT  AND 
OUTPUT  ARRAYS. 

-  ARRAY  OF  FLOATING-POINT  KEYS  WHOSE  CORRESPONDING 
INDICES  ARE  TO  BE  SORTED. 

-  (WRAY  OF  INDICES,  1,...,NRECS  ,  SORTED  IN  DESCENDING 
ORDER  ACCORDING  TO  THE  CORRESPONDING  RKEY  ENTRY. 
THUS,  RKEY(  IffflEXd)  )  IS  THE  LARGEST  ENTRY,  AND 
RKEYf  INDEX(NRECS)  )  IS  THE  SMALLEST. 


C — 

DIMENSION  RKEYS(NkECS),  IMOYNRECS) 


C— 

C —  +1NITIALIZE  INDEX  ARRAY  TO  NATURAL  ORDER 
DO  100  1=1,IW£CS 
INDEX ( I )= J 
100  CONTINUE 


C— 


C —  »IF( THERE  IS  MORE  THAN  1  RECORD) THEN 
IF(NRECS.LE.1)G0  TO  500 


C— 

C—  »D0  UNTILXNO  INTERCHANGES  OCCUR) 

LIMIT=fRECS-l 
200  CONTINUE 

C— 

C —  ^INITIALIZE  INTERCHANGE  FLAG  TO  NONE 

LASTX=0 

C — 

C —  *D0  F0R1EACH  CNSECUTIVE  PAIR  OF  RECORDS) 

DO  400  J=l, LIMIT 

C— 

C —  *IF<  INDICES  OF  THESE  RECORDS  SHOIU)  BE  INTERCHANGED) 

IFIRST  =  INDEX (J) 

I SEC  =  INDEX (J+l) 

IF(RKEYSdFIRST)  .GE.RKEYS!  ISEC) )  GO  TO  300 

C — 

C —  +INTERCHANGE  IWJICES  AND  NWK  LAST  EXCHANGE  LOCATION 

INDEX! J)  =  ISEC 
INDEX! J+l)  =  IFIRST 
LASTX  =  J 

C— 

C —  *END  IF  (INTERCHANGE  TEST) 

300  CONTINUE 
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C—  *END  DO  (RECORD  PAIR  LOOP) 

400  CONTINUE 

♦END  DO  (INTERCHANGE  PASS  LOOP) 
LIMIT  =  LASTX 
IF(LASTX.NE.O)  GO  TO  200 

C — 

C—  *END  IF  (SINGLE  RECORD  TEST) 

500  CONTINUE 
C — 

RETURN 

EM) 


1 
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SUBROUTINE  TBITSLINONES.  IFROM,  ITO) 

CHiinniHiHHHUinmiuiiumtiMmiminmiKimimiiHi 

C++  TBITSl  -  TRANSFER  1-BITS  FROM  LEF  OF  A  BIT-VECTOR. 

C***  TBITSL  IS  A  FORTRAN  SUBROUTINE  WHICH  WILL  TRANSFER 
C++*  A  SPECIFIED  NUMBER  OF  1-BITS  IN  TIE  LEFTMOST  PORTION  OF 
C+++  A  BIT-VECTOR.  'IFROM',  TO  THE  CORRESPONDING  POSITIONS 
C*M  OF  A  BIT-VECTOR,  ITO".  TIE  TRANSFERRED  BITS  ME  THEN 
C++*  ZEROED  OUT  IN  'IFROH'.  THESE  BIT-VECTORS  ME  KEPT  IN 
C**+  ARRAYS,  ORGANIZED  IN  THE  FOLLOWING  MANNER  -  THE  FIRST  WORD 
C«*  OF  TIE  ARRAY  CONTAINS  THE  CURRENT  NUMBER  OF  1-BITS  IN  THE 
Cm*  BIT-VECTOR,  AM)  THE  REMAINING  WORDS  OF  THE  ARRAY  CONTAIN  THE 
C«*  ACTUAL  BIT-VECTOR. 

Cm-* 


C+++  INPUT  - 
C+«  NONES 
C*** 

C*+* 

C»*  INPUT/OUTPUT 
C***  I FROM 
C*+* 

C*** 

C*** 

C*« 

C*m 

C*M 

C*M 

C*M 

C*M  ITO 

C*M 

C*M 

C*M 

C*** 

C*M 

C*M 

c*« 

C*M  CONMON  INPUT  ■ 
CM*  length 

C+M 


NUMBER  OF  1-BITS  TO  BE  TRANSFERRED.  NOTE  THAT 
NONES  MUST  NOT  BE  GREATER  T>W  TIE  NUMBER  OF 
IS  IN  TIE  BIT  STRING,  'IFROH. 

ARRAY  CONTAINING  TIE  BIT-VECTOR  FROM  WHICH 
THE  LEFTMOST  1-BIT  POSITIONS  ME  TO  BE 
TRANSFERRED.  NOTE  THAT  THESE  LEFTMOST 
1-BITS  ARE  ZEROED  OUT  IN  'IFROM'  UPON 
COMPLETION  OF  THIS  ROUTINE.  IFROM(l)  IS 
A  COUNTER  WHICH  INDICATES  THE  CURRENT 
NUMBER  OF  1-BITS  IN  THE  BIT-VECTOR.  THE 
ACTUAL  BIT-VECTOR  IS  THE  CONSECUTIVE  BITS 
CONTAINED  IN  THE  WORDS  IFRUM(2)-IFR0M( LENGTH) 
ARRAY  CONTAINING  THE  BIT-VECTOR  TO  WHICH 
THE  LEFTMOST  1-BITS  IN  'IFROM'  ARE  TO  BE 
TRANSFERRED.  TIE  RIGHTMOST  POSITIONS  IN  ITO 
REMAIN  AS  BEFORE.  ITO(l)  IS  A  COUNTER  WHICH 
INDICATES  TIE  CURRENT  NUMBER  OF  i-BITS  IN 
THE  BIT-VECTOR.  THE  ACTUAL  BIT  VECTOR  IS  THE 
CONSECUTIVE  BITS  CONTAINED  IN  TNE  WORDS  - 
ITO<2)  -  ITO(LENGTH) 

LENGTH  (IN  COMPUTER  WORDS)  OF  THE  ARRAYS 
CONTAINING  TIE  VMIOUS  BIT-VECTORS!  IFLYVt,  ETC 


C— 


PMANETER  NAXAC=10e,MAXBIT=36,HAXVEC=2+(HAXAC-l  >/HAXBIT 
COrtlON  /ACSTATE/  LENGTH, NACVCl MAX VEC),  IFIYVC(MXVEC). 

&  MAINVC(MAXVEC),NORSVC(MAXVEC),  LOSTVC(NWXVEC) 

DINENSION  IFROM! 1),  ITO(l) 

C — 

C—  ♦INITIALIZE  TRANSFER  LOOP 
NLEFT  =  NONES 
INDEX  =  2 

C— 

C —  *00  WHILEIALL  APPROPRIATE  WORDS  MVE  NOT  BEEN  MODIFIED) 

1000  CONTINUE 

IF(NLEFT.LE.O)  GO  TO  4000 
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IFdNDEX.GT.LENGTH)  GO  TO  4000 

C— 

C—  ♦COUNT  NUMBER  (if  IS  IN  t£XT  UORD 

NXT1S  =  N1BITS<  IFRCMdWCX)  ) 

C— 

C —  ♦IF(NOT  ALL  1-BITS  IN  THIS  WORD  SHOULD  BE  TRANSFERRED) 

IF(NXTlS.LE.MfFT)  GO  TO  2000 

C — 

C —  ♦TRANSFER  TO  PROPER  NUMBER  OF  IS 

L TRANS  =  LBITS<  IFROHdWEX)  ,NLEFT) 

I  T0(  INDEX )  =  OR(ITOdNDEX).LTRANS) 

IFROHt INDEX)  =  XOR<  IFROMdNDEX)  .LTRANS) 

C — 

C —  ♦ELSE  (MJL  1-BITS  IN  THIS  WORD  ARE  TO  BE  TRANSFERRED) 

GO  TO  3000 

2000  CONTINUE 

C — 

C —  ♦TRANSFER  TO  ENTIRE  UORD 

ITO(  INDEX )  =  GR(ITO(IW)EX>,IFRCH(INDEX)) 

IFROMdNDEX)  =  0 

C— 

C—  «£ND  IF  (Ml  IS  TEST) 

3000  CONTINUE 

C— 

C —  ♦UPDATE  NUMBER  OF  IS  LEFT  TO  TRANSFER 

fLEFT  =  N.EFT  -  NXT1S 
C— 

C—  ♦  INCREMENT  INDEX  FOR  NEXT  WORD  OF  BIT-VECTOR 

INDEX  =  INDEX  ♦  1 

C— 

C —  *END  DO  (UORD  LOOP) 

GO  TO  1000 
4000  CONTINUE 
C— 

C —  ♦UPDATE  IS  COUNTER  FOR  THESE  BIT-VECTORS 

1T0(1)  =  ITO(l)  +  NONES 

IFROH(l)  =  IFROnil)  -  NOTES 

C — 

C—  ♦PERFORM  ERROR  0£CX  TO  ENSURE  APPROPRIATE  NUMBER 
C— -  OF  1-BITS  HAS  TRANSFERRED 

IF(TEEFT.GT.O)  HR I TE (0,9001) NONES. I FROM (1),NLEFT 

C— 

RETURN 

9001  FORMAT (“0*«f*«*tt  TBITSL  ERROR  -  TOO  FEU  1-BITS  TO  TRANSFER" • 
V  /,  •  tUtitii  NONES,  IFROM(l),  NLEFT  =  ’,315) 

END 


SUBROUTINE  TBITSRtlMES.  I  FROM.  ITO) 

CHHii>n»mumnniiiiiiiiiiiWHwiniinnimiinii[inini> 

C++  TBITSR  -  TRANSFER  1-BITS  FROM  RIGHT  OF  A  BIT-VECTOR. 

C**+  TBITSR  IS  A  FORTRAN  SUBROUTIIC  WHICH  HILL  TRANSFER 
C***  A  SPECIFIED  NUMBER  OF  1-BITS  IN  THE  RIGHTMOST  PORTION  OF 
C+*+  A  BIT-VECTOR.  'IFROM',  TO  THE  CORRESPOMJINO  POSITIONS 
C+»*  OF  A  BIT-VECTOR,  /ITO/.  THE  TRANSFERRED  BITS  ARE  TTCN 
C**+  ZEROED  OUT  IN  'IFROM'.  THESE  BIT-VECTORS  ARE  KEPT  IN 
C**+  ARRAYS,  ORGANIZED  IN  THE  FOLLOWING  MANNER  -  THE  FIRST  WORD 
C++*  OF  TIC  ARRAY  CONTAINS  THE  CURRENT  NUMBER  OF  1-BITS  IN  THE 
C**+  BIT-VECTOR,  AND  TTC  REMAINING  WORDS  OF  THE  ARRAY  CONTAIN  THE 
C***  ACTltt.  BIT-VECTOR. 

C++* 


C++*  INPUT  - 
C++*  NONES 
C++* 

c+*+ 

C**+  INPUT/OUTPUT  - 
C++*  IFROM 
C++* 

c**+ 

C**+ 

C++* 

C++* 

C**+ 

c+*+ 

c**+ 

C**+  ITO 
C**+ 

c+*+ 

C++* 

c+*+ 
c**+ 

C++* 

C++* 

c*++  COMMON  INPUT 
C++*  LENGTH 
C+*+ 


-  NUMBER  OF  1-BITS  TO  BE  TRANSFERRED.  NOTE  THAT 
NONES  MUST  NOT  BE  GREATER  THAN  THE  NUMBER  OF 
l'S  IN  THE  BIT  STRING,  'IFROM. 

-  ARRAY  CONTAINING  THE  BIT-VECTOR  FROM  WHICH 
TTC  RIGHTMOST  1-BIT  POSITIONS  ARE  TO  BE 
TRANSFERRED.  NOTE  THAT  THESE  RIGHTMOST 
1-BITS  ARE  ZEROED  OUT  IN  IFROM'  UPON 
COffLETION  If  THIS  ROUTINE.  IFROH(l)  IS 

A  COUNTER  WHICH  INDICATES  THE  CURRENT 
NUMBER  OF  1-BITS  IN  THE  BIT-VECTOR.  THE 
ACTUAL  BIT-VECTOR  IS  THE  CONSECUTIVE  BITS 
CONTAINED  IN  TTC  WORDS  IFROM(2)-IFROM<LENGTH) 

-  ARRAY  CONTAINING  THE  BIT-VECTOR  TO  WHICH 
TIC  RIGHTMOST  1-BITS  IN  'IFROM'  ARE  TO  BE 
TRANSFERRED.  TIC  LEFTMOST  POSITIONS  IN  ITO 
REMAIN  AS  BEFORE.  ITO(l)  IS  A  COUNTER  WHICH 
INDICATES  THE  CURRENT  NtffflER  OF  1-BITS  IN 
THE  BIT-VECTOR.  THE  ACTUAL  BIT  VECTOR  IS  THE 
CONSECUTIVE  BITS  CONTAINED  IN  THE  WORDS 

IT0<2)  -  ITO(LENGTH) 


-  LENGTH  (IN  COMPUTER  WORDS)  OF  THE  ARRAYS 
CONTAINING  THE  VARIOUS  BIT-VECTORS?  iaYVC,  ETC 

CHHUimuiiiiii[nminininiim»mnnmmwmiiHiiin4 

C — 


PARAMETER  NAXAC=108,HAXBIT=36,NAXVEC=2+(HAXAC-l)/)ttXBIT 
COMMON  /ACSTATE/  LENGTH, NACVCIMAXVEC),  IFLYVC(MAXVEC), 
i  MAINVC(MAXVEC).NORSVC(MAXVEC),  LOSTVC(fttXVEC) 

DlfCNSION  IFROH(l),  ITO(l) 


C—  UNITIALIZE  TRANSFER  LOOP 
NLEFT  =  NOtCS 
INDEX  =  LENGTH 

C — 

C—  +DO  WHIL£(ALL  APPROPRIATE  WORDS  HAVE  NOT  BEEN  MODIFIED) 
1000  CONTINUE 

IF(NLEFT.LE.O)  GO  TO  4000 
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IF (INDEX  ,L£.  I)  GO  TO  4000 

C— 

C —  *COKT  NUWER  OF  1'S  IN  NEXT  WORD 

NXT1S  *  MBITS)  IFROM(IWEX)  ) 

C — 

C —  *IF(NOT  ALL  1-SITS  IN  THIS  WORD  SHOULD  BE  TRANSFERRED) 

IF(NXTIS.LE.NLEFT)  GO  TO  2000 

C — 

C —  *  TRANSFER  THE  PROPER  NUMBER  OF  I  S 

NTRANS  =  LBITS1  I  FROM  INDEX)  .NXTIS-ftEFT) 
ITOUNOEX)  *  OR(ITO< INDEX ),X0R(IFRQM( INDEX). NTRANS)) 
IFROM(IWEX)  -  NTRANS 

C — 

C —  *€LSE  (ALL  1-BITS  IN  THIS  WORD  ARE  TO  HE  TRANSFERRED) 

GO  TO  3000 

2000  CONTINUE 

C — 

C —  * TRANSFER  T)C  ENTIRE  WORD 

ITO< INDEX)  =  OR  ( ITO( INDEX ) , IFROM( INDEX) ) 

IFROH< INDEX)  =  0 

C— 

C —  »0C  IF  (ALL  TS  TEST) 

3000  CONTINUE 

C— 

C—  •UPDATE  NUMBER  OF  TS  LEFT  TO  TRANSFER 

N-LFT  =  HEFT  -  NXT1S 

C— 

C —  * DECREMENT  INDEX  FOR  NEXT  WORD  OF  BIT-VECTOR 

index  =  ihex  - 1 

c— 

c —  *END  DO  (WORD  LOOP) 

GO  TO  1000 
4000  CONTINUE 

C — 

C—  »UPDATE  I  S  COUNTER  FOR  TFCSE  BIT-VECTORS 
ITO(l)  =  ITO(l)  ♦  NOHS 
IFROH(l)  =  IFROM(l)  -  NMCS 

C— 

C—  •PERFORM  ERROR  CHECK  TO  ENSURE  APPROPRIATE  NUMBER 
C —  OF  1-BITS  HAS  TRANSFERRED 

IF(HEFT.GT.O)  WRITER, 9001)NONES.IFROM(1>,NL£FT 

C — 

RETURN 

9001  FORMAT)  "Om***M  TBITSR  ERROR  -  TOO  FEW  1-BITS  TO  TRANSFER". 
i>  /,  "  IHIIIH  NONES,  IFR0M(1>.  HEFT  =  ",3I5) 

END 
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SUBROUTINE  UEUPDAT(NAC) 


C++  UEUPOAT  -  UPDATE  UE-STRENGTH  FOR  SCENARIO. 

Cm  THIS  ROUTINE  IS  USED  TO  UPDATE  THE  UE-STRENGTH  OF  THE 
Cm  SCENARIO.  TIE  UE-STRENGTH  IS  INITIALIZED  AT  THE  START  OF 
Cm  THE  aYING  SCENARIO  AND  NORMALLY  DOES  NOT  CHANGE  THROUGHOUT 
Cm  THE  SIMULATION.  HOWEVER,  IF  THE  USER  HAS  SELECTED  TIE 
Cm  AUGMENTATION  MODE  FOR  RESERVE  AIRCRAFT,  AND  ENOUGH  RESERVES 
Cm  ARE  AVAILABLE  TO  MORE  THAN  REPLACE  COMBAT  LOSSES,  THEN  THE 
Cm  EXCESS  RESERVES  ARE  USED  TO  INCREASE  THE  CURRENT  UE-STRENGTH. 
Cm  INCREASING  THE  UE-STRENGTH  REQUIRES  RECALCULATION 
Cm  OF  TIC  LENGTH  OF  THE  AIRCRtfT-STATUS  BIT-VECTORS,  AM)  ALSO 
Cm  REINITIALIZATION  OF  THE  TOTAL-AIRCRAFT-POPULATION  VECTOR. 

Cm 

Cm  INPUT  - 

Cm  NAC  -  NEW  UPDATED  UE-STRENGTH. 

Cm  MAXAC  -  PARAMETER  IWICATING  MAXIMUM  ALLOWABLE 

Cm  AIRCRAFT  IN  THE  CURRENT  SGM  CONFIGURATION. 

Cm  MAXBIT  -  NUMBER  OF  BITS  PER  COMPUTER  WORD. 

Cm  COMMON  OUTPUT  - 

Cm  LENGTH  -  LENGTH  (IN  COMPUTER  WORDS)  OF  THE  AIRCRAFT 
Cm  BIT-VECTORS.  LENGTH  MUST  BE  LARGE  ENOUGH  SO  THAT 

Cm  THE  BIT-VECTORS  ARE  AT  LEAST  ‘NAC’  BITS 

Cm  LONG  AT  ’MAXBIT*  BITS  PER  COMPUTER  WORD. 

Cm  NAC(I)  -  (1*1,2,..., LENGTH)  BIT-VECTOR  WITH  FIRST  *NAC* 
Cm  BITS  SET  TO  1.  THIS  BIT-VECTOR  INDICATES  THE 

Cm  TOTAL  AIRCRAFT  POPULATION  FOR  THE  SCENARIO. 


C— 

PARAMETER  MAXBIT=36,  NAXAC=10e.  HAXVEC=2+(HAXAC-1> /MAXBIT 
COMMON  /BITS/HASK0,HASK(35).  HLEFT0,HSKLFT(36),IZC0UT, 

V  IC0UNT<63) 

COMMON  /ACSTATE/  LENGTH.  IWCVC I HAXVEC),  IFLYVC(MAXVEC), 
fc  HAINVC(MAXVEC),NORSVC(MttXVEC),  LOSTVC(MAXVEC) 

C — 

C —  (PERFORM  ERROR  CHECK  TO  ENSURE  NO  UE-OVERFLOW 
IF ( NAC. LE. MAXAC) GO  TO  100 
WRITE  U>,  9001  )NAC,  MAXAC 
NAC  *  MAXAC 
100  CONTINUE 
C — 

C —  (RECOMPUTE  LENGTH  OF  AIRCRAFT-STATUS  BIT-VECTORS 
LENGTH  =  2  ♦  (NAC-1) /MAXBIT 

C— 

C—  (INITIALIZE  TOTAL-AIRCRAFT-POPULATION  BIT-VECTOR 
CALL  SPRAY (MSKLFT( 36 ),NACVC( 2) .LENGTH- 1) 

NACVC(l£NG TH)*M$KLFT(H0D(NAC-1. MAXBIT) >+l 
NAC  VC  ( 1 )  =NAC 

C— 

RETURN 

9001  FORMAT COfW*m  UEUPDAT  EJPOR  -  UE  OVERFLOW*,/, 
t  *  mmt(  NAC,  MAXAC  *  *,215,/. 
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fc  •  $»$»$<**  me  HILL  BE  TRUNCATED  TO  HAXAC*) 

END 
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SUBROUTINE  HCDIST(NBRKAC,PBRKSE8> INUXUC, IACVC) 


C++  UCDIST  -  DETERMINE  BREAK  DISTRIBUTION  INTO  WORK  CENTERS. 

C++*  UCDIST  IS  A  FORTRAN  SUBROUTINE  UHICH  SIMULATES  THE  PROCESS 
C++*  OF  AIRCRAFT  BREAKING  INTO  UORKCENTERS.  GIVEN  A  NUMBER  OF  AIRCRAFT 
Cm  UHICH  MVE  BROKEN  INTO  AT  LEAST  ONE  UORKCENTER  -  1,2,...,NUC,  THIS 
Cm  ROUTINE  DETERMINES  (BY  SIMULATION)  UHICH  PART I CULM  UORKCENTERS 
C++*  TIC  AIRCRAFT  BROKE  INTO. 

Cm 


C++*  IHPUTS  — 

C++*  IACVC 

Cm 

Cm 

Cm 

Cm 

Cm 

C++* 

Cm  NBRKAC 

Cm 

Cm 

Cm 


-  BIT  VECTOR  INDICATING  AIRCRAFT  FROM  UHICH  BROKEN 
AIRCRAFT  ARE  TO  BE  SELECTED.  A  1-BIT 
INDICATES  AN  AIRCRAFT  UHICH  IS  A  CANDIDATE  FUR 

ONE  OF  TIC  BROKEN  AIRCRAFT.  THIS  ROUTINE  ARBITRARILY 
SELECTS  TFE  LEFTMOST  1-BITS  AS  THOSE  AIRCRAFT  UHICH 
WILL  BREAK  INTO  MAINTENANCE.  NORMALLY,  IACVC 
IS  'IRYVC'  OR  'NORSVC'. 

-  NUMBER  OF  BROKEN  AIRCRAFT  UHICH  ARE  TO  BE  BROKEN 
INTO  THE  DIFFERENT  UORKCENTERS.  THE  LEFTMOST 
'NBRKAC'  1-BITS  IN  'IACVC'  ARE  SELECTED  AS  THE 
AIRCRAFT  UHICH  BROKE. 


Cm  COMMON  IMVTS  - 

Cm  NIC  -  TOTAL  NUMBER  OF  UORKCENTERS 

Cm  HASK(I)  -  CONTAINS  A  1  IN  THE  ITH  BIT  (COUNTING  FROM  THE  LEFT) 

Cm  m  ZEROES  EVERYUHERE  ELSE.  1=0 . 35 

Cm  COtllON  INPUTS/OUTPUTS  - 

Cm  INREPR(J)  -  NUTffiER  OF  AIRCRAFT  IN  UORKCENTER- J. 

Cm  LISTRP(I.J)  -  LISTRP(  .  ,J)  IS  A  LIST  OF  AIRCRAFT  NUMBERS 
Cm  INDICATING  THOSE  AIRCRAFT  REQUIRING  MAINTENANCE  IN 

Cm  THE  JTH  HORK-CENTER  (J=1,2,...,NUC).  THIS  LIST 

Cm  CONTAINS  EXACTLY  IMREPR(J)  AIRCRAFT  NUfflERS.  TO  SAVE 

Cm  SPACE,  TTESE  LISTS  HAVE  BEEN  PACKED  INTO  BIT-FIELDS 

Cm  INS1EAD  OF  NORDS.  EACH  NUTfBER  IS  STORED  IN  A  BIT-FIELD 

Cm  "LFLD*  BITS  HIDE?  HENCE,  IF  ’HAXBIT"  IS  THE  LENGTH 

Cm  OF  A  COMPUTER  HURD  ON  THIS  SYSTEM,  THEN  THERE  ARE 

Cm  (MAXBIT/LFLD)  BIT-FIELDS  STORED  PER  UORD.  THE  AIRCRAFT 

Cm  NUMBERS  STORED  IN  TVESE  BIT-FIELDS  INDICATE  A  UNIQUE 

Cm  BIT-POSITION  IN  THE  VMIOUS  AIRCRAFT-STATUS  BIT- 

Cm  VECTORS.  THE  AIRCRAFT  ARE  NUMBERED, LEFT-TO-RIGH), 

Cm  0, 1,2,...,( MAXAC— 1 )  .  TO  GET  THE  ITH  AIRCRAFT  NUMBER 

Cm  IN  A  HORK-CENTER  LIST,  T)€  CORRESPONDING 

Cm  BIT-POSITION  AM)  HORD-INDEX  MUST  BE  COMPUTED. 

c— 

PARAMETER  M)XUC=25,fttXBIT=3Z>,MAXAC=i0e,MAXVEC=2+()«XAC-l  J/MAXBIT 
PARAftTER  LFLft=7,NPERURD=MAXBIT/LFU),HXINUC=l+(MAXAC-t)/NPERURD 
COMMON  /RSEED/  SEED 

COMMON  /UCINPUT/  NMC,  NCREUS(MAXUC),  SRATE(MAXUC) 

COMMON  /UCMAINT/  LISTRP(HXINMC,lttXHC>.  INREPR(HAXUC) 

COtWN  /BITS/  HASK0,MASK(35),  MLEFT0,MSKLFT(36), 
fc  IZCOUT,ICOUNT(S3> 

COTWN  /ACSTATE/  LENGTH, NACVC(MAXVEC),  IFLYVC(NAXVEC), 
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na  i  m;  ( haxvec ) ,  mm  (haiveu.  lqstvci haxvec) 
DIMENSION  PBRKSEQ(2,HAXWC),  INDXHC(l),  IACVC(l) 

♦IF (THERE  ARE  ANY  BRCKEN  AIRCRAFDTHEN 
IF(NBRKAC.EQ.O)  GO  TO  700 

♦INITIALIZE  NUMBER  OF  SELECTED  AIRCRAFT  TO  NONE 
NSELEC  =  0 

♦DO  FOR (EACH  WORD  OF  THE  IWUT  AIRCRAFT  VECTOR) 

DO  500  IHORD  =  2, LENGTH 

♦INITIALIZE  DO 
IACSIT  -  IACVCI IUORD) 

IGRKVC  «  0 

♦DO  FORI  EACH  BIT  OF  THIS  BIT-^CTOR  WORD) 

DO  400  IBIT  *  l.MAXBIT 

♦IF< THIS  BIT  INDICATES  AN  ELIGIBLE  AIRCRAFT)TI£N 
MASKAC  »  HASK(IBIT-l) 

IF(ANDUACBITiMASKAC)  .EQ.  0)  GO  TO  300 

♦SELECT  THIS  AIRCRAFT  AS  BROKEN 
IBRKVC  =  OR (IBRKVC, MASKAC) 

NSELEC  =  NSELEC  ♦  1 

*amJTE  AIRCAFT  * 

I AC  *  <IH0RD-2)«MAXBIT  ♦  (IBIT-1) 

♦DO  FORIEACH  WORKCENTER) 

DO  100  J  =  liNHC 

♦DRAW  RANDOM  SAMPLE  FROM  UNIFORM  (0,1) 

RDRAM  =  UNIFHl(SEED) 

♦CONTINUE  LOOP  WITH  NEXT  HORKCENTER  IF 
DRAU  INDICATES  NO  BREAK  INTO  THIS  UC 
IF(RDRAU  .GT.  PBRKSEQ(2,J))  GO  TO  100 

♦UPDATE  NUHBER/DISTRIBUTION  FOR  THIS  HC 
JREAL  =  IW)XWC(J) 

NTOREP  =  INREPR(JREAL)  +  1 

♦MAKE  'IRAND'  A  RAWKM  INTEGER  BETWEEN  1  AW) 
THE  NO.  OF  A/C  THAT  HILL  BE  IN  THIS  U/C 
I  RAW)  =  INT  ( UNIFHM  SEED  )*FLOAT<  NTOREP))  +  1 


♦MOVE  T)€  A/C  CURRENTLY  AT  SPOT  'I RAW)'  IN  THE 
LIST  TO  THE  RIGHTMOST  SPOT. 
FLD(H0D(NT0REP-1tWERVRD)HFLDiLFLD, 

LISTRPl  l+INTOREP-l  (/WERWRD.  JREAL) ) 

=  FLD(MOO(IRAWM,W)ERHRD)<HJLD.LFLD, 
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t.  LISTRP(  l+(  IRAKM  >/NPERWRD,  JREAL) ) 

—  *INSERT  THE  SELECTED  A/C  INTO  SPOT  'IRAK)' 

FLD(M0D(  IRAKH,NPERNRD)*LFLD,LFLD, 

LISTRP( l+( IRAWH ) /NPERURD, JREAL) ) 

=  I  AC 

—  ♦INCREMENT  THE  NO.  GF  A/C  IN  THIS  W/C 

INREPR(JREAL)  =  INREPR(JREAL)  +  1 

—  <€XIT  WGRKCENTER  LOOP  IF  DRAW  ALSO  INDICATES 

—  NO  BREAKS  INTO  REMAINING  WORKCENTERS 
IF(RDRAU  .LE.  PBRKSEQ(1,J))  GO  TO  200 

—  *£ND  DO  (WORKCENTER  LOOP) 

100  CONTINUE 

—  *****  ERROR  *****  IF  THIS  STATEMENT  IS  REACHED, 

—  THEN  EITHER  PBRKSEQU.NWC)  DOES 

—  NOT  EQUAL  1.0,  OR  THE 

—  RANDOM  DRAW  IS  GREATER  THAN  1.0 
WRITE<6,9001)PBRKSEQ(1,NWC),RDRAW 

—  «EXIT  DO  (WORKCENTER  LOOP)  -  THIS  IS  THE 

—  NORMAL  EXIT  FROM  TIC  HORKCENTER  LOOP 

200  CONTINUE 

—  *IF  ALL  T)C  BROKEN  AIRCRAFT  HAVE  BEEN  SELECTED 

—  THEN  ALL  LOOPS  ARE  TERMINATED 

IFtNSELEC  .GE.  NBRKAC)  GO  TO  600 

—  *£ND  IF  (FLYA8LE  AIRCRAFT  TEST) 

300  CONTINUE 

—  *END  DO  (BIT  LOOP) 

400  CONTINUE 

—  HJPDATE  IfFUT  BIT-VECTOR 

IACVC(IWORD)  =  X<JR<IACVC<  IWORDMBRKVC) 

—  ♦£»  DO  (WORD  LOOP) 

500  CONTINUE 

—  *****  ERROR  *****  IF  THIS  STATBCNT  IS  REACICD, 

—  T)€N  NOT  ENOUGH  BROKEN  AC  WERE  FOUND 
WRITE!  6, 9002)  WRKAC,  I  ACVC<  1 ) ,  NSELEC 

—  *£XIT  -  ALL  BROKEN  AC  HAVE  BEEN  SELECTED 

600  CONTINUE 

—  HJPDATE  LAST  WORD  OF  INPUT  VECTOR 

IACVCIIWORD)  »  XOR<IACVC(IWORD),IBRKVC) 

C— 
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C—  «UPOATE  NUMBER  OF  AIRCRAFT  IN  WVT  VECTOR 

IACVC(l)  =  IACVC(l)  -  NBRKAC 

C — 

C—  *£ND  IF  (ZERO  BREAKS  TEST) 

700  CONTINUE 
C — 

RETURN 

9001  FORMAT  COWttm*  MCDIST  ERROR  -  SEQUENTIAL  SAMPLING  ERROR*./. 

4  *  ««$»(  PBRKSEQ(l.NUC).  RDRAM  =  ‘.2F10.4) 

9002  FORMAT ( *0$$$<$$M  MCDIST  ERROR  -  INCONSISTENT  BROKEN  AIRCRAFT*, 

4  *  tttttttt  NBRKAC.  IACVC(i).  NSELEC  =  *.315) 

END 
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. . . 

SUBROUTINE  WCPR08<NWC,PBRKWC,PBRKSEQ.INDXWC,PWCPR0D) 

CHIIHHIHHHHIHmWtHHIIIIIlHHH»HIIIHI»HHHHH»HH» 


C++  WCPROB  -  INITIALIZE  WORK-CENTER  SEQUENTIAL  BREAK  PROBABILITIES. 
Cm  THIS  ROUTINE  CALCULATES  THE  PROBABILITIES  NECESSARY  TO 
C++*  SIMULATE  THE  DISTRIBUTION  OF  AIRCRAFT  BREAKS  INTO  THE  VARIOUS 
Cm  WORKCENTERS. 

Cm 


Cm  INPUTS  - 
Cm  NUC 

Cm  PBRKHC(J) 

Cm 

C++* 

Cm 

Cm  OUTPUTS  - 


-  NUMBER  OF  WORKCENTERS  BEING  MODELED. 

-  PROBABILITY  THAT  AN  AIRCRAFT  WILL  BREAK  INTO 
WORKCENTER-J.  NOTE  THAT  THIS  BREAK  MAY  BE  DUE 
TO  SORTIE  BREAKS  OR  OROUMJ-ABORTS.  DEPENDING 
ON  HOW  THIS  ROUTIfE  IS  CALLED. 


C++* 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 

Cm 


PBRKSEQU.JI-  PROBABILITY  THAT  AN  AIRCRAFT  BREAKS  INTO  THE 
WORKCENTER  INDICATED  BY  'INDXWCUI',  AND  DUES 
NOT  BREAK  INTO  ANY  OF  TIC  WORKCENTERS  - 
INDXWC(U+l),INDXWC(J+2),  ....  INDXWC(NWC) 

GIVEN  THAT  TIC  AIRCRAFT  MS  BROKEN  INTO  AT  LEAST 
0f£  OF  THE  WORKCENTERS  - 
INDXWC(J).  INDXWCtJ+l),  ....  INDXWC(NWC). 

THUS,  PBRKSEQU.(NHC)  MUST  EQUAL  1.0  . 

PBRKSEQ(2,J1-  PROBABILITY  THAT  AN  AIRCRAFT  HAS  BROKEN  INTO  THE 
WORKCENTER  INDICATED  BY  'INDXWCU)',  GIVEN  THAT 
TIC  AIRCRAFT  HAS  BROKEN  INTO  AT  LEAST  ONE  OF  THE 
WORKCENTERS  IMICATED  BY  - 
IWXUC(J),  INDXWCIJ+l),  ....  INDXWC(NWC). 

INDXWC(J)  -  INDICATES  THE  INDEX  OF  THE  WORKCENTER  WITH  THE 
JTH  LARGEST  BREAK  PROBABILITY.  THUS,  IM)XWC(l) 
INDICATES  TIC  WORKCENTER  WITH  TIC  LARGEST 
BREAK  PROBABILITY,  INDXWC(NWC)  IWIICATES  THE 
ONE  WITH  TIC  SMALLEST.  ETC.  . 

PWCPROD  -  PRODUCT-FORMULA  OVERALL  WORK-CENTER  BREAK-RATE. 

COfCUTED  FROM  TIC  INDIVIDUAL  WC  BREAK-RATES. 


Cmmmmmmmmmmmmmmmmmmmmmmm 

C— 


DIMENSION  PBRKWC(NWC),  PBRKSEQ(2,NWC),  IW)XWC(NHC) 


C — 

C—  +C0MPUTE  SORTED  ARRAY  OF  WORK-CENTER  INDICES  ACCORDING  TO 
C —  LARGEST -TO-SMALLEST  BREAK-RATE. 

CALL  SORTDSINWC+O.PBRKWC, INDXWC) 

C — 

C —  +INITIALIZE  END-POINT  PROBABILITIES 
PBRKSEQU.NUC)  =  1.0 
PBRKSEQ(2,NWC)  =  1.0 


C— 

C —  +D0  UNTIL (ALL  PROBABILITIES  HAVE  BEEN  CALCULATED) 

J  =  NWC  -  1 

POLO  =  1.0  -  PBRKWC< DOXWC(NWC) I 
100  CONTINUE 
C— 

C—  +C0WVTE  NEXT  PROBABILITIES 
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PROS  =  PBRKUC(IH)XWC(J)> 

PNEW  =  POLD  *  (1.0  -  PflOB) 

PBRKSEQ < 2 > J )  =  PROB/d.O  -  PNEW) 

PBRKSEQd.J)  =  PBRKSEQ<2» J)  *  FOLD 
POLD  =  PfCW 

—  *€ND  DO  (WC  LOOP) 

J  =  J  -  1 

lF(J.GT.O)  GO  TO  100 

—  *SAVE  PRODUCT-fORfUA  OVERALL  WORK-CENTER  BREAK-RATE 

PWCPROD  =  1.0  -  POLD 

C— 

RETURN 

END 


B-73 


SUBROUTINE  WCREADUFILE, MAXWC, NMC.PBRKWC, NCREWS, SRATE) 


C*+  WCREAD  -  READ  AND  INITIALIZE  WORK  CENTER  DATA. 

C***  WCREAD  READS  WORK-CENTER  DATA  FROM  THE  MAINTENANCE 
C***  HWEONER  IffUT  FILE.  THIS  DATA  IS  ASSUCD  TO  BE  ON  UNIT 
C"»  ‘IFILE".  ONE  FREE-FORMAT  RECORD  PER  WORK-CENTER. 

C***  THIS  ROUTINE  RETURNS  THE  NUMBER  OF  WORK-CENTERS  LOADED. *NWC‘; 

C**+  THE  BREAK-RATE  ARRAY.  ■PBRKWC1 :  TIE  SERVICE-RATE  ARRAY,  ‘SPATE* ; 
C***  AN©  THE  SERVERS  ARRAY,  “NCREWS*. 

C***  THE  SERVERS  ARRAY,  ‘NCREWS*,  REPRESENTS  THE  NUMBER  OF  CREWS 
C**+  AVAILABLE  PER  12-HOUR  SHIR.  THE  INPUT  FILE  CONTAINS  THE 
C***  TOTAL  NUMBER  OF  CREWS  AVAILABLE  FOR  TIE  PARTICULAR  WORK-CENTER. 
C**t  THE  12-HOUR  SHIFT  NUMBER  IS  COMPUTED  BY  DIVIDING  THIS  AVAILABLE- 
C***  SERVERS  IN  HALF  AND  ROUNDING  TO  THE  NEAREST  INTEGER  NUMBER. 

C***  IN  ADDITION.  IT  IS  ASSUMED  THERE  IS  ALWAYS  AT  LEAST  ONE  CREW 
C***  PER  SHIFT. 

C*»  THE  ffSC  DESCRIBING  THE  WORK-CENTER,  AND  THE  TOTAL  NUMBER  OF 
C***  SERVERS  ARE  ECHO-PRINTED,  BUT  ARE  NOT  SAVED  FOR  FUTURE  USE. 

C***  THE  MAXIMUM  NUMBER  OF  WORK-CENTERS  WHICH  CAN  BE  LOADED  IS 
C***  SPECIFIED  BY  TIC  ‘MAXWC*  INPUT  PARAMETER.  IF  MORE  THAN 
C**+  THIS  NUMBER  IS  READ,  AN  ERROR  ICSSAGE  IS  PRINTED,  AND  THE  SGM 
Cw  RUN  CONTINUES  WITH  ONLY  TIC  FIRST  ‘IWXWC‘  WORK-CENTERS. 

C***  INPUTS  - 

Cm*  IFILE  -  IWUT  FILE  NUMBER  FROM  WHICH  MAINTENANCE  MAICOWER 

C***  IWVT  DATA  IS  TO  BE  READ 

C***  MAXWC  -  MAXIMUM  NUMBER  OF  WORK-CENTERS  WHICH  CAN  BE  LOADED 

C***  OUTPUTS  - 

C***  NUC  -  NUMBER  OF  WORK-CENTERS  LOADED 

C***  PBRKWC  -  ARRAY  OF  WORK-CENTER  BREAK  RATES 

C**#  NCREWS  -  ARRAY  OF  WORK-CENTER  CREW  NUMBERS 

C***  SRATE  -  ARRAY  OF  WORK-CENTER  SERVICE  RATES 


C— 

OIIENSION  PBRKWC (MAXWC I,  NCREWS (MAXWC),  SRATE (MAXWC) 
CHARACTER  CAFSC*5 

C— 

C —  •PRINT  HEADER  FOR  ECHO-CICCK  OF  INPUT  DATA 
WRITE(6,9001) 

C — 

C —  ^INITIALIZE  NUMBER  OF  WORK  CENTERS 
NWC=1 

C — 

C —  *DO  UNTILINO  MORE  WORK  CENTERS  TO  LOAD) 

100  CONTINUE 
C — 

C—  *READ  fCXT  WORK-CENTER  RECORD 
READ( IFILE, 9000, END=200) 

k  CAFSC, PBRKWC (NWC), SERVERS, SRATE (NWC) 

C — 

C —  iPERFORM  ERROR  CHECK  ON  IHFUT  DATA 

IF( (PBRKWC (NWC). GT. 0.0). AND. (PBRKWC(NWC).LE. 1.0) ) 
Jr  GO  TO  150 
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IF (SERVERS. GE. 0. 0 )G0  TO  150 
IF  i'SRATE'.NWO.&E.  0.0)00  TO  150 

WR I TE ( 6 , 0004 ) CAPS C . PBRKWC I NWC ) . SERVERS , SPATE i NWC ) 
60  TO  100 
150  CONTINUE 

C —  tCOnPUTE  INTEGER-NUMBER  OF  CREWS  PER  12-HGUR  'SHIFT 

NCREWSI NWC )=MAXO(l,INTi SERVERS*. 5  *  .5)) 

C—  *ECHO-PRINT  WORK -'CENTER  INFORMATION 

WR I  TE  i  o,  ^>002 )  NWC ,  CAFSC ,  FBRKUC  i  NWC ) .  SERVERS ,  SPATE  ( NWC ) 

C —  ♦ INCREMENT  WORK -CENTER  INDEX 

NWC  -NWC  *  1 

L 

l —  *END  10  'WORK CENTER  LOOP) 

IF’NWC.LE.MAXWCiGO  TO  100 

*FRINT  ERROR  MESSAGE  IF  STILL  MORE  DATA  ON  THE  FILE 
S£AO(  IFILE.  *000, END=200 ) PBRKWC i NWC ) ,  SERVERS. '3RATE ( NWC 
WRITER.  «003>MAXWC 

:oo  continue 

C—  ♦hDJJST  NUMBER  OF  WORK-CENTERS  TO  ACCOCNT  FOR  EOF  READ 
NWC  -NWL " I 

C—  *CLOSE-OUT  MANPOWER  INPUT  FILE 
CALl  FCLOSEI IFILE' 


RETURN 

v»X)  FORMAT 1 1X.A5.  ll.F*.4,  l  x.FR.2«  IX-F'5. 4) 

*001  FORMAT' *1*./,’  .'X.  '**♦*♦**♦**♦♦♦♦*♦****+♦*+♦♦*♦♦+*♦***♦+**♦«♦♦". 

V  .7X, ’♦♦♦♦♦♦♦♦*♦  AIRCRAFT  MAINTENANCE  *«*+******j, 

i  ’X.  '**♦*-***♦♦♦*♦*♦*♦**♦**♦***-**♦*-****♦«*****♦♦♦", 

V  21X.  *  BREAK.".  IX.  *  TOTAL". 3X,  •SERVICE  RATE*./. 

4  ’X.'WC  r'.2X."  AFSC",  3X,  *  RATE* •  2X - '  SERVERS ” .  2X , * ( ACFT /HOUR )".//) 
*002  FORMAT ( 7 X . 1 3 . 3 X . A5 . 3 X . F&. 4 , 1 X . F7 . 2 , 3X , F* . 4  > 

<W3  FORMAT  rO»K«S«*t  WCREAD  ERROR  -  TCO  MANY  WORK  ENTERS  ". 
t  "IN  THE  INPI.IT  FILE"./. 

V  "  (M«HU  ONLY  THE  FIRST  M3,*  WORK -CENTERS  WERE  USED*,'./. 

4  '  mmti  INCREASE  -MAXWC-  PARAMETER  IF  YOU  WANT  MORE  WC-S*) 

’••04  FORMAT WNCREAD  ERROR  -  INVALID  WORK  CENTER  DATA*,/, 

V  '  PBRKWC.  SERVERS.  SPATE  -  *.JF8.3) 


END 


REtf.  FUNCTION  XNORM  (XMEAN.STDEV.SEED) 


C++  XNORM  -  KWH  RANDOM  SMPLE  FROM  A  NORMAL  DISTRIBUTION. 
Cw*  THIS  IS  A  REAL-VALUED  FORTRAN  FUNCTION  HHICH  G£)£RATES 
Cm  A  RANDOM  SAMPLE  ACCORDING  TO  A  NORMAL  PROBABILITY 
Cm  WITH  THE  GIVEN  INPUT  fEM  MD  STANDARD  DEVIATION. 

Cm  THE  TECHNIQUE  IS  TO  APPROXIMATE  A  NORMAL  DISTRIBUTION  USING 
Cm  THE  CENTRAL  LIMIT  THEOREM.  12  INDEPENDENT  SAULES  ARE 
Cm  DRANN  FROM  A  UNIFORMS.  1)  DISTRIBUTION  AND  THEN  ADDED. 

Cm  THE  RESULT  IS  APPROXIMATELY  NORMALLY  DISTRIBUTED  HITH  t*AN  b 
Cm  AND  STANDARD  DEVIATION  1.  THE  SAMPLE  IS  THEN  TRANSLATED  TO 
Cm  OBTAIN  A  SAMPLE  FROM  A  DISTRIBUTION  HITH  THE  GIVEN 
Cm  INPUT  fCAN  AND  STMDMD  DEVIATION. 

C++* 

Cm  INPUTS  - 

Cm  XMEAN  -  MEAN  OF  THE  NORMAL  DISTRIBUTION  FROM  WHICH 
Cm  THE  SAMPLE  IS  TO  BE  GENERATED. 

Cm  STDEV  -  STANDARD  DEVIATION  OF  NORMAL  DISTRIBUTION  FROM 
Cm  HHICH  SAMPLE  IS  TO  BE  GOtRATED.  IF  THIS  VALLE 

Cm  IS  NEGATIVE,  M  ERROR  MESSAGE  IS  PRINTED. 

Cm  INPUT /OUTPUT  - 

Cm  SEED  -  CURRENT  SEED  OF  RANDOM  NUMBER  GENERATOR. 

Cm  OUTPUT  - 

Cm  XNGRH  -  RANDOM  SAMPLE  FROM  A  NORMAL  DISTRIBUTION  HITH 
Cm  GIVEN  f€M  AND  STANDARD  DEVIATION. 


C — 

C—  *IF(STANDARD  LCVIATION  IS  LEGITIMATE) THEN 
IF(STDEV,LT.0.0)G0  TO  100 

C — 

C—  +0RAU  SAMPLES  FROM  12  UNIFORM  (0,1)  DISTRIBUTIONS  AND 
C —  ADJUST  THE  f€AN  OF  THIS  RANDOM  SAIfLE  TO  ZERO 

XN0RM4JNIFN1(SEED)+U«IFM1(SEED)+UNIFM1(S£ED)+)JNIFM1(SEED) 

&  HWIFM1  (SEED)+UN1FM1  (SEED)+UN1FH1  IS£ED)+UNIFM1  (SEED) 

S,  HWIFH1 (SEED) +UNIFH1 (SEED  >+UNlFHl(SEED)+UNIFMl(SEED)-6. 

C — 

C —  iCONVERT  TO  A  SAMPLE  FROM  DISTRIBUTION  HITH  APPROPRIATE 
C —  MEAN  AND  STANDARD  DEVIATION. 

XNORM  =  THEM  +  STDEV+XNORM 

C— 

C —  tELSE  (NEGATIVE  STANDARD  DEVIATION) 

00  TO  200 
100  CONTINUE 
C— 

C —  *SET  RETURN  VALUE  TO  ZERO  AND  PRINT  ERROR  MESSAGE 

XNORM  =  0.0 
WRITE (6. 9001 1STDEV 

C — 

C —  »£ND  rp  (STD  DEV  TEST) 

200  CONTINUE 
C — 

RETURN 
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9001  FCRfttTI "OW*****  XNORH  ERROR  -  fCGATIVE  STANDARD  DEVIATIONS/, 
&  •  $$*«<*««  STDEV  =  “iF10.5) 

END 
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C«HHiiimiHmHiiHininmmiHHHHimnmimimnmm 

SUBROUTINE  ZBITSL (MONESi IARRAY) 

ClHWHUHIHHIIWWUmmitHHHmUHHUHWtmOHlHIl 

C++  ZBITSL  -  ZERQ-OUT  1-BITS  IN  LEFTMOST  PORTION  OF  A  WORD. 

C++*  ZBITSL  IS  A  FORTRAN  SUBROUTINE  WHICH  WIU.  ZERO-OUT 
C»**  A  SPECIFIED  NUMBER  OF  1-BITS  IN  THE  LEFTMOST  PORTION  OF 
Cm  A  BIT-VECTOR.  THIS  BIT-VECTOR  IS  KEPT  IN  AN  ARRAY, 

Cm  ORGANIZED  IN  TIC  FOLLOWING  FASHION  -  THE  FIRST  WORD  OF 
Cm  THE  ARRAY  CONTAINS  THE  CURRENT  NUMBER  OF  1-BITS  IN  THE 
Cm  BIT-VECTOR,  AND  TIC  REMAINING  WORDS  OF  THE  ARRAY  CONTAIN 
Cm  TIC  ACTUAL  BIT-VECTOR.  THIS  ROUTINE  ZEROES  OUT  THE  PROPER 
Cm  1-BITS,  AND  UPDATES  THE  1-BIT  COUNTER  IN  TIC  FIRST  WORD 
Cm  OF  THE  ARRAY. 

Cm 

Cm  I  (CUT  - 

Cm  NONES  -  NUMBER  OF  1-BITS  TO  BE  ZEROED-OUT.  NOTE  THAT 
Cm  NONES  MUST  BE  .LE.  NUMBER  OF  IS  IN  THE  BIT- 

Cm  VECTOR. 

Cm  IICUT/OUTPUT  - 

Cm  [ARRAY  -  ARRAY  CONTAINING  THE  BIT-VECTOR  TO  BE 
Cm  MODIFIED.  I  ARRAY  U)  IS  A  COUNTER  WHICH  INDICATES 

Cm  TIC  CURRENT  NUMBER  OF  1-BITS  IN  TIC  BIT-VECTOR. 

Cm  THE  ACTUAL  BIT-VECTOR  IS  THE  CONSECUTIVE  BITS 

Cm  CGNTAIICD  IN  THE  WORDS  IARRAY (2)- IARRAY (LENGTH) 

Cm  COniON  I  (CUT  - 

Cm  LENGTH  -  LENGTH  (IN  COMPUTER  WORDS)  OF  THE  ARRAYS 
Cm  CONTAINING  TIC  VARIOUS  BIT-VECTORS;  I  FLY  VC,  ETC 

Cmm«mmmmmmmmmmmmmmmmmmmm+t 

C— 

PARAMETER  (WXAC=108,  MAXBIT=34,  HAXVEC=2+(I¥HAC-1)/HAXBIT 
C0CI10N  /ACSTATE/  LENGTH,  NACVCUWXVEC),  IFLYVC(NAXVEC), 

&  MAINVC!HAXVEC>,NORSVC(NAXVEC),  LOSTVC(MAXVEC) 

DIMENSION  IARRAY! 1) 

C— 

C —  # INITIALIZE  DO 

HEFT  =  NONES 

INDEX  =  2 

C— 

C —  *D0  WHIL£(ALL  APPROPRIATE  WORDS  HAVE  NOT  BEEN  MODIFIED) 

1000  CONTINUE 

IF(HEFT.LE.O)  GO  TO  WOO 
IFIIWEX.GT.LENGTH)  GO  TO  4000 

C — 

C —  COUNT  NUMBER  OF  IS  IN  NEXT  WORD 

NXT1S  =  N1BITS!  IARRAY(IWEX)  ) 

C— 

C —  *IF (NOT  «_L  1-BITS  IN  THIS  WORD  SHOULD  BE  ZEROED) 

IF(NXTlS.L£.ft£FT)  GO  TO  2000 

C— 

C —  +ZERO-OUT  THE  APPROPRIATE  NUWER  OF  IS 

IARRAYY INDEX )  =  XOR(  IARRAY! INDEX)  . 

I  LB1TS! IARRAY (IH€X), HEFT)  ) 

C— 
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♦ELSE  (Ml  1-BITS  IN  THIS  WORD  ARE  TO  BE  ZEROED) 
00  TO  3000 
CONTINUE 


C — 

2000 

—  ♦ZERO-OUT  THE  ENTIRE  WORD 

I ARRAY (INDEX)  =  0 

—  ♦END  IF  (ALL  IS  TEST) 

3000  CONTINUE 

—  ♦UPDATE  NUMBER  OF  IS  LEFT  TO  ZERO-OUT 

NLEFT  =  NLEFT  -  NXT1S 

—  ♦IMCREWXT  INDEX  FOR  NEXT  WORD  OF  BIT-VECTC* 

INEEX  =  INDEX  +  1 

♦OC  DO  (WORD  LOOP) 

GO  TO  1000 
(000  CONTINUE 

—  ♦UPDATE  IS  COUNTER  FOR  THIS  BIT  VECTOR 

I ARRAY! 1 )  =  I ARRAY! 1 )  -  NONES 

—  •PERFORM  ERROR  CHECX  TO  ENSURE  APPROPRIATE  NUMBER 

—  OF  1-BITS  WAS  ZEROED-OUT 

IF(HEFT.GT.O)  WRlTE(6,9001)NClNES,  I  ARRAY  (1).  NLEFT 

RETURN 

9001  FORMAT! ‘omttttt  ZBITSL  ERROR  -  NOT  ENOUGH  IS  TO  ZERO’,/, 
8.  ’*«***«  NOICS.l  ARRAY!  1),  NLEFT  =  *,315) 

END 
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